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SECTION  I 


INTRODUCTION 

For  the  past  five  years,  Harris  SAI,  Inc.,  Ann  Arbor, 

Michigan,  has  been  developing  advanced  computer-aided  simulations 
of  crossed-f ield  amplifiers.  Since  March  1,  1977,  Harris  SAI  has 
been  supported  in  this  effort  by  the  Air  Force  Office  of  Scientific 
Research  under  Contract  No.  F49620-77-C-0091.  This  Final  Technical 
Report  covers  the  entire  period  of  the  contract:  March  1,  1977 
through  February  20,  1980. 

The  report  is  in  two  parts.  Part  I  (this  volume)  describes 
the  work  on  the  distributed-emission  crossed-field  amplifier. 

Part  II  describes  the  parallel  effort  on  the  in jected-beam 
crossed-field  amplifier. 

A.  Research  Objectives 

Crossed-field  amplifiers  are  commonly  used  to  produce 

microwave  power  in  transportable  systems  where  light  weight, 

compactness,  high  efficiency,  and  peak  power  of  kilowatts  or 
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megawatts  are  required.  '  Radar  and  electronic  countermeasure 
systems  are  examples.  In  existing  distributed-emission  CFA's,  it 
is  desirable  particularly  to  raise  the  gain  and  efficiency  and 
extend  the  bandwidth  over  which  the  tube  will  both  start  and 
operate  in  the  desired  mode. 

The  present  research  program  involves  the  development  of 
advanced  computer-aided  design  techniques  and  experimental 
verification  of  the  theory.  This  work  will  lead  to  a  better 
understanding  of  crossed-field  interaction  mechanisms  and  provide 
a  guide  for  design  improvements.  The  computer  will  then  be 
available  to  provide  valuable  design  information  with  a  minimum 
of  costly  experimentation. 


Harris  SAI  has  studied  the  distributed-emission  crossed-f ield 
amplifier  (DECFA)  and  injected-beam  crossed-f ield  amplifier 
( IBCFA)  in  parallel  efforts.  The  primary  objectives  of  the  DECFA 
work  were  as  follows. 

(1)  To  develop  a  computer  simulation  of  the  distributed- 
emission  crossed-f ield  amplifier  applicable  to  both  forward-  and 
backward-wave  devices. 

(2)  To  verify  the  program  using  experimental  data  for 
production  tubes  under  normal  operating  conditions. 

(3)  To  test  the  capability  of  the  model  for  predicting  the 
operational  limits  of  a  tube,  such  as  the  bandwidth  and  the 
upper  and  lower  current  mode  boundaries. 

(4)  To  investigate  transient  behavior  and  possible  models 
for  the  modulator  to  improve  understanding  of  starting  conditions. 

(5)  To  incorporate  the  option  of  an  asymmetrical  cathode. 

B .  Status  of  the  Research  Effort 

A  computer  model  has  been  developed  for  both  forward-wave  and 
backward-wave  DECFAs.  This  program  is  believed  to  be  the  first 
to  simulate  successfully  a  backward-wave  CFA.  This  cylindrical 
model  of  the  reentrant  CFA  cross-section  includes  the  vane 
RF  circuit,  a  secondary-emitting  cathode,  and  an  external  modulator 
The  modulator  has  proven  to  be  important  in  stabilizing  the  model. 
Earlier  work  showed  a  "runaway"  effect  of  current  and  RF  power, 
a  problem  that  has  now  been  overcome.  An  improved  cathode  emission 
model  appears  to  reproduce  well  the  observed  phenomena  of  backbom- 
bardment  and  emission  limitation  of  the  tube  current. 

Experimental  data  have  been  obtained  for  production  tubes, 
the  backward-wave  Raytheon  QKS1842  and  QKS1705,  the  forward-wave 
Varian  SFD-261  and  Raytheon  QKS1319,  and  the  cathode-driven  QKS1840 
The  measurements  on  the  QKS1319  were  made  by  Raytheon  Company  under 
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a  subcontract  to  Harris  SAI.  Simulations  of  the  QKS1842,  SFD-261 
and  QKS1319  tubes  have  been  performed  at  their  midband  frequencies. 
Generally  there  is  reasonable  agreement  with  measurement, 
although  the  power  balance  is  not  yet  satisfactory  for  the  forward- 
wave  tubes.  The  anode  current  and  RF  output  power  are  now 
typically  7  percent  and  41  percent,  respectively,  higher  than 
measured  in  the  forward-wave  tubes  and  up  to  25  percent  below 
the  measured  values  in  the  backward-wave  tube.  However,  it  is 
believed  that  only  minor  adjustments  to  the  RF  network  model  and 
the  trajectory  calculations  are  now  needed  to  improve  these 
quantitative  results. 

It  is  noteworthy  that  the  multiple  operating  modes  and  high 
cathode  backbombardment  power  actually  observed  in  the  QKS1842 
are  predicted  by  the  simulation.  The  observed  oscillation  of  the 
QKS1319  at  current  levels  above  30A  is  also  correctly  predicted. 

In  late  1979  the  program  was  installed  on  the  Harris  550 
computer  which  made  large  blocks  of  CPU  time  economically  available. 
The  program  was  tested  with  up  to  27,900  rods  and  the  results 
showed  consistency  with  those  of  the  usual  7,000-rod  program. 

It  was  concluded  that  higher  resolution  than  provided  by  the 
7,000-rod  model  is  not  necessary.  The  off-center  cathode 
(as  in  the  production  SFD-261  tube)  will  be  incorporated  in  a 
future  model  by  making  a  perturbation  to  the  direct  radial  electric 
field.  Time  restrictions  also  precluded  a  study  of  transient 
effects,  but  the  present  program  has  the  capability  for  simulating 
starting  and  turn-off  with  a  modulator  of  either  line  or  hard- 
tube  type . 

This  study  has  succeeded  in  establishing  the  basic  modeling 
techniques  for  the  beam,  the  cathode,  the  RF  network  in  forward- 
wave  or  backward-wave  tubes,  and  the  modulator.  Only  minor 
refinements,  chiefly  to  the  trajectory  calculations,  are  now 
needed  to  provide  a  program  of  use  to  the  engineer  for  design 
purposes . 


C.  Summary 

Existing  DECFA  programs  and  their  limitations  are  compared 
in  Section  XI.  Section  III  and  Appendices  A-G  give  the  details  of 
the  tube  data  that  Harris  SAI  has  compiled  during  this  study,  and 
that  are  needed  to  perform  a  simulation.  The  major  assumptions  and 
approximations  of  the  model  are  listed  in  Section  IV,  and  Sections 
V  through  IX  provide  a  description  of  the  model.  The  computations 
involve,  principally: 

•  primary  and  secondary  emission  at  the  cathode  (Section  V) ; 

•  particle  trajectories  (Secton  VI); 

•  space-charge  fields  (Section  VII) ; 

•  the  RF  network  (Section  VIII); 

•  the  modulator  (Section  IX) ; 

Details  of  the  circuit  equations  are  given  in  Appendices  H-J. 
Section  X  gives  a  more  detailed  description  of  the  numerical 
procedure.  The  results  of  the  simulations  of  three  tubes  are 
presented  in  Section  XI,  which  includes  a  brief  discussion  of  the 
design  of  the  QKS1842.  The  computer  requirements  are  explained 
in  Section  XII.  A  typical  run  now  requires  about  20  CPU  hours 
on  the  Harris  550  minicomputer,  but  more  efficient  trajectory 
calculations  are  planned  which  may  halve  this  time.  A  brief 
explanation  of  the  program  output  information  is  given  in  Appendix 
K,  and  Appendix  L  reproduces  the  printed  results  from  a  typical 
run. 


The  professional  personnel  who  assisted  in  this  work  are 
acknowledged  in  Appendix  M. 


SECTION  II 

STATE  OF  THE  ART  OF  SIMULATION  OF 
THE  DISTRIBUTED-EMISSION  CFA 


For  over  twenty  years  various  workers  have  attempted  to  model 
the  distributed-emission  crossed-f ield  amplifier  with  only 
partial  success.  Like  the  magnetron,  this  tube  (Figure  1)  presents 
modeling  difficulties  because  of  its  interdependent  parts  (cathode, 
reentrant  beam,  RF  circuit,  modulator)  all  of  which  must  be 
treated  simultaneously. 
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Rigid-beam  models  of  CFA's  and  magnetrons  '  provide  useful 

qualitative  pictures  but  limited  quantitative  design  information. 

Of  the  two  types  of  computer  simulation  now  in  use,  the  rectangular 

single-wavelength  model^-17  treats  only  forward-wave  tubes,  while 

the  present  Harris  SAI  model  treats  the  entire  tube  as  a  cylindrical 

interaction  region  and  can  treat  both  forward-wave  and  backward- 

wave  devices.  It  is  based  on  an  earlier  program  by  Dombrowski 
18 

and  Price,  but  the  cathode  model,  space-charge  field  calculations, 
trajectory  calculations  and  the  RF  network  all  have  been  exten¬ 
sively  developed  by  Harris  SAI  and  the  variable  anode  voltage 
from  the  modulator  has  been  added.  A  similar  cylindrical  model 

is  being  used  to  study  high-power  relativistic  magnetrons  at 

19 

the  Naval  Research  Laboratory. 

Although  the  cylindrical  coordinate  system  requires  a 
more  complex  trajectory  calculation  than  that  for  rectangular 
geometry,  a  clear  advantage  is  that  the  measured  magnetic  field 
can  be  used  directly.  The  rectangular  model  requires  an  empirical 
reduction  of  the  magnetic  field  to  lower  the  Hartree  voltage. ^ 
Significant  advantages  of  the  Harris  SAI  cylindrical  model  are 
(1)  the  RF  circuit  model  which  can  treat  both  forward-wave 
and  backward-wave  networks  with  equal  ease,  and  (2)  the  realistic 
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Figure  1.  Continuous-cathode  distributed-emission 
crossed-f ield  amplifier.  (After  Skowron) 


angular  boundary  condition  which  represents  the  reentrant  beam 
as  a  continuous  electron  stream  as  in  the  actual  tube.  Both 
features  are  essential  in  simulating  a  backward-wave  amplifier 
correctly . ^ 

In  the  rectangular  model  study  performed  by  Harris  SAI  for 
12  13 

Lincoln  Laboratory,  '  backward-wave  interaction  was  simulated 
(a)  by  moving  the  reference  frame  in  the  direction  opposite  to 
the  beam  from  RF  input  to  RF  output,  and  (b)  by  moving  the 
reference  frame  with  the  beam  from  RF  output  to  RF  input. 

It  was  shown  that  method  (a)  gives  the  beam  dynamics  incorrectly 
because  the  beam  modulation  at  the  input  is  computed  immediately 
following  the  sever  region,  whereas  the  actual  beam  has  the 
modulation  produced  by  the  RF  interaction  between  the  output  and 
the  input.  Hence  the  computed  current  and  RF  output  power  are 
too  low.  Method  (b)  requires  a  prior  estimate  of  both  the 
RF  output  power  and  the  "hot"  phase  delay  of  the  RF  voltage. 

Trials  showed  that  the  behavior  of  the  beam  is  sensitive  to  both 
of  these  quantities.  The  value  of  RF  drive  power  derived 
varied  by  a  factor  of  2  or  more,  with  corresponding  phase 
variations,  in  successive  passes  around  the  tube.  In  addition, 
that  older  model  needed  double  the  interaction  impedance  relative 
to  the  measured  value  to  obtain  sufficient  RF  interaction. 

In  the  present  cylindrical  model,  the  network  equations  are 
identical  for  both  forward-wave  and  backward-wave  circuits:  only 
the  signs  of  the  voltages  on  even-numbered  vanes  are  changed  to 
enable  the  beam  to  "see"  a  backward-traveling  wave.  Since  the 
entire  cylindrical  region  is  viewed  at  every  time  step,  the  particles 
move  in  the  realistic  transient  RF  fields  from  all  the  vanes. 

The  only  difference  between  the  forward-wave  and  backward-wave 
calculations  is  the  larger  number  of  RF  periods  required  for  the 
latter  to  reach  saturation  while  RF  current  induced  by  the  beam 
at  the  RF  input  (where  the  beam  modulation  is  greatest)  is  fed 
back  along  the  RF  circuit  to  the  RF  output.  Transient  calculations 
with  ideal  spokes  of  charge  demonstrate  this  effect. 


The  features  of  the  cylindrical  and  rectangular  DECFA  models 
are  compared  in  Table  1. 

The  RF  network  is  a  unique  feature  of  the  Harris  SAI  CFA 
analysis.  While  the  NRL  magnetron  program  is  relativistically 
correct  and  uses  a  full  solution  of  Maxwell's  equations,  the  simple 
vane  boundaries  used  do  not  reproduce  the  measured  phase  delay 
or  interaction  impedance  of  a  CFA.  Typical  CFA  circuits  such 
as  the  backward-wave  strapped  vane  line  and  interdigital  line  or 
the  forward-wave  stub-supported  helix  and  meander  line  require 
a  three-dimensional  analysis  to  give  their  total  electromagnetic 
field.  However,  the  present  quasistatic  equivalent  network 
is  fully  adequate  for  a  two-dimensional  model. 

Although  all  previous  simulations  have  assumed  a  fixed 
anode  voltage  (the  "short-circuited  diode"  approximation) ,  the 
nonzero  impedance  of  the  modulator  allows  the  CFA  voltage  to 
fluctuate  with  the  beam  current.  These  fluctuations  are  too 
rapid  or  too  small  to  be  observed  experimentally,  but  are  essential 
in  the  simulation.  Both  the  Harris  SAI  CFA  model  and  the  NRL 
magnetron  model  now  impose  a  load-line  constraint  that  limits 
the  anode  current  by  reducing  the  voltage  as  the  current  rises. 

This  has  proved  necessary  to  ensure  a  stable  RF  solution. 

In  the  "cathode-driven  CFA"  a  gain  of  up  to  30  dB  has  been 

obtained  by  forming  the  secondary-emission  cathode  into  a  slow- 
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wave  structure,  on  which  the  RF  drive  signal  is  applied.  (See 

Appendix  F.)  Litton  Industries  have  developed  a  rectangular 

15-17 

model  of  this  device.  However,  that  program  gives  limited 

resolution.  The  model  uses  a  space  charge  array  having  26  axial 
points  and  11  radial  points,  and  allows  only  192  particles  per 
RF  wavelength.  It  also,  more  significantly,  neglects  recirculated 
charge  at  the  RF  input  and  ignores  backward  power  flow  on  the 
RF  circuits.  The  realistic  direct  field  is  included  from  the 
cathode  vanes,  but  only  a  single  space  harmonic  of  the  RF  field 
is  present  in  the  Litton  program. 


TABLE  1 


COMPARISON  OF  DISTRIBUTED- EMISSION 
CROSSED-FIELD  AMPLIFIER  MODELS 


I 
» 

i 

|  Cylindrical  Model 

1.  Model  interaction  region 
covers  entire  cylindrical 
reentrant  tube. 

2.  Interaction  region  is 
stationary. 

3.  Reentrancy  of  the  tube 
supplies  the  boundary 
condition  at  the  ends  of 
the  interaction  region. 

The  mode  of  beam  inter¬ 
action  (number  of  spokes 
of  charge  around  the  tube) 
develops  automatically: 
no  one  mode  is  selected  by 
the  user. 

4.  The  trajectory  equations 
use  cylindrical  polar 
coordinates. 

5.  The  measured  magnetic 
field  is  supplied 
directly  to  the  model. 


Rectangular  Model 

1.  A  rectangular  interaction 
region  covers  a  single 
retarded  RF  wavelength 
only. 

2.  The  interaction  region 
moves  around  the  tube  at 
the  cold  RF  phase  velocity 

3.  The  beam  distribution  is 
forced  to  be  periodic  over 
the  interaction  region. 

Only  a  single  beam  mode  is 
treated  with  the  periodicity 
of  the  user-selected  inter¬ 
acting  space  harmonic. 


4 .  The  curvature  of  the  anode 
is  neglected. 

5.  The  magnetic  field  used  in 
the  model  is  artificially 
set  lower  than  the 
measured  value  to  maintain 
the  ratio  of  the  anode 
voltage  to  the  Hartree 
voltage. 


TABLE  1  (cont.) 


The  RF  fields  are  given 
by  a  Green's  function 
for  the  actual  anode  vane 
structure.  Network  circuit 
equations  give  the  transient 
solution  for  the  vane 
voltages.  All  space 
harmonics  and  all  time 
harmonics  are  included 
automatically. 


The  solution  for  a 
backward-wave  network  is 
obtained  from  that  for  a 
forward-wave  network 
just  by  changing  the 
signs  of  the  voltage 
and  current  on  alternate 
vanes. 


At  present,  10,000  rods 
of  charge  and  a  257 
(angular)  x  33  (radial) 
space-charge  array  cover 
15  RF  wavelengths  around 
the  tube. 


Only  a  single  interacting 
space  harmonic  is  included 
at  the  drive  frequency. 

The  vane  circuit  is 
replaced  by  an  equivalent 
transmission  line.  The 
time  distribution  of  the 
induced  current  is 
approximated  by  the 
spatial  distribution  over 
the  interaction  wavelength. 

For  backward-wave  interaction 
the  reference  frame  must 
move  with  the  beam.  The  RF 
output  power  is  supplied 
and  the  drive  power  derived. 
The  signs  of  the  impedance 
and  attenuation  in  the  trans¬ 
mission  line  equations  are 
opposite  to  the  forward-wave 
signs. 

About  4,000  rods  of  charge 
and  a  99  x  51  space-charge 
array  cover  one  RF  wave¬ 
length.  The  resolution 
of  the  rectangular  model 
is  greater  than  in  the 
cylindrical  model  because 
of  the  larger  array. 

Cost  considerations, 
not  inherent  model  limita¬ 
tions,  are  the  reason  for  the 
difference . 


[•] 


TABLE  1  (cont.) 


The  secondary-emission  9. 

coefficient  can  be  spec¬ 
ified  as  a  function  of 
primary  impact  energy. 

Up  to  120  emission  sites 
cover  the  circular  cathode 
(over  15  to  20  RF  wave¬ 
lengths)  .  Options  of 
time  smoothing  of  the 
secondary  charge  and 
space-charge  limitations  are 
included  to  reduce  the 
fluctuations  at  each 
emission  site. 

The  anode-cathode  voltage  10. 

is  variable  in  time  and 
controlled  by  an  external 
modulator  circuit  which 
reduces  the  voltage  in 
response  to  an  increase 
of  current. 

The  present  program  11. 

requires  about  260,000 

32-bit  words  of  CPU 

memory.  One  pass  around 

the  tube  (15  RF  periods) 

with  7,000  particles 

requires  about  13,500  CPU 

seconds  on  an  IBM  360/67 

computer. 


There  are  96  secondary- 
emission  sites  per  RF 
wavelength  in  the  Varian 
model,  where  time  smoothing 
is  performed  in  the  moving 
reference  frame,  and  Child's 
Law  limiting  is  included. 
The  Litton  model  reemits 
every  incident  primary 
rod  with  a  new  charge 
obtained  from  the  secondary- 
emission  coefficient. 


The  anode-cathode  voltage 
is  constant. 


The  Harris  SAI  version  of 
the  rectangular  model  uses 
an  estimated  60,000  32-bit 
words  of  CPU  memory.  One 
pass  around  the  tube  requires 
approximately  1,200  CPU 
seconds  on  an  IBM  360/67 
computer. 


In  order  that  the  computer  program  be  useful  for  design  it 
is  important  to  correlate  the  theoretical  results  with  experimental 
measurements  for  a  range  of  conditions  in  different  tubes.  Fairly 
good  agreement  (anode  current  about  20  percent  below  the  measured 
value  and  RF  output  power  about  15  percent  too  low)  has  been 
reported  by  Varian  Associates^  for  a  space-charge-limited 
cathode  concentric  with  the  anode.  The  Varian  model/  however, 
shows  variations  of  +  8  percent  of  the  mean  anode  current  and 
RF  power  over  successive  passes  around  the  tube.  That  model 
has  also  been  tested  for  the  production  SFD-261  CFA  in  which  the 
cathode  is  noncylindrical,  and  on  a  tube  with  a  tapered  vane 
pitch.  The  computed  power  grows  less  rapidly  along  the  circuit 
than  observed  in  the  instrumented  tubes.  For  the  Litton  program 
there  are  no  published  comparisons  with  measurement. 

The  primary  aim  of  these  simulations  is  to  predict  the 
operating  anode  current  and  RF  output  power.  However,  the  tube 
designer  also  needs  to  know  the  range  of  conditions  under  which 
the  CFA  can  start  and  operate  without  oscillation  or  interfering 
modes,  and  under  which  the  cathode  will  supply  sufficient  current. 

Starting  conditions  have  been  studied  using  ballistic 
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theory  '  '  and  with  space  charge  but  the  steady- state 

assumption  of  the  latter  model  is  probably  not  realistic.  The 

14 

space-charge  forces  modify  the  trajectories  significantly, 
making  the  full  simulation  necessary.  The  Harris  SAI  program 
can  include  both  the  formation  of  the  beam  and  the  buildup  of 
voltage  due  to  a  line- type  modulator,  on  an  accelerated  time  scale. 

Useful  tests  of  the  cathode-emission  model  are  its  ability 
to  predict  the  backbombardment  power,  to  form  a  stable  steady- 
state  charge  hub,  to  predict  the  backbombardment  power  with  RF 
drive,  and  to  demonstrate  the  limits  of  available  cathode  current 
as  the  anode  current  is  increased.  Both  Harris  SAI  and  Varian 
originally  used  an  unsmoothed  emission  model  in  which  each  primary 
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charge  generated  secondary  charge  computed  from  the  impact 
energy.  However,  subsequent  anomalous  results  have  shown  the 
need  for  some  form  of  smoothing  to  suppress  excess  charge, 
since  the  space-charge  mesh  cannot  resolve  the  local  potential 
minimum.  The  Varian  model  uses  time  smoothing  and  applies  a 
simple  Child’s  Law  emission  limit  using  the  first  mesh  interval 
as  an  elemental  diode.  It  has  not  succeeded  in  reproducing 
emission  limits.  Harris  SAI  uses  an  improved  model,  originally 
developed  for  a  thermionic  cathode  in  the  NRL  magnetron  study. 

Here  both  the  local  electric  field  and  the  charge  near  the  cathode 
limit  the  emission.  The  present  application  to  secondary  emission 
is  new.  Initial  results  described  in  this  report  are  qualitatively 
well  supported  by  observations  of  backbombardment  power  and 
emission  limitation. 

Mode  boundaries  have  not  been  predicted  in  any  published 
analysis.  As  the  first  step,  the  Harris  SAI  program  now  shows 
the  threshold  RF  output  power  for  oscillation  due  to  reverse- 
directed  power  flow  on  the  circuit.  It  also  gives  a  reasonable 
estimate  of  the  RF  drive  power  required  to  suppress  this  oscillation 


SECTION  III 


TUBE  DATA 

A  considerable  amount  of  experimental  tube  data  has  been 
obtained  by  Harris  SAI  during  the  present  and  earlier  studies. 

Table  2  compares  the  six  tubes  for  which  data  have  been  obtained; 
those  data  are  listed  in  Appendices  A-F.  The  secondary-emission 
characteristics  of  the  four  cathode  materials  are  given  in 
Appendix  G. 

The  following  quantities  must  be  supplied  by  the  engineer 
for  performing  a  simulation: 

•  Circuit  type  (forward-wave  or  backward-wave) ,  and 
number  of  active  vanes; 

•  Magnetic  field; 

•  Tube  dimensions,  i.e.  width  in  the  magnetic-field  direction, 
anode  and  cathode  radii,  shape  of  cathode  if  eccentric  or 
not  cylindrical,  pitch  of  vane  circuit,  vane  spacing  on 

the  anode; 

•  Frequency  and  operating  band; 

•  Operating  voltage; 

•  RF  drive  power; 

•  Circuit  cold-test  measurements  of  phase  delay  per  pitch, 
interaction  impedance  at  the  vane  tips  and  attenuation 
per  pitch; 

•  Secondary  emission  of  cathode  as  a  function  of  primary 
impact  energy; 

•  Estimate  of  operating  current; 

•  Type  of  modulator  (line  or  hard-tube) . 

These  quantities  are  listed  on  the  forms  used  for  Appendices 

B-F. 
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SECTION  IV 


ASSUMPTIONS  AND  APPROXIMATIONS 

The  following  are  the  main  assumptions  and  approximations 
of  the  model. 

1.  The  magnetic  field  is  axial,  uniform,  and  constant  in 

time . 

2.  Motion  in  the  magnetic-field  direction  is  zero;  this  is 
a  two-dimensional  model  only,  using  rods  of  charge  to  represent 
the  beam. 

3.  The  electric  field  is  treated  as  an  average  value  over  the 
width  of  the  tube  in  the  magnetic-field  direction.  (Actual  beams 
show  maximum  anode  bombardment  and  sole  emission  in  the  center  of 
the  transverse  plane.) 

4.  The  beam  velocity  and  fundamental  circuit  phase  velocity 
are  much  less  than  the  velocity  of  light  (v2/c2<<  1)  such  that 
motion  is  nonrelativistic ,  forces  due  to  RF  magnetic  fields  are 
negligible  and  Poisson's  equation  replaces  the  wave  equation  for 
the  electric  field. (The  self-magnetic  field  due  to  the  rotating 
beam  can  be  included  as  a  small  correction  in  future  simulations.) 

5.  For  computing  the  electric  field,  the  electrode  surfaces 
are  treated  as  smooth  and  perfectly  conducting.  However,  the  RF 
network  accounts  for  measured  skin  losses. 

6.  The  cathode  is  a  cylinder  concentric  with  the  anode. 

An  "off-center  cathode"  will  be  treated  in  future  versions  of 
the  program. 

7.  The  slow-wave  structure  is  replaced  by  a  low-frequency 
RF  network.  Each  node  of  this  network  represents  one  of  the 
active  vanes  of  the  CFA. 


8.  The  total  voltage  seen  by  the  beam  at  any  vane  is  the 
sum  of  the  modulator  voltage  and  an  RF  voltage.  The  modulator 
voltage  is  assumed  to  vary  more  slowly  in  time  than  the  RF  voltage, 
in  order  that  the  RF  network  equations  can  be  solved  independently 
of  the  modulator. 

9.  At  any  instant  the  average  RF  voltage  over  all  the 
vanes  is  much  less  than  the  modulator  voltage,  so  that  the 
modulator  circuit  can  be  treated  independently  of  the  RF  network. 

10.  The  cold  RF  network  is  perfectly  matched  at  each  end. 

11.  The  anode  vane  segments  are  identical  and  equally 
spaced.  However,  a  "tapered-pitch"  circuit  will  be  included 
in  future  versions  of  the  program. 

12.  Between  the  anode  vanes  the  electric  field  is  uniform 
around  the  circumference.  (This  field  would  be  insufficiently 
accurate  for  a  cathode  circuit  but  is  a  good  approximation  where 
electrons  move  rapidly.) 

13.  There  are  no  neutral  or  ionized  gas  atoms  present. 

14.  The  electron-electron  collision  time  is  large  compared 

with  the  transit  time  through  the  amplifier.  Nearest  neighbor 

electron  forces  are  small  compared  with  the  effects  of  the  average 

field  due  to  many  electrons,  and  the  collisionless  particle- in-cell 
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model  is  applicable  to  treat  collective  space-charge  interaction. 

15.  The  relative  change  of  the  local  electric  field  during 
one  time  step  of  the  calculations  is  small.  The  effects  of  this 
approximation  have  been  checked  by  running  with  a  reduced  time 
step  and  verifying  that  the  results  are  consistent. 

16.  The  computed  results  are  insensitive  to  the  number  of 
simulation  rods  used.  This  approximation  has  been  tested 
separately  for  each  tube  under  simulation.  About  7,000  rods 
suffice. 


17.  Electrons,  primary  or  secondary,  are  emitted  over 
several  sites  covering  the  cathode.  The  total  emission  current 
is  uniform  over  each  site.  Sufficient  sites  should  be  chosen 
to  give  at  least  six  per  spoke  of  charge  (per  RF  wavelength 
around  the  circumference) .  Usually  120  sites  suffice. 

18.  If  sufficient  charge  is  available  at  any  of  the 
cathode  emission  sites  to  produce  a  local  potential  minimum, 
the  potential  minimum  is  treated  as  coincident  with  the 
cathode  and  its  depth  is  neglected. 

19.  The  time  interval  between  successive  electron  emissions 
is  an  integer  multiple  of  the  time  step. 

20.  The  electron  thermal  velocities  are  negligible  in 
the  RF  interaction;  however,  the  user  may  select  a  single 
initial  energy  for  primary  or  secondary  electrons  emitted  from 
the  sole. 

21.  The  modulator  acts  as  a  constant  voltage  source  in 
series  with  a  resistance,  an  inductance,  and  the  parallel 
combination  of  a  capcitance  and  the  CFA. 

22.  The  modulator  inductance  and  capacitance  are  chosen 
so  that  transient  circuit  voltages  and  currents  decay  over  a 
few  RF  periods,  in  order  to  accelerate  the  development  of  a 
steady  state  in  the  simulation. 


SECTION  V 


CATHODE  EMISSION 


A.  Introduction 

The  operation  of  any  CFA  depends  critically  on  the  ability 
of  the  cathode  to  supply  sufficient  current.  A  realistic 
model  enables  the  tube  designer  to  predict  the  emission  limits 
and  backbombardment  heating  with  thermionic  or  secondary- 
emitting  surfaces  of  various  types.  The  macroscopic  model 
used  here  and  in  the  injected-beam  CFA  study  (Part  II  of  this 
report)  appears  to  reproduce  the  observed  behavior  well. 

The  user  may  simulate  either  thermionic  or  secondary  emission 
or  both.  For  thermionic  emission  a  fixed  number  of  charge 
rods  is  made  available  for  emission  at  each  time  step.  These 
rods  are  distributed  uniformly  around  the  cathode.  However, 
the  actual  number  emitted  is  controlled  by  the  local  electric 
field  and  charge  density  in  a  manner  described  below.  These 
primary  charges  may  be  used  throughout  the  run.  Alternatively 
they  can  be  used  for  the  first  NTINJ  time  steps  only,  to  develop 
a  beam  hub  before  the  secondary-emission  model  takes  over. 

The  thermionic  charges  may  be  combined  with  or  replaced 
by  secondary  charges  produced  by  backbombarding  electrons.  The 
switch  NSECSW  controls  the  inclusion  or  exclusion  of  secondary 
emission. 

B.  Secondary  Emission 

A  full  description  of  the  secondary-emission  process  as 
modeled  by  Harris  SAI  is  presented  here. 

The  computations  of  impact  energy,  backbombardment  current, 
and  the  summation  of  secondary  charge  are  performed  in  the 
"predictor"  section  of  the  trajectory  subroutines  (Section  VI) . 
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A  separate  subroutine  allocates  the  secondary  charge  to  new 
simulation  particles  (rods)  at  the  start  of  the  following  time 
s  tep . 

When  a  particle  is  computed  to  intercept  the  sole  (cathode) , 
its  last  time  step  is  repeated  in  several  subintervals  (typically 
1/15  the  normal  time  step) .  The  angular  electric  field  is 
linearly  interpolated,  as  a  function  of  radial  distance  from 
the  sole,  between  a  zero  value  at  the  sole  and  the  value  at 
the  start  of  the  time  step.  The  radial  space -charge  and  RF 
fields  are  not  changed  for  the  subintervals,  but  the  exact 
direct  radial  field  between  the  cylindrical  anode  and  sole  is 
evaluated  for  each  subinterval.  The  reduced  time  interval 
gives  a  more  accurate  sole  bombardment  energy. 

The  angular  region  of  the  sole  where  the  interception 
occurred  is  determined  as  the  position  for  secondary  emission. 

The  cylindrical  surface  may  be  divided  into  one  to  256  such 
regions . 

According  to  a  user  option,  secondary  emission  may  be  set 
to  zero  and  the  primary  charge  treated  as  if  it  had  not  been 
emitted  in  the  case  where  its  radial  position  at  the  start  of 
the  time  step  was  less  than  1.0002  times  the  sole  radius. 

That  procedure  suppresses  the  emission  of  secondary  charge  due 
to  primary  particles  emitted  and  collected  in  the  same  time 
step.  The  amount  of  such  primary  charge  is  recorded  and  is 
printed  at  each  step. 

The  secondary-emission  coefficient  is  supplied  as  a 
tabulated  function  of  primary  impact  energy,  for  incidence 
normal  to  the  cathode  surface.  For  each  intercepted  particle, 
the  program  interpolates  the  table  linearly  using  the  computed 
impact  energy. 

The  angle  of  primary  incidence,  3,  is  used  to  correct  the 

?  4 

secondary-emission  coefficient"  to  a  value  Sz  using  the 
equation: 


S]eP(i-cose) 
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(1) 


where  6^  is  the  value  determined  for  normal  incidence,  and  the 
input  parameter  p  is  typically  0.55.  An  incident  charge  q 
then  results  in  the  production  of  secondary  charge  60q  at 
the  same  angular  site.  As  all  the  trajectories  are  followed  in 
a  given  time  step,  the  total  secondary  charge  to  be  emitted 
is  summed  for  each  of  the  emission  sites  around  the  cathode. 

Next,  program  options  of  time  smoothing  and  spatial 
smoothing  of  the  charge  to  be  emitted  may  be  invoked  for  each 
emission  site.  Firstly,  the  time-smoothing  option  (subroutine 
TSMSEC)  combines  the  charge  with  that  from  previous  steps  in 
the  manner  described  by  McDowell‘S  with  a  smoothing  parameter, 
TIMSEC,  between  zero  and  unity.  A  value  TIMSEC  of  0.15  causes 
about  73  percent  of  the  secondary  charge  from  a  given  step  to 
be  emitted  over  the  eight  succeeding  steps.  Secondly,  the 
space-smoothing  option  (subroutine  SMSEC)  distributes  the  charge 
determined  for  a  given  emission  site  over  seven  adjacent  sites 
in  the  proportions 

S:S^:S  :S  :S  :S  :S  , 

where  S  is  a  user-chosen  parameter.  This  option  has  not  been 
used  in  major  runs. 

At  this  point  the  total  thermionic  or  secondary  charge 
available  for  emission  in  each  cathode  site  in  the  next  time 
step  is  known.  Next  the  local  electric  field  and  the  charge 
already  at  the  cathode  are  used  to  control  the  actual  emission. 

C.  Control  of  Emitted  Charge 

The  cathode  is  divided  into  between  120  and  256  angular 
regions  for  computation  of  the  emitted  charge.  It  is  assumed 
that  in  each  region  at  a  given  time  step  one  of  two  conditions 
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exists : 


(a)  Space-charge  limitation,  defined  here  as  the  condition 
in  which  the  charge  available  for  emission  during  the  time  step 
is  at  least  sufficient  to  create  a  potential  minimum  within 

the  first  mesh  interval  normal  to  the  surface. 

(b)  Emission  limitation,  in  which  the  normal  electric 
field  is  directed  towards  the  cathode  throughout  the  first 
mesh  interval. 

In  case  (b)  all  the  available  charge  is  emitted;  in  case 
(a)  sufficient  charge  is  emitted  to  reduce  the  local  electric 
field  to  zero  on  the  cathode  surface. 

Hence  the  instantaneous  local  electric  field  normal  to  the 
cathode,  including  direct,  RF  and  space-charge  fields,  controls 
the  emission,  while  the  magnetic  field  is  ignored.  When  a  region 
of  the  cathode  is  space-charge  limited  (with  the  definition  above) 
the  implicit  assumption  is  that  a  sheath  of  charge  covers  the 
surface  there  and  creates  a  "virtual  cathode"  at  the  potential 
minimum. 

The  time-dependent  feature  of  this  model  allows  the  emitted 
charge  to  fluctuate  with  the  instantaneous  RF  and  space-charge 
f ields . 

D.  Approximations  of  the  Emission  Model 

Three  approximations  made  here  are; 

(1)  the  time  required  to  form  the  potential  minimum 

(of  the  order  of  the  local  plasma  oscillation  period)  is  much 
less  than  the  time  step  of  the  simulation, 

(2)  the  distance  from  the  cathode  to  the  potential  minimum 
Srs  much  less  than  the  mesh  interval  normal  to  the  surface  and 

(3)  the  depth  of  the  potential  minimum  is  much  less  than 
the  potential  difference  over  the  first  mesh  interval. 

Strictly,  approximation  (1)  is  violated  in  the  model,  where 


a  time  step  of,  typically,  1/10  the  cyclotron  period  is  also 
1/10  of  the  plasma  oscillation  period  in  a  Brillouin  stream. 
However,  the  error  is  probably  unimportant  for  small  fractional 
changes  per  step  in  the  emitted  charge. 


To  justify  approximations  (2)  and  (3)  consider  the  plane- 
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parallel  diode  with  no  magnetic  field.  Langmuir's  solution 


shows  that  the  distance  of  the  potential  minimum  from  the 


cathode  is  approximately  equal  to  the  Debye  shielding  distance 
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Ad.  This  distance  decreases  as  the  local  charge  density  increases 


and  varies  as  the  mean  normal  velocity  of  electron  emission 
from  the  cathode. 


The  depth  Vm  of  the  potential  minimum  is  given  by 


m 


=  V  £n 


(: 


net  emitted  current 


available  secondary  current 


3 


Volts 


(2) 


where  V  is  the  mean  energy  of  available  secondaries  about  2 
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electron  volts.  Published  estimates  of  the  potential 

minimum  and  its  position  in  a  crossed-f ield  diode  have  the 

same  order  of  magnitude.  An  upper  estimate  of  XD  can  be 

obtained  using  the  Brillouin  charge  density  in  the  applied 

magnetic  field  and  an  initial  secondary  electron  energy  of 

2eV.  In  the  tubes  studied  here,  these  estimates  range 

between  1/103  and  1/163  times  the  cathode-anode  spacing. 

The  radial  mesh  interval  used  in  the  model  is  1/32  times  the 

cathode-anode  spacing.  For  a  saturated  emitter  current  of  ten 

times  the  net  emitted  current  (with  an  effective  secondary- 

emission  coefficient  of  1.11)  the  potential  minimum  is  only 

-4.6  V,  while  the  anode  voltage  is  10  kV  or  more. 
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It  should  be  mentioned  in  passing  that  Fischer  has  made 
a  quite  different  approximation,  namely  that  the  potential 
minimum  equals  the  mean  ballistic  impact  energy  (eV)  due  to 
the  RF  field,  i.e.  over  100  volts.  However,  with  the  resulting 
high  space-charge  field  the  ballistic  equations  are  probably 
not  valid.  Moreover  the  steady  potential  assumed  by  Fischer  is 


inconsistent  with  the  time-varying  RF  field  at  each  cathode  point 

E.  Algorithm  for  Control  of  Emission 

The  steps  followed  by  the  program  are 

(1)  calculation  of  space-charge  and  RF  electric  fields 
at  each  emission  site  (subroutines  ESCCAT  and  RFC AT) , 

(2)  counting  of  the  charge  already  present  above  each 
site  (subroutine  CCGCAT) , 

(3)  limiting  of  the  available  charge  where  necessary 
(subroutines  SECCAT  and  NOEMIT) , 

(4)  printing  of  diagnostic  information  (subroutine  PRTCAT) , 

and 

(5)  initialization  of  new  particles  at  the  cathode 
surface  (subroutine  INSECS) . 

For  solution  of  Poisson's  equation  the  anode-cathode 
region  is  covered  by  a  computational  mesh  of  chosen  radial 
interval  Ar  and  angular  mesh  interval  A0,  which  may  be 
equal  to  or  less  than  the  angular  interval  DTHSEC  of  each 
emission  site.  Figure  2  shows  the  configuration  at  the  cathode 
with  L':  =  0.5  DTHSEC  in  this  example. 

The  region  of  angle  DTHSEC  between  rc  and  rc  +  Ar,  where 
rc  is  the  cathode  radius,  is  used  to  count  the  charge  already 
present.  Applying  Gauss's  Law  and  setting  the  radial  electric 
field  E  to  zero  at  r  =  rc  gives 

(charge  allowed  for  emission) 

+  (charge  already  in  the  region) 


c  _  E  h  Ar  DTHSEC  (3 


Figure  2.  Mesh  elements  at  the  cathode  in  the 
DECFA  simulation. 


where  £q  is  the  permittivity  of  free  space,  Er  is  the  total 
radial  electric  field  evaluated  at  rc  +  Arand  h  is  the  width 
of  the  cathode  in  the  magnetic-field  direction. 

The  RF  electric  field  is  computed  at  the  midpoint  of  the 
angular  segment  DTHSEC .  The  space-charge  fields  are  already 
stored  at  the  Poisson  mesh  points.  Their  values  are  averaged 
over  the  mesh  points  at  radius  Ar  above  each  emission  site 
using  an  area-weighting  procedure. 

The  charge  in  the  region  rQ  £  r  £  rc+1.5Ar  is  assigned 
to  the  angular  emission  segments  by  the  same  area-weighting 
procedure.  A  fraction  RFRAC  of  each  particle  is  considered 
to  lie  within  a  distance  Ar  of  the  cathode,  given  by 

RFRAC  =  1  -  J  r  -  ~  j Ar  (4) 

for  a  particle  at  radius  r. 

If  the  charge  allowed  for  emission  exceeds  the  available 
charge  stored  for  emission  at  any  site,  that  site  is  treated 
as  emission  limited  and  only  the  available  charge  is  emitted. 

F.  Test  Results 
1 .  Child's  Law 

Here  the  magnetic  field  is  set  to  zero  and  a  thermionic 

emission  current  is  chosen  sufficient  to  give  space-charge- 

limited  operation.  With  the  present  emission  limitation  and 

the  predictor-corrector  trajectory  calculation  the  anode 

currents  from  the  model  are  two  percent  greater  than  those 
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predicted  by  Child's  Law  for  the  cylindrical  geometry. 

They  follow  the  three-halves  power  law  as  the  voltage  is  varied. 


2 .  Secondary  Emission 

Before  the  present  cathode  model  was  implemented,  runs  with 
no  limit  imposed  on  secondary  emission  accumulated  too  many 


charge  rods  and  showed  excess  backbombardment  power  at  the 
cathode.  Apparently  the  model  had  insufficient  resolution  to 
return  excess  charge  to  the  cathode, at  which  the  radial  mesh 
interval  is  larger  than  the  distance  to  the  actual  potential 

9 

minimum.  Similar  instabilities  have  been  computed  by  Varian 
without  smoothing. 

Four  runs  were  performed  to  test  the  effects  of  varying 
degrees  of  smoothing  at  the  cathode.  The  results  of  Figures  3 
and  4  and  Table  3  show  the  computed  anode  current  and  RF 
output  power  in  forward-wave  CFA  as  the  simulation  progresses 
over  50  RF  periods.  These  results  were  obtained  without  the 
modulator  and  show  a  runaway  effect  of  current  and  RF  power. 
Notice  that  the  computed  anode  current  and  output  power  are 
higher  than  the  measured  22  A  and  1.65  PQ,  but  decrease  as 
the  degree  of  smoothing  in  the  model  increases.  The  thermionic- 
emission  model  gives  the  closest  agreement,  for  the  anode 
current  is  within  10  percent  of  the  measured  value. 

All  the  secondary-emission  runs  (runs  5,  6,  7  and  8) 
used  a  Child's  Law  limiting  procedure,  which  produces  similar 
effects  to  the  more  refined  emission  limitation  described 
above.  Runs  5,  6,  and  7  used  120  equiangular  regions  subdividing 
the  cathode  for  computing  the  local  secondary  emission.  In 
run  5,  no  additional  smoothing  was  performed.  Runs  6  and  7  used 
a  time-smoothing  method,  following  that  of  Varian,  but  in  a 
fixed,  not  a  moving,  reference  frame.  Complete  spatial 
smoothing  was  applied  in  run  8  by  distributing  the  total 
secondary  charge  evenly  around  the  cathode  in  each  step.  The 
"thermionic  cathode"  of  run  10  was  simulated  by  placing  300 
primary  charge  rods  at  equal  intervals  around  the  circumference 
at  each  time  step  and  suppressing  secondary  emission.  The 
corresponding  total  emitted  current  of  204  A  compares  with  the 
177  A  used  by  Varian  and  is  enough  to  maintain  a  retarding 
average  electric  field  around  the  cathode  for  about  75  percent 
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TABLE 


of  the  time.  The  net  cathode  current  is  the  emission  current 
minus  the  current  of  electrons  returned  to  the  cathode. 

The  details  of  the  behavior  of  these  test  runs  are  not 

fully  understood.  Apparently  time  smoothing  significantly 

reduces  the  noise  in  the  discrete  particle  model.  The  angular 

distribution  of  secondary  charge  around  the  cathode  is  important 

14  21 

in  ballistic  theory  where  the  RF  field  produces  phase  focusing  ' 
but  the  corresponding  effects  under  space-charge-limited  conditions 
in  an  actual  tube  are  not  so  clearly  defined. 


SECTION  VI 


PARTICLE  TRAJECTORIES 

In  each  time  step,  the  radial  and  angular  components  of 
position  and  velocity  of  the  simulation  particles  must  be 
advanced,  with  the  vane  voltages  given  at  the  start  of  the 
step.  The  rapid  angular  variation  of  the  RF  electric  field 
and  the  cylindrical  geometry  present  special  problems  of 
resolution.  Initial  trials  at  Harris  SAI  used  a  simple  one-step 
trajectory  calculation  in  local  rectangular  coordinates. 

The  predictor-corrector  calculation  now  used  is  less  efficient 
but  more  accurate  than  the  simple  one-step  algorithm.  It 
includes  two  electric  field  computations  per  time  step. 

At  the  start  of  the  time  step,  at  time  t,  the  following 
quantities  are  known  for  each  particle: 

1.  Radius, 

2.  Angle, 

3.  Velocity, 

4.  Local  electric  field  (as  radial  and  angular 
components) ,  and 

5.  Charge. 

The  calculation  of  the  position  and  velocity  at  time  t  +  At 
proceeds  by  using  a  basic  algorithm  twice  for  each  particle. 

A.  Basic  Algorithm 

The  equation  of  motion  of  a  particle  of  charge  q  and  mass 
m  in  an  electric  field  E  and  magnetic  field  B  is 


ov  =  a  (E  +  V  x  B) 
at  m 


A  transformation  from  cylindrical  coordinates  (r,8)  to  a  local 
rectangular  system  is  performed.  At  the'  start  of  the  time 
step,  define  local  coordinates  as 


=  0 


y  =  r 
P 


and  the  velocity  as 


v  =  r9  , 
xp 


v  =  r  , 

yp 


with  the  local  electric  field  components 


E  =  E 
x  9 


=  E. 


Next,  the  particle  equation  of  motion  Equation  1  is 

solved  for  the  new  position  (x  ,y  )  and  velocity  (v  ,v  ). 

n14n  xn  Yn 

The  scheme  described  by  Vaughan^  is  used  to  give  cycloidal 
trajectories  which  would  be  exact  in  a  uniform  electric  and 
magnetic  field. 


Finally,  the  new  cylindrical  coordinates  are  computed 
using  the  equations 
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(8) 


and 


P 

n 


t®n"1,xn/yn)  +9p 


/ 


where 

veloci 


9 

P 

ty 


is  the  previous  angular  position, 
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(rn,  r ®n)  is  given  by 


and  the  new 
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and 


X  V 

r  *  =  -v  —  +  v  — 

n  n  yn  r„  xn  r 

J  n  n 


(11) 


a  rotation  of  the  coordinate  system  through  the  angle  tan 

(x  /y  )  . 
n  1  n 

The  direct  electric  field  is  radial  and  equal  to  -V  / 

(r»:n(r  /r  )  )  at  radius  r,  where  V„  is  the  anode-sole  voltage, 
as  a 

r,  is  the  anode  radius  and  r  is  the  sole  radius.  The  electric 

cl  S 

field  changes  with  position  both  because  of  the  cylindrical 
anode  and  sole  and  because  of  the  RF  field.  The  refinement 
of  this  calculation  accounts  for  this  change. 


B .  Predictor 

For  each  particle,  the  stages  of  the  calculations  are 


as  follows: 


1.  Store  the  initial  position  and  velocity  in  separate 
arrays  for  use  by  the  corrector. 

2.  Advance  the  trajectory  using  the  basic  cylindrical 
algorithm  and  the  field  E ^ ,  check  for  electrode  interception, 
and  store  the  predicted  positions  and  velocity  for  time  t  +  At 

3.  Use  the  predicted  position  to  evaluate  the  local 
charge  density  at  time  t  +  At. 

Solve  Poisson's  equation  for  the  space-charge  potential 
at  time  t  +  At,  derive  the  space-charge  field,  and  advance 
the  network  voltage  solution  from  time  t  to  time  t  +  At. 

At  this  stage  the  integrals  are  incremented  for  the 
Fourier  components  of  the  induced  currents  and  vane  voltages. 
These  components  are  used  for  output  information  only  at 
intervals  of  one  RF  period. 

C.  Corrector 

For  each  particle: 

1.  Compute  the  local  electric  field  E2  at  the  predicted 
position,  using  the  space-charge  field  and  the  resultant  RF 
field  from  all  the  vane  voltages. 

2.  Repeat  the  trajectory  advance  from  time  t  to  time 
t  +  At,  but  now  using  the  arithmetic  mean  of  the  electric 
fields  E-^  and  E^,  to  derive  the  corrected  position. 

3.  Recompute  the  electric  field  E^  at  the  corrected 
position,  increment  the  induced  current  on  each  vane  (at 
time  t  +  At) ,  and  store  the  local  electric  field  E^  for  use 
as  the  field  E^  of  the  next  time  step. 

The  time-consuming  element  of  this  procedure  is  the 
summation  of  the  RF  electric  field  on  all  particles  due  to 
the  separate  vanes  of  the  tube.  This  summation  is  performed 


twice  for  each  time  step.  However,  trials  have  shown  that  an 
accurate  resolution  of  the  spatially  varying  electric  field  is 
essential  for  realistic  RF  interaction  to  appear  in  the  overall 
solution. 


D.  Test  Results 

The  program  has  been  tested  using  ballistic  motion  of  a 
single  particle  in  the  direct  field  only  for  one  RF  period. 
Results  have  been  compared  as  the  time  step  is  successively 
reduced.  In  the  absence  of  RF  fields,  the  simple  "predictor- 
only"  calculation  gives  errors  of  4.04,  2.07,  0.839,  0.420, 
0.204  and  0.09  percent  in  total  energy  for  time  steps  of 
respectively  1/20,  1/40,  1/100,  1/200,  1/400,  and  1/800  the 
cyclotron  period,  whereas  including  the  corrector  reduces 
the  error  to  only  0.04  percent  for  a  step  of  1/20  the 
cyclotron  period. 

The  charge  distributions  of  Figures  5  and  6  have  been 
obtained  for  the  QKS1300  amplitron.  With  only  one  field 
evaluation  per  step,  prior  to  the  corrector  calculation,  little 
interaction  appears  (Figure  5) ,  but  with  the  local  field 
recalculated  at  the  end  of  the  step,  the  charge  spokes 
develop  (Figure  6) . 

When  the  time  step  is  halved,  the  induced  vane  currents 
and  sole  backbombardment  power  agree  to  within  a  few  percent 
and  the  plotted  beam  profiles  appear  consistent. 

E .  Proposed  Development 

The  present  procedure  has  two  disadvantages,  namely 

(1)  the  "corrected"  particle  position  is  consistent 
with  the  derived  RF  fields,  but  not  with  the  space-charge  field 
unless  a  second  Poisson  solution  is  performed  after  the 
corrector  step,  and 


4  4  4  4 
4  *  I 
4  4  4  4 
4  4  4 
4  *  *  I 
4  4  4  4 
4  *  ( 

*  * 

4  * 

*  I 

4  4 

4  4  I 
4  4  4 

4  4  4 

4  4  4 


4  I  *\4 
4  4  -4 
4  4  •— 


* 

X 

3 

X 

* 

X 

X 

X 

X 

1 

% 

X 

X 

X 

X 

—* 

X 

X 

1 

* 

X 

X 

* 

■X 

X 

X 

X 

X 

— 

* 

X 

X 

•v 

X 

1 

X 

X 

X 

X 

X 

X 

\ 

N  ^ 

X 

X 

X 

X 

t 

X 

X 

X 

X 

X 

X 

X 

X 

— « 

• 

X 

X 

X 

X 

•* 

X 

X 

X 

X 

* 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

! 

V* 

X 

> 

X 

X 

1 

X 

X 

■> 

X 

< 

X 

X 

X 

X 

X 

1 

X 

X 

•:j. 

X 

■  < 

X 

■r 

-> 

X 

X 

X 

> 

X 

X 

_n 

X 

X 

■X 

X 

—4 

• 

V- 

■> 

V 

•> 

■»<■ 

1 

X 

>r 

■** 

X 

X 

X 

X 

X 

■r 

•<r 

e  5.  Charge  distribution  computed  with  one  electric  field 
ation  per  time  step. 


(2)  the  RF  field  is  time-consuming  to  evaluate  twice 
per  step. 

In  further  development  of  this  program  it  is  planned  to 
test  an  alternative  time-centered  trajectory  calculation, 
which  may  prove  more  efficient  by  eliminating  one  of  the  RF-field 
calculations  per  time  step.  If  this  proves  successful  an 
alternative  RF  field  calculation  will  then  be  tested  for  a 
further  increase  of  efficiency. 

1 .  Time-Centered  Trajectory  Calculation 

The  steps  of  the  time-centered  trajectory  procedure  are  as 
follows  for  each  active  particle: 

(1)  Start  with  the  particle  velocity  v  at  time  t  -  it/2 
and  its  position  r  at  time  t. 

(2)  Evaluate  the  local  electric  field  E  (r,t),  separating 
the  space-charge  field  and  the  r  and  •  components  due  to  the 
voltage  on  each  vane.  Store  these  components  in  two  work 
arrays  of  dimension  NVANES. 

(3)  Advance  the  velocity  v  to  time  t +  it/2  as  in  the 
present  cylindrical  model. 

(4)  For  each  vane  .increment  the  induced  current  by  the 
amount  due  to  the  particle  (See  Appendix  H) .  The  stored  field 
components  at  (r,t)  are  used  here,  and  the  velocity  v(t)  is 
approximated  as 

v ( t )  =  \  (v(t  -  4~)  +  v  (t  +  ~  ) )  .  (12 

2  r  2 

The  work  arrays  are  then  free  to  be  reused  for  the  next 
particle . 

(5)  Advance  the  position  r(t)  to  r(t  -r  it)  and  test  for 
interception  on  the  anode  and  cathode. 


(6)  Store  the  charge  density  at  time  t  +  At  for  the 
next  solution  of  Poisson's  equation. 


2 .  RF  Field  Calculation 

Because  the  details  of  the  circuit  {which  is  generally 
three-dimensional)  are  modeled  with  an  equivalent  network, 
it  is  necessary  to  separate  the  space-charge  field  from  the 
field  due  to  the  vane  voltages  in  order  to  compute  the 
induced  current  on  each  vane.  The  form  of  network  used  for 
the  backward-wave  device  requires  also  that  the  RF  voltages 
be  defined  separately  from  the  direct  voltage  due  to  the 
modulator . 

At  present,  the  program  computes  the  RF  field  at  each 
particle  by  summing  the  components  from  all  the  vanes  within 
a  certain  angular  interval  centered  at  the  particle.  For 
each  vane  the  RF  voltage  is  multiplied  by  the  Green's 
functions  for  the  radial  and  angular  components.  The  induced 
current  is  computed  simultaneously.  This  procedure  requires 
only  two  41  x  41  RF  field  arrays,  but  is  time-consuming. 

1c  may  be  practicable  to  replace  the  direct  summation  of 
the  RF  field  on  each  particle  with  a  solution  of  Laplace ' s 
equation.  Then  the  induced  charge  and  current  on  each  vane 
will  be  obtained  from  the  normal  component  of  the  space-charge 
field  as  in  Appendix  H.  However,  this  modification  must  be 
carefully  tested  to  ensure  that  it  can  reproduce  the  RF  field 
with  sufficient  resolution  with  the  available  machine  storage. 
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SECTION  VII 


SPACE- CHARGE  FIELDS 

The  electrostatic  field  is  obtained  by  solving  Poisson's 
equation  in  finite-difference  form  at  each  time  step.  The 
axial  self-magnetic  field  is  now  neglected.  It  is  estimated 
below  as  less  than  one  percent  of  the  external  magnetic  field. 

A.  Solution  of  Poisson's  Equation 

In  this  two-dimensional  model  each  particle  is  treated 
as  a  rod  of  charge,  of  arbitrary  length  in  the  magnetic-field 
direction.  The  charge  density  and  electric  fields  are  functions 
of  the  cylindrical  coordinates  (r,f)  only. 

The  cylindrical  interaction  region  is  subdivided  with  a 
mesh  of  points  for  the  finite-difference  solution  of  Poisson's 
equation.  Charge  from  the  simulation  particles  is  assigned 

to  the  mesh  by  an  "area-weighting”  or  "cloud-in-cell"  procedure. 

2  8 

Poisson's  equation  is  solved  by  standard  numerical  methods 
for  the  cylindrical  geometry  with  the  anode  and  cathode  voltages 
specified.  Only  the  natural  reentrancy  boundary  condition  is 
imposed  in  the  angular  direction. 

The  radial  and  angular  electric  fields  are  computed  and 
stored  in  separate  arrays  prior  to  the  calculation  of  the 
trajectories.  The  area-weighting  method  is  again  used  to 
interpolate  the  space-charge  field  to  the  position  of  a  given 
electron.  For  greater  accuracy,  the-  space-charge  field  is 
dor  it  *d  for  zero  anode  v<.  1  c  age  md  the  exact  external  electric 
5 1 e  1  u  is  added  a:  the  parti  1  e  position. 

Ir.  or  J-  r  t  .1  mach  ire  o'  <.  r  f  low  during  the  Poisson 

so:  u*  ion,  the  vh  ■  an  ...n  e;  actually  *s  stored  in  normalized 
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(the  number  of  units  of  initial  charge 
assigned  to  the  mesh  point) 


x  (radius  at  the  mesh  point) / (cathode  radius) 

This  number  is  of  the  order  of  unity.  A  factor,  VSCCON, 
multiplies  the  normalized  potential  solution  to  give  the 
actual  ratio  of  voltage  to  external  magnetic  field.  (The 
magnetic  field  is  included  in  this  term  only  for  greater 
efficiency  in  the  trajectory  calculations.) 

Most  of  the  simulations  described  here  use  a  mesh  with 
256  angular  intervals  around  the  tube  (for  up  to  18  RF 
wavelengths,  or  18  charge  spokes) .  There  are  32  radial  intervals 
between  the  anode  and  the  cathode. 

A  useful  check  of  the  space-charge  calculation  is  obtained 
by  allowing  a  hub  of  charge  to  develop  with  uniform  cathode 
emission  and  without  the  RF  field,  the  solution  for  the  "smooth¬ 
bore  magnetron."  The  analysis  and  test  results  are  given  in 
Appendix  J. 

B .  Charge  Per  Rod 

The  initial  charge  per  rod  is  determined  from  the  input 
variable  NBRILL .  This  is  the  number  of  rods  required  to  fill 
the  cylindrical  interaction  region  with  a  Brillouin  charge 
density  of  Cgoc:/n,  where  ^ / 2 v  is  the  cyclotron  frequency  and 
n  the  electron  charge- to-mass  ratio.  Usually  NBRILL  is  taken 
as  10,000  or  20,000:  the  smaller  value  where  the  anode 
voltage  is  further  below  cutoff.  Section  XII  shows  consistent 
test  results  with  NBRILL  as  high  as  45,000. 


Because  of  the  sensitivity  of  the  anode  current  (computed 
or  measured)  to  the  magnetic  field  it  is  useful  to  estimate 
the  correction  that  would  be  added  by  including  the  axial 
field  due  to  the  rotating  beam.  The  estimates  of  Table  4  are 
obtained  from  Ampere  '  s  Law  as 
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self 


=  y A  I  /h 
0  circ 


(1) 


for  a  circulating  current  I  over  a  tube  width  h  in  the 

11  c  ire 


The  values  of  I  .  are  computed 
circ  ^ 


magnetic-field  direction 
as  in  Appendix  J.  The  self-magnetic  field  opposes  the  applied 
field,  but  the  correction  is  less  than  one  percent.  With  the 
approximation  of  Equation  1  it  can  readily  be  included  in  the 
computation . 
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THE  EQUIVALENT  NETWORK  FOR  THE  RF  CIRCUIT 

A.  Introduction 

The  slow-wave  vane  structure  of  the  DECFA  is  represented 
by  an  equivalent  low-frequency  network  for  which  electrostatic 
node  voltages  are  derived  at  each  time  step.  With  approximations 
discussed  with  the  detailed  equations  in  Appendix  H,  the  RF 
and  direct  fields  are  treated  independently  and  the  vane-to- 
cathode  capacitances  are  included  as  parts  of  the  network. 

The  network  capacitances,  inductances  and  resistances  are  computed 
from  the  measured  cold-test  values  of  phase  delay,  interaction 
impedance  and  attenuation  at  the  drive  frequency. 

B.  Comparison  with  Other  Models 

The  present  formulation  has  definite  advantages  over  the 
continuous-transmission-line  RF  interaction  model  used  by  previous 
workers  (e.g..  Refs.  11  and  16). 

1.  By  creating  the  total  field  due  to  all  the  vanes, 

the  model  automatically  includes  all  space  harmonics.  In  fact, 
the  spatial  Fourier  components  of  the  RF  field  have  angular 
periodicity,  so  the  space  harmonics  appropriate  to  an  infinite 
linear  structure  are  not  required  for  computing  the  total  RF 
field . 

2.  The  fully  transient  network  equations  automatically 
include  RF  signals  generated  at  frequencies  other  than  the 
drive  frequency.  The  behavior  of  these  signals  depends,  of 
course,  on  the  values  of  the  network  elements  and  the  resulting 
dispersion  and  impedance  characteristics.  Only  a  single-frequency 
RF  drive  signal  is  applied  at  present. 
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3.  Since  the  entire  cylindrical  tube  is  treated  through¬ 
out  the  calculation,  only  the  natural  reentrancy  boundary 
condition  is  imposed  on  the  charge  distribution.  Moreover,  for 
a  given  anode  voltage  and  magnetic  field,  the  charges  are  free 
to  move  into  any  spoke  configuration  or  interaction  mode. 

4.  Power  may  flow  in  either  direction  along  the  network 
and  reflected  waves  due  to  "hot"  mismatch  at  the  input  and 
output  are  automatically  included.  Hence  the  program  can 
predict  RF  oscillation. 

In  a  third  model,  a  full  time-varying  solution  of  Maxwell's 
equations  is  obtained  for  the  electric  and  magnetic  fields  in 
a  magnetron  using  the  actual  "toothed"  vane  shape  as  the  anode 
boundary.  No  cold-test  dispersion  or  impedance  data  are  used. 
However,  that  model  may  be  difficult  to  adapt  to  the  complex 
forward-wave  and  backward-wave  circuits  of  crossed-f ield 
amplifiers,  such  as  the  strapped  vane  line  of  the  QKS1842  and  the 
stub-supported  helix  of  the  SFD-261  tube. 

C .  Forms  of  Network 

The  general  form  of  the  low-frequency  equivalent  network 
is  shown  in  Figure  7.  Nodes  1,  2,  3,  4  ...,  represent 
the  vanes  of  the  CFA  circuit  and  their  voltages  produce  quasistatic 
RF  fields  seen  by  the  electrons.  Only  the  driving  current  is 
given  a  sinusoidal  variation  in  time;  all  other  currents  and 
voltages  may  have  arbitrary  time  dependence,  periodic  or  non¬ 
periodic.  Appendix  I  gives  the  general  equations  and  the  calcu¬ 
lations  used  for  the  network  elements. 

The  network  used  in  the  present  calculations  is  a  simple 

case  of  the  general  strucutre  of  Figure  7.  This  network,  shown 

for  the  forward-wave  tube  in  Figure  8,  uses  only  one  driving 

current  source,  and  the  admittance  Y „  between  alternate  nodes 

c 

is  zero.  The  three  elements  G  (conductance),  C  (capacitance) 

(3  S 

and  l/rfa (inductance)  are  chosen  automatically  in  the  program, 
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wave  circuit 


jure  8.  Simple  equivalent  network. 


N 


using  the  values  of  delay  per  section  (pitch) ,  characteristic 
2 

impedance  V^/(2  x  power  flow),  and  attenuation  (dB)  per  pitch 
for  a  given  frequency.  The  program  also  computes  the  matching 
load  elements  and  F^,  and  the  driving  current  1^  to  give 
a  particular  RF  drive  power. 

In  a  backward-wave  tube,  the  direction  of  beam  motion 
is  still  defined  as  from  vane  1  to  vane  NVANES ,  while  the  RF 
power  flow  and  measured  external  phase  delay  are  in  the 
reverse  direction.  The  form  of  network  used  is  the  same  as 
in  Figure  3,  but  three  modifications  are  made: 

1.  The  driving  current  is  applied  at  the  last  vane 
(number  NVANES  in  Figure  8) . 

2.  The  voltage  seen  by  the  network  at  the  even-numbered 
vanes  is  defined  as  (-1)  times  the  voltage  acting  on  the  beam. 

For  this  transformation  it  is  necessary  to  separate  the  RF 
vane  voltages  from  the  modulator  voltage. 

3.  The  current  induced  in  the  even  vanes,  as  seen  by  the 
network,  is  defined  as  (-1)  times  the  value  computed  from  the 
beam  motion. 

With  this  transformation  applied  to  the  instantaneous  voltages 
and  currents,  the  phase  delay  per  pitch  of  each  Fourier  component 
of  the  wave  seen  by  the  beam  equals  r  minus  the  delay  seen  by 
the  network,  and  the  directions  of  the  waves  seen  by  the  beam 
and  network  are  opposite,  as  required.  However,  the  network 
time  delay  remains  equal  to  the  externally  measured  value, 
and  the  transient  response  time  of  the  numerical  solution  is 
correct,  at  least  at  the  drive  frequency. 

The  simple  network  of  Figure  8  has  the  advantage  of  elimi¬ 
nating  spurious  reflected  signals.  These  appeared  during  trial 

comoutation  with  non-zero  values  of  the  admittance  Y  between 

c 

alternate  vanes.  Harris  SAI  has  attempted  to  use  the  general 
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network  of  Figure  7  to  fit  experimental  values  of  phase  delay 
or  impedance  at  more  than  one  frequency,  or  to  fit  a  given 
value  of  the  cutoff  frequency.  However,  this  has  not  proved 
straightforward,  as  there  is  in  general  no  physical  solution 
for  the  network  elements  for  an  arbitrary  dispersion  or 
impedance  variation.  At  present,  instead,  the  network  of 
Figure  8  is  used,  and  its  elements  are  changed  by  the  program 
as  the  drive  frequency  is  changed. 

D .  Oscillation  Conditions 

The  network  of  Figure  8  is  sufficiently  general  to  include 
the  backward  traveling  wave  due  to  reflection  at  the  ouput 
or  to  the  induced  currents  from  the  beam.  This  wave,  identical 
to  the  interacting  wave  apart  from  its  direction,  causes  feedback 
of  power  to  the  RF input.  If  the  resultant  power  flow  at  the 
first  node  is  in  the  backward  direction,  the  value  of  RF  drive 
power  printed  by  the  program  is  negative  (Appendix  I)  and  the 
network  is  producing  an  oscillation.  This  is  not  a  numerical 
effect:  the  network  voltages  are  stable  and  well-defined 

throughout  the  calculation.  This  oscillation  condition  occurs 
only  when  the  induced  RF  current  exceeds  a  threshold  value, 
determined  in  turn  by  the  RF  voltage  driving  the  charge  to  the 
anode  and  hence  by  the  tube  output  power.  Computed  and  measured 
threshold  powers  for  oscillation  are  compared  in  Section  XI,  with 
order-of-magnitude  agreement  for  two  tubes.  In  the  actual  CFA 
the  oscillation  shows  as  a  mode  boundary  and  may  induce  spurious 
output  signals  at  other  frequencies. 


SECTION  IX 


THE  MODULATOR 

A.  Line  Type  and  Hard  Tube  Modulators 

Most  practical  distributed-emission  crossed-f ield  amplifiers 

are  operated  in  pulses  each  covering  several  thousand  RF  periods. 

Two  types  of  modulator,  "line-type"  or  "hard-tube",  are  used  to 

supply  the  anode-cathode  voltage  and  the  direct  current  during 
i  29 

the  pulse.  ' 

In  the  SFD-261  and  QKS1319  forward-wave  CFAs  a  hard-tube 
modulator  supplies  an  open-circuit  voltage  between  the  pulses, 
and  the  tube  is  started  in  each  pulse  by  applying  the  RF  drive 
power.  The  anode-cathode  voltage  falls  to  the  operating  value 
as  the  anode  current  rises.  A  quench  electrode  removes  the 
charge  from  the  CFA  at  the  end  of  each  pulse.  In  the  QKS1842 
backward-wave  CFA  the  anode-cathode  voltage  is  pulsed  by 
discharging  a  line-type  modulator  while  the  RF  drive  power  is 
being  fed  through  the  circuit. 

One  purpose  of  the  present  simulations  is  to  reproduce  the 
anode  current  and  RF  output  power  under  the  dynamic  steady-state 
operating  conditions  during  a  single  pulse.  A  simple  equivalent 
circuit  is  used  to  represent  the  modulator. 

B .  Equivalent  Circuits 

In  Figure  9,  VgP£j.  is  a  constant  open-circuit  voltage,  L 
and  R  are  series  inductanceandresistar.ee,  and  C  is  the  CFA 
capacitance  together  with  external  leakage  capacitance.  The 
induced  current  ICFA(t)  is  obtained  from  the  charge  motion  in 
the  CFA  (Appendix  H) .  The  following  ordinary  differential 
equations  are  solved  for  the  voltage  V^p^(t)  seen  by  the 


electrons  and  for  the  current  I  supplied  by  the  modulator: 


d  V 


CFA 


dt 


JCCT  '  ICFA 


(1) 


d  I 


CCT 


dt 


VOPEN  ~  VCFA  ~  RICCT 


(2) 


Here  the  direct  voltage  V(-.pA  is  treated  as  independent 
of  the  RF  voltages  on  the  vanes  by  neglecting  the  sum  of  the 
RF  voltages  at  each  instant  {Appendix  H) . 


C .  Typical  Parameters 

For  a  hard-tube  modulator,  v0pen  typicaHy  about  1.1 


times  the  mean  operating  voltage  V^p^.  For  a  line-type  modulator 
it  is  about  twice  the  operating  voltage.  This  value  and  the 
circuit  elements  L,  R  and  C  are  supplied  to  the  present  program. 
In  practice  only  the  CFA  voltage  VCFA  and  current  ICCT  are 
accurately  measured,  as  average  values  over  the  steady-state 
portion  of  a  pulse. 


To  run  a  simulation  the  user  must  first  estimate  the 
desired  operating  voltage  and  current  I  (given  the 

magnetic  field) .  Next  the  open-circuit  voltage  vqPEN  is 
chosen  for  the  modulator  type.  The  resistance  R  is  then 
derived  for  a  steady  state  using  Equation  2  with  d  I^CT/dt  =  0. 
Table  5  shows  the  values  used  in  present  runs. 

Using  typical  tube  values  (e.g.  L  =  1.0x10  H,  R=  1421:.,  an 
_  1  o 

C  =  1.0x10  F  on  the  QKS1842)  give  characteristic  times, 

L/R  and  RC,  for  decay  of  transients  that  are  unacceptably  long 
for  the  computer  simulation.  (These  tines  are  about  1400  RF 
periods,  a  fraction  0.08  of  the  pulse  length  of  1.7  ;,s  in  the 
QKS1842  and  fourteen  times  the  length  of  a  typical  run.) 


Therefore  the  inductance  and  capacitance  must  be  chosen  arbitrarily 
by  the  user.  Present  simulations  use  zero  inductance  and  use 
for  C  the  capacitance  of  the  cylinders  that  form  the  anode  and 
cathode  of  the  CFA  (ignoring  the  detailed  structure  of  the 
vanes) . 

With  the  line-type  modulator  in  the  QKS1842  the  CFA 
voltage  starts  at  its  full  open-circuit  value  instead  of  at  zero 
as  in  the  actual  tube.  Hence  transient  values  of  the  current 
and  power  appear  on  the  computed  curves  for  the  first  50  RF 
periods.  The  alternative  procedure  is  to  choose  a  small  but 
nonzero  inductance,  and  to  start  V^FA  at  or  near  zero.  This  should 
give  an  accelerated  simulation  of  the  starting  phenomena. 


D .  Stabilizing  Effect  of  the  Load  Line 

Both  the  measured  and  the  computed  values  of  current  and 
voltage  in  Table  5  satisfy  the  constraint  of  Equation  2  with 
d  ICCT/dt  =  0.  This  constraint  serves  to  stabilize  the  calcu¬ 
lation.  In  earlier  simulations  without  the  modulator  the  current 
and  RF  power  could  rise  indefinitely  at  a  fixed  anode  voltage. 
However,  the  modulator  reacts  by  reducing  the  voltage  if  the 
current  exceeds  the  mean  operating  value. 

The  present  modulator  model  is  applicable  also  to  studies 
over  a  frequency  band,  when  the  modulator  setting  is  fixed  and 
only  the  drive  frequency  varies .  The  computed  current  and  voltage 
should  adjust  automatically  for  best  synchronism.  Of  course, 
the  correct  RF -circuit  phase  velocity  and  impedance  must  be 
provided  for  each  frequency  used. 
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DESCRIPTION  OF  THE  NUMERICAL  SOLUTION 


The  major  elements  of  the  calculation  are  shown  in  the 
flow  diagram.  Figure  10.  The  following  are  the  principal  steps 
in  advancing  the  particle  coordinates  and  the  RF  and  direct 
voltages  by  one  time  step,  At.  Auxiliary  computations  not 
listed  here  provide  output  information  such  as  the  Fourier 
components  of  the  vane  voltages  and  currents,  and  the  time- 
averaged  anode  current,  cathode  emission,  electron  bombardment 
power  on  the  electrodes,  and  the  RF  output  power. 

1.  Compute  the  cathode  emission  at  each  secondary-emitting 
site  using  the  available  charge,  the  local  electric  field,  and 
the  coordinates  of  the  charge  in  the  anode-sole  region  (see 
Section  V) . 

2.  Place  the  emitted  charge  on  the  cathode  surface  with 
a  user-chosen  initial  energy  (normally  2eV) . 

3.  For  each  particle  in  the  anode-cathode  region: 

(a)  Store  the  coordinates  (r,S)  in  arrays  RP  and 
THP .  The  velocity  (dr/dt,  dl/dt)  is  already  ir. 
arrays  VRP  and  VTP.  Advance  the  position  by  one 
time  step.  Use  the  local  electric  field  from 
arrays  ERDBP  and  ETDBP  which  contain  the 
components  (E  /B,  E^/'B),  where  3  is  the  magnetic 
field. 

(b)  Test  for  interception  on  the  anode  and  cathode. 

If  the  particle  has  intercepted  the  cathode, 
recompute  its  final  velocity  using  a  smaller  time 
step  to  obtain  an  accurate  impact  energy. 

Evaluate  the  local  secondary  charcre  produced. 


57 


(c)  For  charges  not  intercepted,  store  the  predicted 
position  (r,9)  for  time  t  +  At  in  arrays  R  and 
TH,  and  increment  the  induced  direct  current. 

Also  assign  the  local  space-charge  density  to 
array  WSC. 

4.  Solve  Poisson's  equation  for  the  space-charge  potential 
in  the  same  array  WSC  (actually  normalized  to  avoid  machine 
overflow) . 

5.  Derive  the  electric-field  components  over  the  space- 

charge  mesh  in  arrays  ERSC  and  ETSC,  which  contain  the  component 

(E  /B,  E./B). 

r  o 

6.  Store  the  radial  space-charge  electric  field  at  the 
cathode  and  one  mesh  interval  above  the  cathode  to  control  the 
emission  in  the  next  step,  using  array  SCCAT.  This  array,  but 
not  the  full  space-charge-field  array,  is  included  in  the 
continuation  data  file  used  to  run  a  simulation  in  several 
stages . 

7.  Derive  the  convection  currents  and  electrode  impact 
powers.  Store  the  induced  direct  current  at  time  t  +  At 

in  the  scalar  CURCCT. 

8.  Advance  the  anode  voltage  VAVARY  and  the  modulator 
current  DCTOT  from  time  t  to  time  t  +  At  using  the  induced 
current  CURCCT  in  a  predictor-corrector  computation. 

9.  The  RF  network  equations  are  advanced  in  NCSTEP 

time  steps  subdividing  the  main  time  step  At.  The  following 

steps  are  performed  to  advance  the  vane  voltages  ( l<n<NVANES ) 

and  their  derivatives  dV  /dt  (arrays  WRF  and  DVRF)  by  one 

n 

such  intermediate  step. 

(a)  Advance  the  voltage  Vn  and  its  derivative  dVn/dt 
in  the  predictor  step. 


(b)  For  the  network  currents,  obtain  the  derivatives 

dl  /dt  at  vane  n  using  the  stored  values  for  I  at 
n  n 

times  t,  t  -  At,  t  -  2At,  and  t  -  3At  (arrays 
CURIND,  CURP,  CURPP  and  CURPPP) .  Here  the  value 
of  dl  /dt  at  time  t  +  (NC/NCSTEP) At  ( l£NC<NCSTEP) 
is  obtained  by  polynomial  extrapolation. 

(c)  Obtain  the  second  derivative  d2V  /dt2  from  the 

n 

network  matrix  equation  (Appendix  I) .  Apply  the 

corrector  to  obtain  dV  /dt  and  V  . 

n  n 

10.  For  a  backward-wave  amplifier,  change  the  sign  of 

the  RF  voltage  on  the  even-numbered  vanes  to  obtain  the  voltage 
seen  by  the  beam  (in  array  WBMDB  as  Vn/B)  . 

11.  For  each  particle  in  the  anode-cathode  region: 

(a)  Evaluate  the  local  electric  field  (direct,  RF, 
and  space-charge  components)  in  temporary  scalars 
ERDBA  and  ETDBA,  using  the  predicted  position 

(r ,3)  at  time  t  +  at. 

(b)  Average  this  field  with  the  field  stored  in 
arrays  ERDBP  and  ETDBP  from  time  t  (before  the 
predictor  step) . 

(c)  Advance  the  trajectory  again  using  arrays  RP,  THP, 
VRP  and  VTP,  repeating  the  predictor  step  but  with 
the  newly  averaged  electric  field.  No  test  for 
interception  is  necessary  since  most  intercepted 
particles  have  already  been  eliminated.  Any 
further  interception  is  computed  in  the  succeeding 
time  step. 

(d)  Recompute  the  local  electric  field  at  the  "corrected 
position  and  simultaneously  compute  the  induced 

current  I  on  each  vane  (array  CURINB) .  Store  the 

n 

final  position,  veioctiy  and  electric  field  in 


arrays  R,  TH,  VR,  VT,  ERDBP  and  ETDBP  for 
time  t  +  At. 

12.  For  the  backward-wave  tube,  change  the  signs  of 

the  RF  induced  current  on  the  even  vanes  to  obtain  the  current 
(array  CURIND)  seen  by  the  network. 

13.  If  more  time  steps  are  required,  return  to  stage  1 
for  the  next  step. 


SECTION  XI 


RESULTS  OF  THE  SIMULATIONS 

The  distributed-emission  CFA  program  has  been  applied 
to  study  three  existing  tubes,  and  the  results  are  described 
here.  The  corresponding  experimental  data  are  in  Appendices 
B-D  for  the  Raytheon  QKS1842  backward-wave  tube  and  for  the 
Varian  SFD-261  and  the  Raytheon  QKS1319  forward-wave  tubes. 

There  is  reasonable  agreement  between  the  computations  and  the 
measurements.  The  numerical  results  are  summarized  in  Table 
6  and  more  details  are  given  below. 

A.  QKS1842  Backward-Wave  CFA 

Typical  charge  distributions  in  the  QKS1842  backward-wave 
amplifier  are  shown  in  Figures  11  and  12.  The  cylindrical 
anode-cathode  region  is  displayed  for  convenience  as  two 
rectangular  plots,  from  angle  0  to  r  radians  and  tt  to  2n  radians 
around  the  tube.  The  RF  output  is  at  an  angular  position  of  0 
radians  and  the  sever  region  at  the  input  end,  near  2 tt  radians. 
In  the  actual  calculation  the  charge  distribution  is  of  course 
continuous.  The  electron  stream  moves  from  left  to  right, 
close  to  synchronism  with  the  interacting  space  harmonic  of 
the  RF  wave,  but  in  the  opposite  direction  to  the  power  flow 
in  this  backward-wave  amplifier. 

There  are  two  distinct  modes  of  operation.  The  charge 
hub  covering  most  of  the  cathode  and  the  distinct  charge  spokes 
near  the  anode  (Figure  11)  form  as  expected  in  the  high-power 
mode.  However,  at  power  levels  of  about  168  k'W,  as  in  normal 
tube  operation,  the  beam  develops  an  irregular  lumped  charge 
distribution  (Figure  12).  The  high-power  mode  was  achieved 
in  practice  with  a  more  highly  emitting  cathode  (gold  and 
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Figure  11.  Beam  profile  computed  in  QKS1842  high-power  mode. 


Charge  distribution  after  step  IS84 


Figure  12(a).  Beam  profiles  computed  in  QKS1842  low-power  mode 


CHAME  DISTRIIUTION  AFTER  STEP  «»• 


12(b).  Beam  profiles  computed  in  QKS1842  low-power  mode 


magnesium  oxide)  than  the  cermet  or  platinum  surface 
normally  used.  The  RF  drive  power  was  raised  from  5.5  kW 
to  19  kW  to  suppress  oscillations  and  the  magnetic  field 
lowered.  RF  powers  up  to  500  kW  were  obtained  at  about  50 
percent  efficiency,  compared  with  only  168  kW  at  39  percent 
efficiency  in  the  low-power  mode.  The  existence  of  these  two 
modes  was  unexpected  before  these  simulations  were  performed. 

They  provide  an  explanation  for  design  problems  encountered 
in  the  QKS1842  tube,  such  as  delayed  starting  and  high  cathode 
heating. 

Previous  simulations  without  the  modulator  have  consistently 
shown  these  two  types  of  mode.  When  the  spokes  of  charge 
develop  they  transfer  to  the  RF  field  a  high  instantaneous 
power  ranging  between  168  kW  and  770  kW.  The  irregular  groups 
of  charge  deliver  only  about  150  kW  or  less.  Figure  13  is  a 
further  example,  where  the  beam  is  transferring  only  44  kW  to  the 
wave.  These  groups  show  a  tendency  to  break  into  spokes  near 
to  the  anode,  suggesting  that  the  same  cold  RF-network  wave  is 
being  amplified  in  both  cases. 

Apparently  in  the  low-power  mode  the  spokes  of  charge 
develop  intermittently  and  then  coalesce  at  irregular  intervals. 


1 .  High-Power  Mode 

Numerical  results  for  the  high-power  mode  are  summarized 
in  Figure  14  and  Table  7. 

The  modulator  open-circuit  voltage  of  44,000  V  and  short 
resistance  of  458  Cl  are  chosen  so  that  the  load  line  passes 
through  the  operating  point  of  48  A  at  22,000  V.  In  fact  in 
all  three  cases  a  higher  voltage  (e.g.  27,900  V)  and  lower- 
current  (e.g.  35  A)  develop  in  the  simulation.  The  measured 
magnetic  field  of  0.377 T was  obtained  using  a  solenoid  calibratio 


Fiqure  13.  Irrecjular  charge  distribution  computed  in  QKS1842  with 
unsmoothed  secondary  emission  and  9,576  interacting  particles. 
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made  without  the  CFA  and  is  believed  to  be  correct.  Reducing 
the  magnetic  field  by  four  percent  (point  3)  reduces  the  voltage 
by  only  2.6  percent  to  27,180  V. 

The  RF  output  power  is  displayed  as  a  function  of  time  in 
Figure  15.  The  "1-period  mean"  is  the  fundamental  Fourier 
component  from  a  single  RF  period.  The  "15-period"  mean  is  the 
moving  average  covering  one  transit  time  around  the  tube  for 
the  cold-circuit  backward  wave.  The  high  transient  power  at 
the  start  is  due  to  the  open-circuit  modulator  voltage.  After 
this  power  has  been  dissipated  in  the  external  load  terminating 
the  network,  a  dynamic  equilibrium  develops.  The  fluctuations 
are  apparently  due  to  the  reverse-directed  circuit  wave.  They 
appear  in  both  forward-wave  and  backward-wave  tubes. 

The  corresponding  current  delivered  by  the  modulator 
(Figure  16)  is  almost  constant.  This  external  current  has  been 
smoothed  by  the  shunt  capacitance  of  the  CFA.  Counting  the 
actual  charges  emitted  and  collected  by  the  cathode  gives  a 
nee  current  that  fluctuates  between  26  A  and  44  A.  However, 
over  15  RF  periods  the  average  values  from  these  two  alternative 
current  calculations  agree  to  within  1.3  percent. 

The  upper  lineprinter  plots  of  Figure  17  show  the  RF  power 
at  each  of  the  35  active  vanes,  averaged  over  the  preceding 
15  RF  periods.  The  standing-wave  pattern  towards  the  output 
(vanes  1-18)  is  due  to  the  reflected  wave.  With  a  higher  output 
power  or  a  lower  RF  drive  power  the  tube  would  approach osci 1 lation 

The  lower  plot  of  Figure  17  shows  the  phase  lag  of  the 
induced  current  from  the  peak  voltage  on  each  vane,  again 
averaged  over  15  periods.  At  vane  1  (RF  output)  the  induced 
circuit  current  and  the  RF  voltage  are  117°  out  of  phase 
because  the  beam  contains  spokes  of  charge  recirculating  after 
passing  the  drift  region  (beyond  vane  35).  Clearly  the  modulatic 
fed  through  on  the  beam  is  appreciable  (s.e  also  Ficur-.  ID. 
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Figure  17.  RF  power  and  phase  computed  around  the  QKS1842 
in  the  high-power  mode. 


In  this  cylindrical  model  the  anode-cathode  region  seen 
by  the  beam  is  continuous.  In  the  rectangular  model,'*''*'  the 
recirculating  rods  must  be  shifted  in  phase  to  make  the  transit 
time  of  the  interaction  region  an  integer  multiple  of  the  RF 
period.  No  such  correction  is  needed  here  because  the  entire 
RF  network  is  viewed  at  every  time  step. 


2 .  Low- Power  Mode 

In  this  mode  the  RF  output  power  fluctuates  because  of  the 
irregular  beam  shape  (Figure  18)  ,  while  the  anode  current  is 
almost  constant  (Figure  19) .  The  numerical  results  (Table  8 
and  Figure  20)  are  well  within  the  range  of  the  measurements 
especially  as  the  magnetic  field  is  known  only  within  about  5 
percent . 

The  peaks  of  the  RF  power,  such  as  at  period  138, 
correspond  to  formation  of  the  spoke  mode  of  charge  distribu¬ 
tion.  The  spokes  appear  only  intermittently  in  the  charge 
distributions  displayed  by  the  program.  For  most  of  the  run 
the  particles  form  an  irregular  pattern  such  as  in  Figure 


3.  Results  at  Constant  Voltaqe 


Before  the  modulator  circuit  was  included  in  the  model 
the  anode  voltage  was  fixed  and  the  anode  current  was  unconstrained. 
The  plotted  power  and  current  of  Figures  21  and  22  demonstrate 
clearly  the  "runaway"  effect  that  can  be  caused  by  even  a  small 
numerical  error.  Excess  charge  driven  to  the  anode  leads  in 
turn  to  higher  RF  fields.  The  modulator,  however,  acts  by 
reducing  the  direct  voltage  to  oppose  growth  of  the  current. 


B .  SPD-261  Forward-1, 1 Vave  CFA 

This  tube  was  selected  for  study  because  of  the  experimental 
data  available  for  instrumented  tubes,  and  for  comparison  of 
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power  computed  over  250  periods  in  QKS1842  low- power  mi 


Figure  22.  Anode  current  computed  in  the  QKS1842  CFA  at  24.3  kV  and 
0.4114  T  with  10, 000-rod  model  and  full  secondary  emission.  Demonstrate 
runaway  effect  in  absence  of  modulator. 


COMPUTED  AND  MEASURED  RESULTS  FOR 
QKS1842  BACKWARD-WAVE  AMPLIFIER  IN  LOW-POWER  MODE 
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the  present  computed  results  with  those  of  Varian  Associates . ^ 

The  computed  currents  and  powers  are  higher  than  measured, 
as  Figure  23  and  Table  9  show.  The  Varian  computer  program  has 
predicted  an  output  power  of  1.4  at  an  anode  current  of 
17.7  A,  compared  with  the  measured  1.65  PQ  at  22  A  (Point  C, 
on  Figure  23) .  These  values  are  averaged  over  five  passes  around 
the  CFA.  The  Varian  results  for  a  single  pass  vary  between 
1.33  PQ  at  16.3  A  to  1.53  Pq  at  19.1  A.  They  lie  above  the 
measured  V-I  curve  while  the  present  results  lie  below  it. 
Increasing  the  magnetic  field  by  3  percent  (the  limit  of 
experimental  error)  reduces  the  computed  current  from  27.0  A 
to  23.5  A  but  the  output  power  of  2.33  Pq  is  still  above  the 
measured  1.65  PQ. 

As  in  the  QKS1842  simulations,  the  modulator  constrains 
a  dynamic  steady  state.  Figures  24  and  25  show  the  evolution  of 
a  simulation  over  time.  The  transit  time  per  pass  around 
the  tube  with  the  cold  RF  phase  velocity  is  about  18  RF  periods, 
so  that  5 u  passes  are  covered  here. 

The  instantaneous  charge  distribution  (Figure  26)  shows 
that  the  spokes  of  charge  formed  in  the  RF  field  are  passed 
through  the  drift  space  and  cause  anode  bombardment  even  in  the 
region  of  low  RF  power,  the  first  20  vanes.  In  experimental 
tubes  the  anode  vane  temperature  is  observed  to  peak  near  vane 
40  (3.5  radians  from  the  input)  when  the  cathode  is  mechanically 
centered.  Moving  the  cathode  closer  to  the  anode  at  the  RF 
input  shifts  the  peak  temperature  back  to  vane  27  of  the 
instrumented  tube.  Future  simulations  will  display  both  the  RF 
power  and  the  anode  heating  or.  each  vane. 

Variar.  Associates  have  concluded  that  amplification  m  the 

collection  of  recirculated  charge  near 
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SFD-261  is  due  first  to 


■u tod  in  SFD-261  forward-wave  CFA.  Modulator 


RUN  NO.  23029 


computed  11  CFD--261  forward-wave 
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Figure  26.  Beam  profile  computed  in  SFD-261 
forward-wave  DECFA. 
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the  RF  input  and  then  to  the  current  extracted  directly  from 
the  hub.^  However,  it  is  difficult  to  distinguish  these  two 
regions  in  the  present  calculations,  perhaps  because  of  the 
limits  of  the  numerical  accuracy.  In  fact  Figure  27,  which 
shows  the  RF  power  and  phase  distributions  over  the  62  vanes, 
predicts  a  more  rapid  rise  over  the  first  half  of  the  circuit 
than  measured  (see  Figure  18,  p.  58  of  ref.  10). 

As  the  output  power  is  raised,  the  measured  V-I  curve 
bends  upwards  and  the  tube  becomes  more  noisy.  Point  7  on 
the  extreme  right  of  Figure  23  was  obtained  in  an  attempt  to 
reproduce  this  behavior.  Here  the  cathode  is  believed  to  be 
emission- limited,  and  unable  to  supply  sufficient  current  by 
secondary  emission.  The  actual  limit  of  current  measured  can 
be  as  high  as  40  A  at  the  low  end  of  the  frequency  band,  but 

depends  on  the  cathode  in  the  particular  tube  under  study. 

The  most  significant  feature  here  (see  Table  9)  is  the  computed 
mean  impact  energy  and  effective  secondary-emission  coefficient, 
both  of  which  are  higher  than  at  the  lower  currents  in  the 
other  runs.  Apparently  there  is  a  reduced  charge  density  near 
the  cathode  as  more  of  the  emitted  charge  is  pulled  towards  the 
anode,  and  hence  less  charge  is  available  to  produce  secondaries 
and  maintain  the  beam. 

C .  QKS1319  Forward-Wave  CFA 

This  tube,  which  has  not  previously  been  simulated  by 

computer,  has,  of  all  the  tubes  studied,  the  highest  ratio  of 

RF  voltage  to  anode  voltage  (a  range  of  0.12  to  0.65)  and  the 
lowest  ratio  of  anode  voltage  to  cutoff  voltage  (only  0.33). 

A  typical  beam  profile  (Figure  28)  shows  distinct  but  uneven 
charge  spokes  and  a  narrow  hub  of  charge. 


The  cathode  does  not  emit  secondary  charge  for  the  final 
l/9th  of  its  circumf erence ,  but  sufficient  charge  is  generated 
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Figure  27.  RF  power  and  phase  computed  around  the 
SFD-261  DECFA. 
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before  this  region  to  cross  the  gap.  The  beam  feeds  significant 
RF  current  through  to  the  input,  as  in  the  other  tubes  studied. 

The  numerical  results  are  summarized  in  Figure  29  and 
Table  10.  Computed  points  1,  3  and  4  lie  on  the  same  load 
line  as  the  measured  point  A  but  the  computed  currents  and  powers 
are  high.  This  tube  clearly  provides  a  severe  but  useful  test 
of  the  accuracy  of  the  computations.  At  present  an  increase  of 
the  magnetic  field  from  the  actual  0.1150 T to  0.1300  T 
(point  1  on  Figure  29)  is  necessary  for  close  agreement. 

However,  the  anode  voltage7  of  10,100  V  is  then  only  0.99  times 
the  Hartree  voltage,  and  the  large  fluctuations  of  computed 
output  power  (Figure  30)  are  probably  unrealistic.  The  compu¬ 
tation  should  be  repeated  after  adjustments  have  been  made  to 
the  model. 

D.  Cathode  Phenomena 

The  QKS1842  is  normally  operated  in  the  lower  power  mode. 

Here  as  much  as  14  percent  of  the  input  power  appears  as 
cathode  heating  due  to  backbombardment ,  whereas  5  percent  is 
typical  of  other  CFAs.  The  simulation  provides,  for  the  first 
time,  an  explanation  of  this  observed  anomaly.  Notice  in  Table 
11  that  the  QKS1842  low-power  mode  gives  the  highest  cathode 
impact  energies  and  the  highest  relative  values  of  backbombardment 
power.  Apparently  two  factors  are  involved. 

First,  the  normal  effect  of  space  charge  is  to  reduce  the 

distance  of  penetration  of  electrons  into  the  RF  field  and  hence 

21 

reduce  the  impact  energy  from  the  ballistic  value.  In  the 
QKS1842  low-power  mode,  however,  the  reduced  space  charge  at 
the  voids  near  the  cathode,  increases  the  impact  energy  of 
returning  electrons  by  allowing  them  to  experience  the  higher 
RF  field  closer  to  the  anode  circuit.  Second,  the  angular 
space-charge  field  near  the  cathode  varies  more  than  in  the 
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Figure  30.  RF  output  power  computed  in  QKS1319  DECFA  with  magnetic  field  raised 
to  0.13  T . 
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Run  .1150  9,520  37.0  3,500  206,000  Complete  results 

>19004;  not  available 

IV  i  1  oils 


CATHODE  CHARACTERISTICS  COMPUTED  IN 
DISTRIBUTED-EMISSION  CROSSED- FI ELD  AMPLIFIERS 
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more  uniform  stream  of  the  high-power  mode. 

The  effect  of  changing  the  cathode  emitter  is  demonstrated 
in  a  second  simulation  with  the  same  modulator,  magnetic  field 
and  RF  drive  but  a  lower  emitting  cermet  cathode.  This  causes 
twice  as  much  backbombardment  power  as  the  more  highly  emitting 
gold  and  magnesium  oxide  surface  and  the  computed  efficiency 
(lower  than  measured)  falls  from  39%  to  36%.  A  similar  effect 
is  observed  in  a  magnetron  with  a  thermionic  cathode,  in  which 
the  backbombardment  energy  rises  as  the  available  cathode  current 
is  reduced  by  lowering  the  cathode  temperature.  Here  too,  it 
is  hypothesized  that  lower  space  charge  density  permits  cycloiding 
electrons  to  gain  more  energy  from  the  RF  circuit  by  moving 
closer  to  the  anode. 

The  program  provides  an  estimate  of  the  maximum  current 
available  from  the  cathode  at  the  emission  boundary. 

The  effective  secondary-emission  coefficient  shown  in 
Table  11  increases  as  the  emission  boundary  is  approached. 

This  coefficient  is  actually  computed  by  counting  the  electrons 
that  escape  from  a  partially  space-charge-limited  cathode. 

Hence  it  is  less  than  that  from  the  emission  table  supplied 
for  the  cathode  material.  As  the  current  drawn  to  the  anode 
is  increased,  this  effective  coefficient,  a,  increases  and  a 
smaller  fraction,  1/a,  of  the  emitted  charge  returns  to  the 
surface.  Thus  a  value  of  a  close  to  unity  implies  that  ample 
charge  exists  to  maintain  the  electron  stream,  whereas  a  larger 
value  indicates  that  the  cathode  is  close  to  emission  limitation. 

The  highest  values  of  a  are  computed  in  the  two  forward- 
wave  tubes  (the  SFD-261  and  QKS1319  tubes) .  These  tubes  also 
have  the  highest  ratios  of  RF  voltage  to  direct  anode  voltage 
(0.34  and  0.65  compared  with  about  0.14  in  the  QKS1842)  so  that 
a  larger  fraction  (1-1/a)  of  the  emitted  charge  is  drawn  to 
the  anode.  In  practice  the  QKS13.9  shows  an  emission  boundary 
at  an  anode  current  of  about  30  A;  the  SFD-261  can  provide  up 
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to  40  A  at  the  low  end  of  its  frequency  band  but  less  at 
higher  frequencies. 


A  further  example  of  the  cathode  modelling  is  the  comparison 
(Table  12)  between  runs  for  the  QKS1842  low-power  mode  with  cermet 
and  platinum  cathodes.  Here  the  anode  voltage  is  a  constant 
24.3  kV  and  cathode  emission  is  limited  by  applying  Child's 
Law  over  the  first  mesh  interval  normal  to  the  surface. 

Replacing  the  platinum  cathode  normally  used  in  production  tubes 
by  the  higher-yield  cermet  cathode  reduces  the  mean  impact  energy 
from  430  eV  to  287  eV.  The  anode  current  falls  (at  constant 
voltage)  from  16  A  to  14.4  A  as  more  current  returns  to  the 
cathode,  keeping  the  net  backbombardment  power  relatively 
unchanged.  In  practice,  however,  with  the  modulator  to  regulate 
the  current,  no  significant  differences  in  tube  performance  have 
been  observed. 

E.  RF  Oscillation  and  Mode  Boundary 

Competing  oscillations  limit  the  output  power,  gain  and 
operating  bandwidth  of  practical  CFA's.  The  situation  is  certainly 
more  complex  in  the  tube  than  in  the  model,  which  now  reproduces 
the  RF  circuit  dispersion  and  impedance  only  at  the  drive 
frequency.  However,  the  model  includes  both  forward  and 
reverse-directed  power  flow,  and  can  produce  standing  waves 
as  a  result.  Figure  31  shows  such  a  standing  wave  between 
vanes  26  and  45  in  the  QKS1319  where  the  computed  output  power 
is  167  kW  with  an  anode  current  of  28. 8A.  The  RF  output  power 
shows  corresponding  fluctuations  over  time  and  approaches  zero 
in  certain  RF  periods  (Figure  30) . 

The  RF  network  elements  are  chosen  so  that  the  cold  circuit 
is  matched  by  the  equivalent  of  an  infinite  line  at  each  end. 
Therefore  these  standing  waves  must  arise  from  power  transferred 
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Figure  31.  RF  power  and  phase  distributions 
computed  in  QKS1319  DECFA  as  averages  over  14 
periods. 


to  the  network  from  the  beam.  A  wave  flows  in  both  directions 
fron  each  node  as  a  result,  but  normally  only  the  forward  wave 
is  in  synchronism  with  the  beam.  The  wave  returning  to  the  input 
is  attenuated.  If  more  power  is  returning  to  the  RF  input 
than  is  being  supplied  by  the  network  driving  current,  the 
computed  net  power  delivered  to  the  network  is  negative  (Appendix 
I).  This  condition  is  interpreted  as  an  oscillation  at  the 
drive  frequency.  The  difference  between  the  cold-circuit 
drive  power  and  the  net  power  delivered  to  the  circuit  in  the 
presence  of  the  beam  is  an  estimate  of  the  minimum  drive  power 
required  for  stability. 

The  numerical  results  (Table  13)  show  good  qualitative 
agreement  with  the  available  measurement  data,  but  quantitative 
estimates  of  the  mode  boundaries  require  a  more  accurate 
calculation  of  the  actual  operating  current  and  voltage . 

Both  in  laboratory  tests  and  in  these  simulations  the 
QKS184  2  high-power  mode  was  stable  with  an  RF  drive  power  of 
19  kw.  Notice  however,  that  Figure  15  shows  significant 
oscillations  of  the  RF  output  power  over  time.  Figure  17  shows 
the  corresponding  standing  waves  on  vanes  1-23  (towards  the  RF 
output  of  this  backward-wave  circuit) .  The  maximum  power  of 
8,500  W  returning  to  the  RF  input  during  the  simulation  is 
above  the  5,500  W  drive  level  that  is  available  for  practical 
operation.  In  the  QKS1842  low-power  mode,  however,  (run  42077, 
Figure  18)  the  backward  power  at  the  input  is  only  about  500  W 
with  an  RF  drive  of  5,500  W.  This  result  implies  that  the 
tube  will  not  oscillate,  at  least  at  the  drive  frequency. 

At  the  two  highest-power  points  computed  for  the  SFD-261 
about  half  the  RF  drive  power  is  returning  to  the  input  vane. 

The  low-power  point  (0.65  Pq  output)  has  a  corresponding 
reverse-directed  power  of  less  than  one  percent  of  the  drive. 
These  numbers  agree  only  qualitatively  with  the  measurements. 


COMPUTED  AND  MEASURED  OSCILLATION  CONDITIONS 


Frei|iionry  reduced  to  1.25  GHz,  but  magnetic  field  too  high  for  synchronism. 


They  predict  stability  with  currents  up  to  39  A,  whereas  the 
instrumented  Varian  tube  can  produce  only  20  A  without  oscil¬ 
lation.  At  the  rated  22  A  an  oscillation  signal  develops 
about  1  GHz  from  the  band  center.  Since  the  tube  is  stabilized 
at  22  A  by  using  a  non-circular  cathode  (as  in  production 
models)  it  appears  that  the  recirculating  charge  is  responsible. 
The  barely  visible  standing  wave  on  vanes  4-9  in  Figure  27 
may  be  a  warning  of  trouble.  It  would  be  interesting  to  repeat 
the  simulation  using  the  actual  cathode  shape  for  the  production 
version  of  the  SFD-261.^ 

In  the  QKS1319  the  measured  upper  mode  boundary  limits 
the  measured  current  to  about  30  A  with  140  kW  RF  output  power. 

In  the  simulation,  oscillations  actually  appear  for  anode 
currents  above  28  A,  at  which  the  returning  RF  power  may  be  as 
high  as  8,700  W  for  3,500  W  drive.  Probably  one  cause  is  the 
high  interaction  impedance  of  84  Q  at  1.3  GHz  and  146  Q  at 
1.25  GHz.  With  incorrect  beam-wave  synchronism  (run  19008) 
the  oscillation  is  computed  with  an  RF  output  power  of  only  33  kW. 
It  is,  however,  still  necessary  to  verify,  with  a  more  accurate 
numerical  model,  that  the  measured  operating  points  are  stable. 

The  Harris  SAI  model  is  significantly  better  than  other 
CFA  models  which  neglect  the  reverse-directed  power  entirely. 
However,  it  predicts  only  oscillations  at  the  drive  frequency, 
while  practical  tubes  produce  spurious  out-of-band  signals 
at  power  thresholds  below  those  now  computed.  A  more  accurate 
correlation  of  theoretical  and  measured  output  power  must  next 
be  obtained,  particularly  for  the  forward-wave  tubes.  Then  the 
computed  standing-wave  pattern  on  the  vanes  near  the  RF  input 
may  be  useful  as  an  additional  indication  of  mode  problems  in 
a  design. 


F.  Sources  of  Computational  Error 

The  preceding  results  have  shown  two  areas,  probably 
related,  where  the  quantitative  results  can  be  improved. 

Firstly,  the  anode  current  and  RF  output  power  computed  for  a 
given  voltage  are  above  the  measured  values  for  the  forward- 
wave  tubes,  and  below  the  measured  values  for  the  backward- 
wave  tube.  Secondly,  the  power  balance  is  unsatisfactory 
at  present. 

Consider  the  input  and  output  powers  averaged  over  several 
RF  periods.  One  or  more  transit  times  around  the  tube  are  used 
in  the  computations.  In  a  steady  state, the  conservation 
equation  is: 


Direct  input  power  from  modulator  +  RF  drive  power 

=  RF  output  power  +  anode  dissipation  power 
due  to  circuit  loss  +  anode  bombardment  power 

+  cathode  backbombardment  power.  (1) 


The  measured  powers  satisfy  this  relation  within  the 
limits  of  experimental  error. 

The  computed  results  of  Tables  7-10  show  a  3  to  12  percent 
deficit  of  total  output  power  in  the  backward-wave  tube  and 
up  to  49  percent  excess  total  output  power  in  the  2  forward- 
tubes,  averaged  over  several  RF  periods  that  show  an  approximate 
steady  state . 

There  are  several  possible  sources  of  these  computational 
errors . 

(1)  The  anode  bombardment  power  is  over-estimated  by 
about  3  percent  of  the  total  power  in  the  backward-wave 
simulations  and  by  about  13  percent  in  the  forward-wave  tubes. 
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At  the  cathode  the  trajectories  of  intercepted  electrons  are 
recomputed  with  a  reduced  time  step  so  that  the  backbombardment 
power  is  accurate.  However,  the  present  program  makes  no  such 
correction  at  the  anode.  Future  versions  of  the  program  will 
compute  the  impact  energy  accurately  at  both  the  cathode  and 
the  anode . 

(2)  It  was  suspected  earlier  in  the  study  that  the  model 
may  be  using  insufficient  simulation  particles.  However,  tests 
have  been  performed  with  as  many  as  26,000  particles  (Section 
XII)  and  the  results  imply  that  the  normal  7,000  particles  are 
sufficient  in  all  three  tubes. 

(3)  The  predictor-corrector  trajectory  calculation,  which 
is  not  time-centered,  may  introduce  small  but  cumulative 
errors  into  the  space-charge  forces.  In  Section  VI  a  modified 
trajectory  calculation  is  proposed. 

(4)  The  induced  RF  current  on  each  vane  includes  terms 
both  from  the  charge  in  the  interaction  region  and  from  the  charge 
collected  on  the  vane.  It  is  now  believed  that  the  collected 
current  should  be  omitted.  This  correction  will  not  affect  the 
observation  of  two  modes  in  the  beam  of  the  QKS1842,  since  the 
"groups  of  charge"  such  as  in  Figure  13  develop  at  the  start  of 

a  run  before  significant  charge  is  collected  on  the  vanes. 
However,  it  could  have  an  important  effect  on  the  quantitative 
results . 

(5)  The  time  step  is  now  1/15,  less  than  1/10,  and 
1/6  cyclotron  period  in  the  QKS1842,  SFD-261  and  QKS1319 
respectively.  Varian's  SFD-261  simulations  use  1/10  cyclotron 
period.  Initial  trials  have  shown  that  the  QKS1842  time  step 
is  sufficiently  small.  For  the  two  forward-wave  tubes  further 
tests  should  be  made  with  a  smaller  step. 


The  QKS1842  was  designed  using  parametric  equations 
to  scale  down  in  size  a  similar  backward-wave  CFA,  the 
QKS1705  which  operates  over  the  same  9.5-10  GHz  frequency  band. 
The  smaller  anode  circumference  and  smaller  vane  spacing  of  the 
QKS1842  give  it  a  reduced  RF  phase  velocity  and  a  reduced 
interaction  impedance.  Consequently,  the  QKS1842  requires  a 
lower  operating  voltage  and  a  higher  magnetic  field  to  provide 
the  beam  drift  speed  for  synchronism.  It  now  supplied  175  kW 
peak  power  at  17.1  A  over  a  reduced  bandwidth,  with  a  5.5  kW 
RF  drive  at  about  39%  efficiency.  This  is  the  low-power 
mode  demonstrated  in  the  simulations.  The  tube  shows  anomalously 
high  cathode  backbombardment  and  has  also  given  difficulties 
with  starting. 


The  parameters  of  the  QKS1842  and  QKS1705  are  compared 
in  Table  14,  together  with  those  of  the  QKS1319  and  Varian 
SFD-261  tubes.  Clearly  the  QKS1705  and  QKS1842  have  the  lowest 
interaction  impedances  and  also  the  lowest  ratios  of  RF  voltage 
to  direct  voltage.  In  addition,  the  QKS1842  is  operating  closest 
to  the  Hartree  threshold  voltage^ ' (even  below  this  if  the 
measured  0.4114  T  magnetic  field  is  accurate).  These  two 
observations  may  account  for  the  irregular  beam  bunches  that 
develop  in  the  low-power  mode  in  the  QKS1842  and  hence  for  the 
high  backbombardment  and  starting  delay.  Detailed  studies  on 
the  SFD-261^1  have  in  fact  shown  that  for  a  given  RF  field 
there  is  a  minimum  ratio  of  anode  voltage  to  Hartree  voltage 
allowing  spokes  of  charge  to  form.  In  the  high-power  mode  of 
the  QKS1842  the  RF  field  is  sufficient  to  produce  the  expected 
behavior,  and  the  efficiency  is  increased  to  50%.  However,  this 
mode  is  unsuitable  for  practical  use  because  of  the  high  RF 
drive  power  (19  kW)  needed  to  suppress  oscillation.  The  5.5  kW 
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actually  used  is  necessitated  by  the  required  gain  and  by 
the  available  driver  tubes. 


SECTION  XII 


PROGRAM  REQUIREMENTS 

The  computing  resources  required  are  controlled  primarily 
by  the  number  of  particles  (rods) ,  the  time  step,  and  the 
length  of  the  run. 

A.  Number  of  Simulation  Particles 

Table  15  summarizes  the  particle  requirements  for  these 
runs.  The  charge  per  rod  and  the  number  of  rods  used  are  controlled 
by  the  input  parameter  NBRILL,  the  number  of  rods  required 
to  fill  the  cylindrical  anode-cathode  region  with  a  density 
such  that  the  electron  plasma  frequency  equals  the  cyclotron 
frequency.  This  is  the  Brillouin-f low  condition  for  a  crossed- 
field  tube. 

In  trials  with  varied  resolution  and  both  with  and  without 
the  modulator  the  number  of  simulation  rods  has  been  increased 
and  the  space-charge  array  size  changed  from  257  (angular) 
by  33  (radial)  to  513  by  33.  Both  in  this  tube  and  in  the 
QKS1842  and  QKS1319  tubes  the  results  are  not  sensitive  to 
the  changes.  Figures  32  and  33  are  good  examples  showing  that 
6,000  rods  are  sufficient,  although  up  to  27,900  are  tested. 

(The  modulator  is  absent  here.) 

These  curves  show  that  the  mean  currents  and  powers  over 
18  periods  agree  within  6  percent,  although  the  finer  details 
of  the  interaction  are  more  sensitive  to  the  model.  The  power 
balance  shows  an  excess  output  power  between  23  and  25  percent 
in  both  cases.  In  the  QKS1319  corresponding  figures  are  53.5  A 
and  392  kW  with  4,100  rods  compared  with  58.4  A  and  402  kW  with 
15,800  rods,  the  mean  current  and  power  over  periods  36  to  52. 

The  results  are  high  compared  with  the  measured  values  because 


TABLE  15 


PARTICLE  NUMBERS  REQUIRED  IN 
CFA  SIMULATIONS 


Number  (NBRILL) 
of  Rods  to  Fill 
Tube  at  Brillouin  Maximum 


Run  Number 

Tube 

Density 

Rods  Used 

42079 

QKS1842 

(high-power  mode) 

10,000 

7,603 

42077 

QKS1842 

(low-power  mode) 

10,000 

7,027 

42073 

QKS1842 

(low-power  mode) 

45,000 

23,471 

23029 

SFD-261 

20,000 

7,330 

23024 

SFD-261 

90,000 

27,886 

19006 

QKS1319 

20,000 

5,602 

Figure  32.  RF  output  power  computed  for  Varian  forward-wave  DECFA .  Modulator  not 
included  in  model. 


Figure  33.  Anode  current  computed  for  Varian  forward-wave  DECFA.  Modulator  not  in  model 


the  modulator  model  was  not  implemented  for  these  runs.  To 
attain  here  the  resolution  of  the  Varian  single-wavelength  model 
would  require  about  60,000  rods  and  a  2049  x  49  space-charge 
array.  Such  detail  appears  not  to  be  necessary  in  the  present 
model . 

B.  Time  Step 

The  three  criteria  used  to  estimate  the  largest  allowable 
step  size  are 

(1)  At  least  16  steps  per  RF  period 

(2)  At  least  8  steps  per  electron  cyclotron  period 

(3)  Not  more  than  one  mesh  interval  to  be  crossed  by 
a  typical  particle  in  one  step. 

The  user  actually  supplies  the  number  of  steps  per  RF 
period.  It  is  recommended  that  consistency  tests  be  performed 
with  a  given  step  and  then  one  half  that  step  whenever  a  new 
CFA  is  simulated. 

The  values  used  for  the  simulations  are  shown  in  Table  16. 
Criterion  (3)  is  satisfied  for  angular  velocities  in  the  charge 
hub  and  for  radial  velocities  less  than  1/6  the  RF  phase 
velocity.  It  is  not  generally  satisfied  in  the  spokes  where 
charge  approaches  the  anode. 

The  QKS1842  interaction  has  been  computed  over  six  RF 
periods  with  32  time  steps  per  period  in  place  ot  the  usual 
16,  and  the  powers  transferred  to  the  RF  wave  agree  within 
9  percent.  This  power  transferred  from  the  beam  to  the  RF 
circuit,  in  one  period  is  a  conveniently  sensitive  measure  of 
the  accuracy.  Similar  tests  should  be  performed  for  the 
QKS1319  and  SFD-261  tubes. 
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In  the  results  described  in  Section  XI,  between  50  and  100 
RF  periods  suffice  to  attain  a  dynamic  steady  state  in  the  two 
forward-wave  tubes.  The  RF  circuit  in  the  backward-wave  tube 
has  a  longer  response  time  to  initial  transients  because  of 
the  feedback  of  power  from  the  RF  input  once  the  beam  is  fully 
bunched  there.  The  QKS1842  high-power  mode  requires  100  RF 
periods  and  the  lower-power  mode  has  been  run  for  250  periods 
because  of  the  larger  fluctuations  in  the  charge  distribution. 

D .  Computing  Time  and  Storage 

The  computing  times  required  for  simulations  of  the  three 
tubes  are  summarized  in  Table  17.  These  figures  are  for  the 
Harris  550  computer.  Relative  values  for  other  machines  are 
estimated  below: 


Computer 


CPU  Time  (sec) 


Harris  550 
IBM  360/67 
IBM  370/148 
CDC  6500 
Amdahl  4 70V/ 6 
Amdahl  470V/7 


1 

1 

0 

2 

0 

0 


16 

51 

28 

124 

083 


The  CPU  time  per  particle  (rod)  per  time  step  is  a  useful 
figure  of  merit  for  computer  programs  cf  this  type.  It  decreases 
as  more  particles  are  introduced.  It  is  almost  independent  of 
the  number  of  vanes  that  produce  the  RF  field  in  the  CFA.  A 
figure  of  6.4  x  10  sec.  is  typical  for  the  550  computer. 

From  Table  17  it  can  be  seen  that  a  typical  100-period 
run  with  16  time  steps  per  RF  period  and  7,000  rods  of 


TABLE  17 


Tube 

QKS1842 
QKS1842 

SFD-261 

SFD-261 

QKS1319 

QKS1319 
I 

I 

I 


COMPUTING  TIME  FOR  DECFA  SIMULATION  ON 
HARRIS  550  COMPUTER 


Number  of 
Active  Vanes 
In  Tube 


Number  of 
Rods  in 
Simulation 


CPU 

Seconds  Per 
Particle  Per 
Time  Step 


CPU  Time 
(hours)  for 
100  RF  Periods 


35 

7,050 

6. 62x10" 3 

20.7 

35 

5,600 

7 . 57x10" 3 

18.8 

62 

5,740 

7 . 24xl0~3 

18.5 

62 

25,500 

5 . 52x10" 3 

62.6 

45 

3,610 

8 . 31x10" 3 

13.3 

45 

12,300 

6 . 40x10" 3 

35.0 
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charge  requires  about  20  hours  of  computing.  However,  more 
efficient  trajectory  and  RF  field  calculations  are  planned 
(Section  VI)  which  should  halve  this  time. 

With  a  maximum  of  20,000  rods  the  program  now  requires 
505,000  24-bit  words  (252,500  48-bit  floating-point  numbers) 
on  the  Harris  550  machine.  The  requirements  on  other  machines 
are  similar,  as  most  of  the  storage  is  for  particle  and  field 
arrays. 


SECTION  XIII 


CONCLUSIONS  AND  RECOMMENDATIONS  FOR  FUTURE  WORK 

The  capabilities  of  the  Harris  SAI  distributed-emission 
CFA  simulation  program  are  summarized  in  Table  18,  which 
shows  the  important  features  of  the  tube  that  a  useful 
simulation  should  reproduce. 

This  AFOSR-sponsored  study  has  advanced  the  understanding 
the  operation  of  the  DECFA,  demonstrating  particularly  the 
following  results: 

•  The  modulator  controls  the  anode  current  and  voltage 
even  over  the  small  time  scale  of  a  few  RF  periods, 
during  which  small  fluctuations  of  voltage  and  current 
occur.  A  rise  of  current  causes  a  fall  in  the 
voltage,  maintaining  stability. 

•  The  circulating  charge  carries  significant  RF 
current  across  the  drift  region  in  both  the  forward- 
wave  and  backward-waVe  tubes. 

•  The  SFD-261  and  QKS1319  forward-wave  tubes  operate 
with  the  characteristic  spoked  beam  expected  from 
analytic  theory. 

•  The  QKS1842  shows  two  types  of  operating  mode,  with 
efficiencies  of  39%  and  56%.  The  lower  efficiency  is 
obtained  in  normal  operation,  where  the  charge  distribu 
tion  is  irregular. 

•  Anomalously  high  backbombardment  power  results  from 
sparse  charge  density  in, the  QKS1842  low-power  mode. 


TABLE  18 


CAPABILITIES  OF  THE  HARRIS  SAI  DECFA  MODEL 


Feature  Result 

Anode  current  at  given  Less  than  measured  value  in  backward- 

voltage  and  given  magnetic  wave  CFA;  above  measured  value  in 

field.  forward-wave  CFA;  in  actual  runs 

the  voltage  is  allowed  to  vary. 

RF  output  power.  About  25%  below  measured  value  in 

backward-wave  CFA;  50%-100%  above 
measured  value  in  forward-wave  CFA. 


Anode  heating. 


Backbombardment 
and  cathode  heating. 


Now  overestimated  by  up  to  30  percent 
but  an  improved  computation  of 
electron  energy  at  the  vanes  will 
be  straightforward.  Power  distribu¬ 
tion  over  vanes  is  computed  and 
can  readily  be  displayed. 

Good  order-of-magnitude  agreement 
with  measurements,  including  the 
anomalously  high  value  of  15%  beam 
power  in  QKS1842  low-power  mode. 


Steady-state  power  balance  After  correction  of  the  anode 
balance  in  model.  heating  a  7%-15%  deficit  of  output 

power  is  computed  in  the  backward- 
wave  tube  and  a  9%-38%  excess  in 
the  forward-wave  tubes. 

Secondary  emission.  The  beam  hub  forms  stably.  Cathode 

emission  model  reproduces  different 
cathodes.  A  lower-yield  surface 
results  in  higher  backbombardment  as 
the  charge  density  falls.  The 
emission-free  drift  region  is  repro¬ 
duced  in  the  QKS1319. 
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TABLE  18  (cont.) 


Feature 


Result 


Cathode  emission  limit. 


Mode  selection  in  beam. 


Oscillation  and  upper 
current  mode  boundary. 


Spurious  RF  signals  away 
from  operating  band. 

Upper  power  limit  due 
to  loss  of  synchronism. 


Approached  but  not  reached  in 
simulations. 

Two  types  of  mode  (lumps  of  charge 
and  spokes  of  charge)  distinguished 
in  QKS1842  tube.  One  mode  (spokes) 
in  each  of  SFD-261  and  QKS1319 
with  the  normal  operating  voltage 
and  magnetic  field. 

RF  network  can  allow  oscillation 
at  drive  frequency.  Order-of- 
magnitude  agreement  with  measured 
RF  output  power  limit  in  QKS1319 
and  with  minimum  RF  drive  in  QKS1842 
high-power  mode,  but  SFD-261 
oscillations  not  reproduced  by 
present  model. 

Not  included  but  may  be  induced  by 
oscillation  due  to  reverse-directed 
RF  power. 

Not  yet  simulated  but  should  be 
predicted  by  model. 


Lower  mode  boundary  due  to  Not  included  with  present  RF  network, 
oscillation  at  lower  cutoff 
frequency  of  backward-wave 
circuit. 


Modulator.  Both  hard-tube  and  line-type 

modulators  are  treated,  but  with 
reduced  inductance  and  capacitance 
for  rapid  decay  of  transient 
voltages.  A  run  requires  a  prior 


TABLE  18  (cont.) 


Feature 


Result 


Starting  a  tube  pulse. 


Turn-off  at  the  end  of  a 
pulse  when  RF  drive  is 
removed  in  forward-wave 
tubes . 

Transverse  effects  (in 
the  magnetic-field 
direction) . 


estimate  of  current  and  voltage 
in  the  steady  state. 

Starting  can  be  simulated  on  an 
accelerated  time  scale.  (About 
1/1000  of  a  pulse  is  treated  in 
a  typical  run) .  Rise  of  line-type 
modulator  voltage  is  not  now  reproduced 
but  is  straightforward  to  include. 

Straightforward  to  reproduce  but 
not  yet  studied. 


All  motion  in  the  magnetic- fie Id 
direction  is  neglected,  but  the 
measured  interaction  impedance 
is  an  average  over  the  width  of 
the  device  and  accounts  for  the 
RF  field  variation  across  the 
vanes . 
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•  Increasing  the  secondary-emission  yield  at  the  cathode 
reduces  the  backbombardment  power  and  raises  the 
efficiency  at  a  fixed  RF  output  power. 

•  Reverse-directed  RF  power  along  the  vane  circuit 
creates  standing  waves  and  may  produce  oscillation 
corresponding  to  a  measured  mode  boundary. 

Consistency  tests  have  shown  that  7000  simulation  rods 
give  sufficient  resolution.  Although  the  steady-state  power 
balance  is  not  yet  satisfactory  in  the  forward-wave  tubes, 
adjustments  to  the  RF  network  and  trajectory  calculations  are 
proposed  to  improve  this.  It  should  be  possible  also  to  halve 
the  present  computing  time  of  20  CPU  hours  per  run  on  the 
Harris  550  minicomputer  (equivalent  to  about  23  hours  on  an 
IBM  360/67) . 

The  following  are  the  priorities  recommended  for  a  continua¬ 
tion  of  this  study: 

1.  Improve  the  trajectory  and  RF  network  models  (as 
described  in  detail  in  this  report)  in  order  to  obtain  results 
that  show  a  better  power  balance,  particularly  in  the  two 
forward-wave  tubes.  Verify  the  program  using  the  measured 
values  for  the  QKS1842,  SFD-261  and  QKS1319  tubes. 

2.  Test  alternative  calculations  of  RF  field  and  induced 
current  to  eliminate  the  RF  Green's  function  and  give  a  more 
efficient  computation. 

3.  Simulate  the  QKS1842  and  SFD-261  tubes  at  band-edge 
frequencies . 

4.  Simulate  the  QKS1705  backward-wave  CFA  to  obtain 
additional  verification  for  the  program. 

5.  Study  the  starting  conditions  in  the  QKS1842  on  an 
accelerated  time  scale.  Choose  the  modulator  inductance  to 
allow  the  anode  voltage  to  rise  from  zero  over  about  one  RF 
transit  time  around  the  tube. 


6.  Reduce  the  RF  drive  to  simulate  the  oscillation 
thresholds  in  the  SFD-261. 

7.  Model  the  production  SFD-261  forward-wave  CFA  with 
a  noncylindrical  cathode,  it  will  suffice  to  include  only 
the  effect  of  the  variable  spacing  on  the  applied  (direct) 
field;  for  the  space-charge  and  RF  fields  the  cathode  will 
be  treated  as  circular. 

8.  Incorporate  an  option  of  a  tapered  pitch  for  the  vane 
circuit.  The  RF  network  will  be  unchanged  if  it  is  assumed 
that  the  phase  delay  between  vanes  and  the  characteristic 
impedance  are  unaltered. 

9.  Test  a  more  general  form  of  RF  network  with  a 
dispersion  characteristic  that  includes  a  cutoff  frequency. 

The  general  equations  for  such  a  network  have  already  been 
included  in  the  model. 

10.  Develop  a  model  for  the  cathode-driven  CFA.  An 
accurate  solution  of  Poisson's  equation  is  required  including 
the  vane  shape  on  the  cathode  as  a  boundary.  A  suitable 
technique  has  been  used  at  the  Naval  Research  Laboratory  for 
magnetron  simulations  and  could  be  applied  to  the  CFA. 

The  basic  elements  of  the  model  have  been  established 
and  thoroughly  tested  in  the  present  study. 

It  is  expected  that  the  above  additional  effort  will 
produce  a  reliable  simulation  program  that  can  be  used  to 
improve  or  scale  existing  crossed-f ield  amplifiers. 
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APPENDIX  A 


DATA  FOR  THE  RAYTHEON  QKS1300  AMPLITRON 


The  parameters  tabulated  below  were  used  in  initial 

development  of  the  DECFA  program.  Some  of  the  operating  values 

were  supplied  by  the  Raytheon  Company;*  the  remainder  are  taken 

18 

from  Dombrowski  and  Price.  A  few  of  the  numbers  have  been 
changed  for  testing  purposes.  For  example,  the  high  RF  drive 
of  43  W  was  used  to  provide  strong  RF  fields  and  ensure  good 
beam-wave  synchronism  with  reduced  computing  time. 


Operating  Value  Used 

Quantity  Value  in  Test  Runs 


Number  of  vanes 

11 

11 

Number  of  active  vanes 

11 

11 

Magnetic  field  (T ) 

CM 

• 

.21 

Beam  width  (mm)  in  the 
magnetic-field  direction 

5.334 

5.334 

Anode  radius  (mm) 

2.337 

2.337 

Cathode  radius  (mm) 

0.991 

0.991 

Drive  frequency  (GHz) 

2.2825 

2.2825 

Anode-sole  voltage  (V) 

1800.0 

1800.0 

RF  drive  power  (W) 

0.4 

43.0 

Cold  Phase  delay  (rad)  per  pitch 
of  network  as  measured  in 
forward-wave  direction 

37/11 

3  TT  /  1 1 

Direction  of  wave  seen  by  beam 

Backward 

Backward 

■k 

Telephone  communication.  Dr.  J. 

Skowron  to  Dr. 

D.  MacGregor 

October  28,  1977. 


Quantity 

Operating 

Value 

Value  Used 
in  Test  Runs 

Interaction  impedance  (ft) 

Not  known 

120 

Characteristic  impedance  (ft) 

120  estimated 

202.3 

Attenuation  (dB)  per  pitch 

Not  known 

0.00312 

Ratio  of  vane  spacing  to  pitch 
(period)  on  the  anode  surface 

0.392 

0.392 

Cathode  emission 

Thermal 

Thermal 

Maximum  current  density  at 
cathode  (A/m2) 

Not  known 

1.068  x  106 

Anode  current  (A) 

0.018 

No  values  have 
been  obtained 

RF  output  power 

25  W 

from  present 
calculations . 

Time  step 

1/64  RF  period 

Pitch  of  circuit  (mm) 

1.335 

1.335 

I .  Basic  Tube  Data 

1.  Manufacturer. 

2.  Tube  identification  number. 

3.  Forward-wave  or  backward-wave 
amplifier  (see  note  1) . 

4.  Total  number  of  anode  vanes. 

5.  Number  of  active  anode  vanes. 

6.  Transverse  magnetic  field  (T) , 
treated  as  uniform. 

7.  Beam  width  (mm)  in  the  magnetic- 
field  direction,  treated  as  constant. 

3.  Cathode  radius  (mm). 

9.  Anode  radius  (mm),  measured  from 
the  axis  to  the  vane  tips . 

10.  Displacement,  if  any,  of  cathode 
axis  from  geometrical  center  (mm) . 

11.  Distance  (mm)  around  the  anode  from 
RF  input  to  RF  output. 

12.  Length  (mm)  around  the  anode  of 
the  RF  circuit  sever  between  the 
RF  output  and  RF  input.  (Items  li 
and  12  should  add  to  equal  the 
anode  circumference.) 

13.  Fitch  (period)  of  the  RF  circuit 
measured  around  the  anode  surface  (mm) 

14.  Ratio  of  vane  spacing  to  pitch 
(period)  on  the  anode  surface. 

15.  Drive  frequency  (GHz). 

a.  Lower  end  of  operating  band. 

b.  Midband. 

c.  Upper  end  of  operating  band. 

(These  frequencies  are  used  in 
following  items . ) 


Raytheon 


QKS1842 


Backward-wave 


40 


0.4115 


14.8 


7.748 


9.708 


None 


52.4 


8.6 


1.52 


0.328 


■  3.«5. 


9.75 


10.0 
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as  a  function 


16.  Anode-cathode  voltage, 
of  frequency: 

Frequency  (GHz)  a. 

b. 

c . 

Anode-cathode  voltage.  a. 

b . 

c. 

17.  RF  drive  power  (W)  for  peak 
output,  tabulated  as  a  function 
of  frequency: 

Frequency  (GHz)  a. 

b. 

c . 

RF  drive  (W) .  a. 

b. 

c. 


Notice  that  a  separata  run  must 
be  made  for  each  frequency  specified. 

18.  Cold-circuit  phase  delay  (degrees) 
per  vane  in  the  direction  of  power 
flow : 

Frequency  (GHz)  a. 

b . 

c . 


Phase  delay  (degrees)  . 


a. 


b. 

c . 


9.5 


9. 

75 

10 

.0 

24 

,300 

25 

,500 

26 

,850 

9.5 

9.75 

10.0 

5,500 

5,500 

5,500 


9.5 

d. 

9.8 

9.6 

e . 

9.9 

9.7 

• 

10.0 

41.0 

d. 

48.5 

43.5 

e  . 

51.0 

46.0 

f . 

53.5 

131 


*.*  V  *•* 

i 

1 

B 

19. 

Interaction  impedance  (2)  as  a 
function  of  frequency  (use  the 
definition  in  note  3  or  supply 
an  alternative  definition  in 

Part  III) : 

Frequency  (GHz) 

a . 

9.5 

d.  9.8 

U- 

b . 

9.6 

e.  9.9 

■ 

c . 

9.7 

f.  io.o 

Impedance  (2) . 

a . 

6.0 

d.  7.8 

b. 

6.6 

e.  8-5 

1 

c . 

7.2 

f.  9.2 

20. 

Cold-circuit  RF  power  attenua¬ 
tion  (dB)  per  anode  pitch  as  a 
function  of  frequency: 

►  -* 

Frequency  (GHz) 

a. 

9.5 

d.  9.8 

•■-. 

b . 

9.6 

e.  9.9 

'  _ 

:  ■ 

5 

c . 

9.7 

f.  10.0 

Attenuation  (dB) . 

a . 

.0187 

d.  .0187 

»•** 

b. 

.0187 

e.  .0187 

L  ‘ 

c . 

.0187 

f.  .0187 

s 

21. 

Type  of  cathode  emitter  surface: 

a.  Thermal  emitter  (yes/no) . 

b.  Secondary  emitter  (yes/no) . 

'-M 

c.  Both  thermal  and  secondary 
emitter  (yes/no) . 

Yes 

■ 

For  thermally  emitting  cathode  (a) 

• 

« 

► 

(1)  Current  density  that  would 
exist  under  temperature- 
limited  conditions  and  with 
no  transverse  magnetic  field 
(A/m2 ) . 

Estimated  as  less  than 

105  of  operating  value. 

« 

» 

4 
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(2)  Average  energy  (eV)  of  thermal 
electrons  on  emission. 


Not  known 


For  secondary  emitter  (b)  : 

(1)  Secondary-emission  coefficient 

(may  be  fractional)  as  a  function 
of  primary  impact  energy  (eV) .  A 
single  constant  value  will  suffice 
if  the  detailed  variation  is  not 
known . 


Primary  Impact 
Enerqv 


Secondary -Emission 
Coefficient 


1.86 


2.36 


2.58 


2.63 


2.63 


1,000 


2.63 


1,200 


2.58 


1,400 


2.53 


Xi ,  Performance  Data  for  Comparison  of  Computation  and  Measurement 

Please  indicate  whether  the  value  supplied  is  a  result  of 
measurement  or  an  analytical  estimate.  If  no  value  is  available, 
please  write  "not  known." 

1.  Drive  frequency  (GHz).  a.  9 . 5 _ d.  9.8 - 


10.0 


Anode  current  (A) . 


15.3  d.  17.1 
16.7  e.  17.1 


17.1 


17.1 


RF  output  power  (W) . 


Efficiency  (see  note  4) . 


159,000  d.  169,000 
184,000  e.  158,000 
181,000  f.  153,000 


40.9 

42.5 


36.6 

33.1 


40.4 


31.7 


Measured  hot  phase  delay 
(degrees)  per  cavity  in 
direction  of  power  :low. 


Not  known 


Hot  RF  phase  delay  (degrees) 
between  input  and  output. 


a .  Not  known, 


Total  anode  power  dissipation 
(W) ,  from  both  RF  attenuation 
and  beam  interception. 


a .  158,000  d.  195,000 

b.  173,000  a .  205,000 

c.  178,000  f.  218,000 
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8. 


Total  cathode  power  dissipation 
(W)  due  to  backbombardment. 


a. 


58.000  d .  62,600 


58,000 

e .  65,800 

58,000 

f.  68,100 

III .  Space  for  Additional  Comments 

The  above  values  were  supplied  in  September,  1976.  The 
following  additional  data  were  supplied  by  Raytheon. 

(1)  The  low  frequency  cutoff  of  the  vane  circuit  occurs 
at  7.88  GHz. 


(2)  The  following  power  measurements  at  9.75  GHz  were 
supplied. 


Date 

Supplied 

Tube 

Number 

Anode 

Voltage 

(V) 

Magnetic 

Field 

(T) 

Anode 

Current 

(A) 

RF 

Output 

(W) 

1/13/77 

8A 

24,300 

0.4114 

17.1 

168,000 

1/13/77 

8A 

22,000 

0.3788 

15.0 

124,000 

4/10/78 

3C 

25,000 

Between 

17.1 

185,000 

0 . 4  and 

0.44 

4/10/78 

A10 

24,600 

Between 

12.9 

132,000 

(platinum 

0.3982 

cathode) 

and 

0.444 

(3) 

Modulator . 

The 

QKS1842  uses  a  line-type  pulser  with 

an  open-circuit 

voltage  approximately  2.2  times 

the  operating 

voltage,  a 

series 

resistance,  a  series 

inductance 

of  about  100 

yH,  and  in  parallel 

with  the 

CFA  a  leakage  capacitance  of  about  100  pF.  The 

CFA 

itself  has  a  capacitance  of  3.65 

pF.  The  pulse  length  is 

1.7  ys . 

(4) 

The  following  data  were  supplied  for  operation 

of  a  tube 

in  the  high-power  mode 

Tube:  QKS1842  X-Band  Test  Vehicle  No.  1. 

Date:  Feb.  21,  1978. 


135 


■ 

4 

%  ;• 

Duty  factor 

0.001 

□ 

Electromagnet  current 

0.64  A 

>V* 

Magnetic  field 

0.3775  T 

«  4  "• 

Frequency 

9.6  GHz 

i 

RF  Drive  Power 

19.0  kW 

.  : 

Cathode:  gold  with 

magnesium  oxide 

,  as  a  cold  secondary 

emitter. 

a 

Measured  results: 

-  ; 

Anode  Current  (A) 

Anode  voltage 

(kV)  RF  Output  Power 

6.0 

18.0 

72.0 

10.0 

19.0 

128.0 

r  ' 

16.0 

20.2 

200.0 

20.0 

20.6 

248.0 

9 

25.0 

21.0 

308.0 

30.0 

21.2 

352.0 

,  _ 

41.0 

21.7 

453.0 

49.0 

21.8 

530.0 

P 

The  high  RF  drive  power 

is 

necessary  to 

suppress  oscillation. 

« 


L 


» 


Notes 


In  a  forward- wave  amplifier,  the  beam  and  the  RF  power  flow 
in  the  same  direction;  in  a  backward- wave  amplifier,  the  beam 
moves  around  the  tube  from  output  to  input  in  the  opposite 
direction  to  the  power  flow.  In  both  tubes,  the  average  beam 
velocity  and  the  phase  velocity  of  the  interacting  space  har¬ 
monic  are  in  synchronism  when  maximum  RF  interaction  occurs. 

The  cold  RF  phase  velocity  v  is  given  by 


where  f  (He)  is  the  drive  frequency  and  3q  (m  )  is  the  phase 

constant  for  the  interacting  space  harmonic,  such  that  the 
phase  delay  per  pitch  is  (3 q  *  pitch  length)  radians. 

The  RF  interaction  impedance  Z^  is  defined  as  the  expression 

IE  !  2 
1  \RF  1 

"0  235P 

where  is  the  peak  RF  electric  field  (V/m)  parallel  to  t 

anode  surface,  evaluated  at  the  anode  surface  for  the  ir.terac 

ing  space  harmonic, 3^ (m  )  is  the  cold-circuit  phase  constant 

for  the  same  space  harmonic  and  Pfiv)  is  the  total  power  flow 
parallel  to  the  anode  (in  all  space  harmonics  and  integrated 
over  the  entire  cross  section  of  the  interaction  region) . 

The  efficiency  is  defined  as  the  expression 


_ (RF  Output  -  RF  Drive) _ 

(Anode  Current) * ^Anode-Sole  Voltage) 


*  100) 


The  chase  velocitv  is  given  bv  one  of  the  following  ex: 


■  s  s  1 0 


for  the  forward -wave 
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DATA  FOR  THE  VARIAN  SFD-261  CFA 


Data  for  two  versions  of  the  SFD-261  are  given:  one 
instrumented  tube  with  a  uniform  anode-cathode  spacing,  and 
the  production  tube  which  has  a  non-circular  cathode. 

Normalizing  parameters  Pq,  fg,  D^,  are  used  to  allow 

an  UNCLASSIFIED  publication.  The  values  are  classified  CONFIDENTIAL. 


Basic  Tube  Data 


r« 


1 

X  • 

Manufacturer . 

Varian 

2. 

Tube  identification  number. 

SFD-261 

3. 

Forward-wave  or  backward-wave 
amplifier  (see  note  1)  . 

Forward-wave 

A 

• 

Total  number  of  anode  vanes. 

Not  known . 

5. 

Number  of  active  anode  vanes. 

62 

6  . 

Transverse  magnetic  field  (T) , 
treated  as  uniform. 

.  2530 

7 . 

Beam  width  (mm)  in  the  magnetic- 
field  direction,  treated  as  constant. 

400  D] 

3 . 

Cathode  radius  (mm) . 

735  Di 

9  . 

Anode  radius  (mm) ,  measured  from 
the  axis  to  the  vane  tips . 

825  Dt 

10  . 

Displacement,  if  any,  of  cathode 
axis  from  geometrical  center  (mm) . 

Zero  (instrumented 
tube  only) 

11. 

Distance  (mm)  around  the  anode  from 

RF  input  to  RF  output. 

4.52  D2 

12. 

Length  (mm)  around  the  anode  of 
the  RF  circuit  sever  between  the 

RF  outpuc  and  PF  inpur. 

0.66  D? 

13 . 

Pitch  (period)  of  the  RF  circuic 
measured  around  the  anode  surface  (mm) 

.  0.072  D2 

14. 

Ratio  of  vane  spacing  to  pitch 
(period)  on  the  anode  surface. 

0.5 

15  . 

Drive  frequency  (GHc) . 

a.  Lower  end  of  operating  band. 

f0  +  0.1 

b.  Midband. 

fn  +  °-3 

c.  Upper  er.d  of  operating  band. 

f0  -  0.5 

(These  frequencies  are  used  in 
following  items.) 
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Instrumented  Production 

16.  Anode-cathode  voltage,  as  a  function  Tube  Tube 

of  frequency: 

Frequency  (GHz)  a.  _ _ 

b .  _ 

c .  _ 

Anode-cathode  voltage.  a.  1.325  _ 1.33  Vi 

b.  1.3  Vx _ 1.304  VL 

c.  1.18  VL _ 1.27  Vj 

17.  RF  drive  power  (W)  for  peak 

output,  tabulated  as  a  function  Instrumented  Production 

of  frequency:  Tube  Tube 


Frequency  (GHz) 


c . 


RF  drive  (W)  . 

a . 

0.08 

P0 

b . 

0.08 

P0 

c . 

0.08 

P0 

Notice  that  a  separate  run  must 
be  made  for  each  frequency  specified. 


13.  Cold-circuit  phase  delay  (degrees) 

per  vane  in  the  direction  of  power  Instrumented 
clow:  Tube 

Frequency  (GHz)  a.  _ 

b .  _ 

c .  _ 

Phase  delay  (degrees) .  a.  77.99 _ 

b.  90.99 _ 

c.  104.49 


Production 

Tube 
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9.  Interaction  impedance  (.3)  as  a 
function  of  frequency  (use  the 
definition  in  note  3  or  supply 
an  alternative  definition  in 
Part  III) : 

Frequency  (GHz)  a. 

b . 


c . 


Impedance  (2).  a. 

b.  45 


c.  30 


Instrumented  Production 
TuLe  Tube 


Cold-circuit  RF  power  attenua¬ 
tion  (dB)  per  pitch  as  a  function  Instrumented  Production 

of  frequency:  Tube  Tube 


Frequency  (GHz)  a. 

b. 

c . 

Attenuation  (aB) .  a. 

b . 

c . 

Type  of  cathode  emitter  surface: 

a.  Thermal  emitter  (yes/no) . 

b.  Secondary  emitter  (yes/no). 

c.  3oth  thermal  and  secondary 
emitter  (yes/no). 

For  thermally  emitting  cathode  (a) : 


.036 

.0252 

.036 

.0252 

.036 

.0252 

No 

Yes 

(Beryllium  oxide 

No 

(1)  Current  density  that  would 
exist  under  cemperature- 
limited  conditions  and  with 
no  transverse  magnetic  field 

/  !!  /  rr  *- 


Not  known . 


Performance  Data  for  Comparison  of  Computation  and  Measursmen 


Please  indicate  whether  the  value  supplied  is  a  result  of 

measurement,  or  an  analytical  estimate.  If  no  value  is  available, 

please  write  "not  known."  Uniform  Off-Center 

r  Instrumented  Cathode 

1.  Drive  frequency  (GHz).  a.  ^0+0  *  ^ _ 


Anode  current  (A) . 
Measured 

RF  output  power  (W) . 
Measured 


Measured  hot  phase  delay 
(degrees)  per  cavity  in 
direction  of  power  flow. 


Hot  RF  phase  delay  (degrees) 
between  inout  ar.d  outout. 


Cathode 


Total  anode  power  dissipation 
(W) ,  from  both  RF  attenuation 
and  beam  in te reaction . 


a . 

r — 1 
. 

o 

+ 

o 

<4-4 

b . 

V0-3 

c . 

V0-5 

a . 

21.4 

22 

b . 

22 

22 

c . 

8.4 

22 

a . 

1.6  P. 

1.2  P1 

b . 

1.65  Px 

1.51  Px 

c  . 

0.62  P 

1.57  p 

a . 

Not  known. 

b . 

c . 

a . 

Not  known. 

b . 

a  . 

Not  known _ 

Estimated  as  0.05xdirect  input 
7.  Total  cathode  power  dissipation  power. 

(W)  due  to  backbombardment .  a.  _ d. _ 

b .  _ e_. _ 

c .  f . 


Ill .  Space  for  Additional  Comments 

All  the  data  summarized  above  are  taken  from  the  report  of 
Varian  Associates,  Inc.* 

The  variation  of  anode-cathode  spacing  in  the  production 
tube  is  tabulated  on  p.  154  of  that  report. 

Experimental  Accuracy 

Magnetic  field:  +  3% 

Anode  voltage:  +  1% 

Phase  delay  per  pitch:  ±  2% 

Modulator 

The  SFD-261  uses  a  hard-tube  modulator  with  an  open-circuit 
voltage  about  1.12  times  the  operating  voltage,  and  a  series 
resistance.  The  pulse  length  is  between  8  ys  and  12  ys.  Except 
for  the  starting  behavior,  the  tube  operation  is  observed  to  be 
insensitive  to  the  modulator,  and  open-circuit  voltages  up  to 
1.38  times  the  operating  voltage  have  been  used. 

The  CFA  capacitance  is  4.88  pF;  any  additional  stray  capaci¬ 
tance  has  not  been  measured. 


•k 

H.L.  McDowell,  "CFA  Design  Improvement  Program,  Volume  II-- 
Computer  Modeling  Studies,"  Final  Report,  Contract  No.  N00173-75- 
C-1294.  Prepared  for  U.S.  Navy  Ocean  Systems  Center,  San  Diego, 
California,  by  Varian  Associates,  Inc.,  Beverly,  Massachusetts; 

2  June  1978. 
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APPENDIX  D 

DATA  FOR  THE  RAYTHEON  QKS1319  CFA 

The  data  presented  here  were  measured  by  Raytheon  Company 
under  Subcontract  No.  009439  from  Harris  SAI,  Inc.  (then 
Shared  Applications,  Inc.).  The  report  provided  by  Raytheon 
is  included  here. 


!'« 
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1.0  PURPOSE 

A  computer  model  and  program  for  distributed  emission  crossed  field 
amplifiers  (DECFA)  that  is  now  under  development  at  Shared  Applications,  Inc. 
(SAI)  is  expected  to  improve  the  techniques  for  designing  these  devices.  Eval¬ 
uations  of  this  program  would  be  enhanced  by  using  design  and  performance 
data  from  an  existing  CFA.  The  purpose  of  this  task  was  to  provide  to 
SAI  certain  technical  information  about  one  of  Raytheon  Company*  s  CFAs,  said 
information  to  be  used  solely  for  the  evaluation  of  the  computer  model. 

2.0  WORKS  TA  TME  N  T 

The  CFA  information  requirements  of  SAI  were  described  in  the  Work 
Statement  delivered  to  Raytheon  as  part  of  this  subcontract  and  they  are  listed  as 
Appendix  I  to  this  final  report.  Data  on  one  tube  type,  the  QKS1319,  was 
provided.  The  QKS1319  is  a  forward-wave  device  whose  development  was  com¬ 
pleted  a  few  years  ago  and  which  is  now  being  used  in  a  production  version  of 
a  radar  system. 

A  photograph  of  the  distributed  emission  crossed  field  amplifier,  QKS13  19, 
is  depicted  in  Figure  1. 

2.  1  Work  Performed  -  General 

The  following  paragraphs  describe  the  work  performed  in  fulfillment 
of  each  item  in  the  work  statement. 

Some  of  the  information  requested  already  existed  and  a  simple 
compilation  was  all  chat  was  necessary.  Other  data  items  required  engineering 
analysis  and  computations.  The  balance  of  the  information  was  obtained 
through  both  iow  level  microwave  rr.easuremer.es  and  full  power  operation. 
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2.  1.  1  SAI  Data  Transmittal  Form 


The  Data  Transmittal  F orm  provided  by  SAI  was  completed 
to  the  extent  possible  using  existing  information  on  dimensions,  cathode 
properties,  circuit  rf  properties,  and  operating  parameters  and  data  under 
standard  conditions. 


report. 


The  form,  as  completed,  is  included  as  Appendix  II  to  this 


2.1.2  RF  Circuit  Interaction  Impedance  Characteristics 


Low  level  microwave  measurements  were  made  with  an 
existing,  typical,  circuit  and  computations  wer  =  made  to  determine  the  circuit 
impedance  characteristics.  These  appear  as  part  of  the  Data  Transmittal  Form, 
Appendix  II. 

2.1.3  Dispersion  Curve 


Low  level  microwave  measurements  were  made  on  the  cold 
rf  circuit  and  computations  were  made  to  determine  the  phase  shift  per  cavity 
as  a  function  of  frequency. 

The  dispersion  curve,  showing  lower  and  upper  cutoff 
frequencies,  is  presented  as  Figure  2. 

A  tabulation  of  numerical  values  of  the  measured  points  is 
included  as  Appendix  III  of  this  report. 
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2.1.4  Full  Power  Operation 


QKS1319,  Serial  No.  24,  was  installed  in  the  high-power 
test  equipment  with  an  electromagnet  and  measurements  were  made  of: 

•  The  minimum  rf  drive  power  for  starting 

•  The  upper  current  mode  boundary 

•  The  power  dissipated  on  the  anode  and  sole  (cathode) 

a.  Minimum  drive  power  for  starting  was  measured  by 
setting  tube  input  conditions  for  normal  operation  and 
then  reducing  the  input  drive  power  to  the  level  where 
the  tube  failed  to  amplify.  As  required,  this  test  was 
performed  at  three  rf  frequencies  and  three  values  of 
magnetic  field. 


« 


b. 


These  data  are  presented  graphically  in  Figure  3  and  a 
tabulation  of  the  numerical  values  for  the  measured 
points  appear  in  Appendix  IV. 

The  upper  current  mode  boundary  was  measured  by 
setting  tube  input  conditions  for  normal  operation  and 
then  raising  the  supply  voitage  to  the  current  level  where 
the  tube  failed  to  amplify. 

As  required,  this  test  was  performed  at  three  rf  frequencies 
and  three  values  of  magnetic  field. 

These  data  are  presented  graphically  in  Figure  4  and  a 
tabulation  or  the  numerical  values  for  the  measured 
points  appear  in  Appendix  IV. 
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c.  The  powet  dissipated  in  the  anode  and  in  the  sole  (cathode) 
was  measured  calorimetrically  using  a  thermocouple  in 
each  of  the  two  cooling  circuits. 

These  data  are  taken  at  each  of  the  test  points  previously- 
described  except  during  test  for  minimum  rf  drive  for 
starting,  and  are  presented  graphically  in  Figure  3  and  4. 
A  tabulation  of  the  numerical  values  for  the  measured 
points  appears  in  Appendix  IV. 

2.2  Conference  at  Raytheon,  June  29,  1978 

After  the  above  design  and  data  information  was  completed,  a 
representive  of  SAI  mec  with  engineers  at  Raytheon  for  the  purpose  of  discussion 
and  analysis  of  measured  and  computed  results. 

The  data  was  judged  to  be  satisfactory  and  in  compliance  with  the 
requirements.  No  additional  action  is  required. 
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CFA  DATA  PROGRAM  FOR  COMPUTER  MODELING 

Subcontract  No.  00943  9 
Terms  and  Conditions 


Statement  of  Work 


Raytheon  shall  provide  SAI  design  and  performance  data  on  one  type 

of  forward-wave  distributed-emission  crossed- field  amplifier,  the 
QKS1319,  as  follows: 

1.  Basic  tube  data  as  requested  in  Section  I  of  the  attached  SAI 
Data  Transmittal  Form.  This  includes  information  on  tube 
dimensions,  cathode  properties,  RF  circuit  properties  and 
standard  operating  parameters. 

2.  RF  circuit  interaction  impedance  characteristics,  as  requested 
on  page  5  of  the  SAI  Data  Transmittal  Form.  If  necessary, 
experimental  measurements  will  be  made  to  determine  the 
circuit  impedance. 

3.  Measurements  of  the  minimum  RF  drive  power  for  starting 
the  QKS13  19  at  three  RF  frequencies  and  three  values  of 
magnetic  field  (nine  measurements  total).  It  should  be  stated 
whether  the  tube  operates  with  constant  RF  drive  and  pulsed 
anode  voltage,  with  a  constant  anode  voltage  applied  and 
pulsed  RF  drive  or  under  other  conditions  to  be  described. 

4.  Measurements  of  the  upper  current  mode  boundary  at  three 
RF  frequencies  and  three  values  of  magnetic  field  (nine  mea¬ 
surements  total). 
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Measurement  of  the  power  dissipated  on  the  anode  and  sole 
(cathode)  at  the  highest  RF  output  power  level  at  each  of  three 
frequencies . 

The  data  requested  by  tasks  3,  4,  and  5  is  to  be  presented 
g  raphically. 

A  dispersion  curve  for  the  cold  RF  circuit  showing  frequency 
plotted  against  phase  shift  per  cavity.  The  curve  should  show 
the  lower  and  upper  cutoff  frequencies,  if  any.  Numerical 
values  of  measured  points  should  be  included  so  that  SAI  may 
derive  suitable  equivalent  network  parameters. 

One  conference  at  Raytheon  is  to  be  held  between  SAI  and 
Raytheon  engineers  to  discuss  and  clarify  the  data  presented 
and  the  measurement  techniques  used. 


DATA  TRANSMITTAL  FORM  FOR  COMPUTER  ANALYSIS  OF 
THE  DISTRIBUTED-EMISSION  CROSSED-FIELD  AMPLIFIER 


JULY  1977 


Nil 


SAI 


SHARED  APPLICATIONS,  INC. 
ANN  ARBOR,  MICHIGAN  48104 
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b 


4 


4 


I.  Basic  Tube  Data 


1.  Manufacturer 

2.  Tube  identification  number 

3.  Forward-wave  or  backward-wave  amplifier 
(see  note  1) 

4.  Total  number  of  anode  vanes 

5.  Number  of  active  anode  vanes 

6.  Transverse  magnetic  field  (T) ,  treated 
as  uniform 

7.  Beam  width  (mm)  in  the  magnetic  field 
direction,  treated  as  constant 

8.  Cathode  radius  (mm) 

9.  Anode  radius  (mm),  measured  from  the 
axis  to  the  vane  tips 

10.  Displacement,  if  any,  of  cathode  axis 
from  geometrical  center  (mm) 

11.  Distance  (mm)  around  the  anode  from  RF 
input  to  RF  output 

12.  Length  (mm)  around  the  anode  of  the  RF 
circuit  sever  between  the  RF  output  and 
RF  input  (Items  11  and  12  should  add 
to  equal  the  anode  circumference) 

13.  Pitch  (period)  of  the  RF  circuit  mea¬ 
sured  around  the  anode  surface  (mm) 

14.  Ratio  of  vane  spacing  to  pitch  (period) 
on  the  anode  surface 

15.  Drive  frequency  (Grit) 

a.  Lower  end  of  operating  band 

b.  Midband 

c.  -Upper  end  of  operating  band 

(These  frequencies  are  used  in  follow¬ 
ing  items) 
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Raytheon 

QKS1319 


Forward  Wave 


52 


45 


1150  Gauss 


35.  56  mm. 


34.  86  mm 


3  9.  70  mm 


Zero 


210.  96  mm 


3  S .  2  5  mm 


4.  79  mm 

0.  424 


1.  250 


1.3  00 


1.350 


16.  Anode-cathode  voltage,  as  a  function  of 
frequency: 


Frequency  (GHz) 

a. 

1.250 

b. 

1.300 

c. 

1.350 

Anode-cathode  voltage 

a. 

10,  200 

b. 

10,  200 

c. 

10,  200 

RF  drive  power  (W)  for 

peak  output, 

tabulated  as  a  function 

of  frequency: 

Frequency  (GHz) 

a. 

1. 250 

b. 

1.300 

c. 

1.350 

RF  drive  (W) 

a. 

3500 

b . 

3500 

c. 

3500 

Notice  that  a  separate  run  must  be  made 
for  each  frequency  specified. 

18.  Cold-circuit  phase  delay  (degrees) 
per  vane  in  the  direction  of  power 
flow: 


Frequency  (GHz) 

a .  _ 

1. 250 

b.  _ 

1.  300 

c .  _ 

1.350 

Phase  delay  (degrees) 

a.  _ 

100° 

b.  _ 

113° 

c . 

1 3  2  c 
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19.  Interaction  impedance  (3)  as  a  function 
of  frequency  (use  the  definition  in 
note  3  or  supply  an  alternative  defini¬ 
tion  in  Part  III): 


Frequency  (GHz) 


Impedance  (£) 


a.  1. 250 

b.  1.300 

c.  1-350 

a.  146.00 


b.  84. 17 

c-  _ 5-5,91 


20.  Cold-circuit  RF  power  attenuation  (d3) 
per  vane  as  a  function  of  frequency: 

Frequency  (GHz)  a.  _ 1.  250 

b .  I.  300 

C.  1.350 

Attenuation  (d3)  a.  _ .  044 

b.  _ .  03  8 

C.  .035 


Type  of  cathode  emit 

ter  surface: 

a . 

Thermal  emitter 

(yes/nc) 

no 

b . 

Secondary  emitte 

r  (yes/no) 

yes  (beryllium) 

c. 

Both  thermal  and 

secondary  emitter 

(yes/no) 

no 

For 

thermally  emitting  cathode  (a) : 

(1) 

current  density 

that  would  exist 

under  temperature-limited  condi¬ 
tions  and  with  no  transverse 
magnetic  field  (A/m2)  _ X 


(2)  average  energy  (eV)  of  thermal 

electrons  on  emission  X 
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For  secondary  emitter  (b)  : 

(1)  Secondary-emission  coefficient  (may 
be  fractional)  as  a  function  of  pri¬ 
mary  impact  energy  (eV) .  A  single 
constant  value  will  suffice  if  the 
detailed  variation  is  not  known 

Primary  Impact  Secondary-Emission 

_ Energy _  _ Coefficient 

50  eV  1.  2 


(2)  Average  energy  (eV)  of  secondary 

electrons  on  emission  Not  known 
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APPENDIX  III 


N 

Frequency  (MHz) 

_N_ 

F reouencv  (MHz) 

5 

1005 

25 

1255 

6 

1021 

26 

1266 

i 

103  8 

27 

1277 

5 

1054 

28 

1289 

Q 

1070 

29 

1301 

10 

1085 

30 

1314 

1  1 

1099 

31 

1329 

12 

1113 

32 

1343 

13 

1126 

33 

1360 

U 

1139 

34 

1377 

1 5 

1150 

3  5 

1396 

1 6 

1  162 

3  6 

i  4 1 6 

17 

1172 

37 

1437 

IS 

1 133 

38 

1460 

19 

1193 

39 

1483 

2  0 

1203 

4  0 

1505 

21 

1213 

41 

153  0 

22 

1223 

42 

15  5  2 

23 

1235 

43 

1573 

24 

1244 

t  < 

15  33 
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APPENDIX  IV 


FULL  POWER  OPERATION  DATA  QKS13  19 
Performance  (§'  20  amps  (peak)  (5  .  01  du 


FI 

F2 

F3 

* 

Power  Output 

115 

111 

113 

B  Normal 

kW  (peak) 

95 

94 

93 

B  Medium 

78 

78 

80 

B  Low 

Prise  Voltage 

10.2 

10.2 

10.  2 

B  Normal 

tkV) 

9.43 

9.43 

9.43 

B  Medium 

S.  98 

8.  82 

3.  68 

B  Low  * 

Anode  Dissipation 

1125 

1230 

1100 

3  Normal 

(watts  i 

995 

1177 

1177 

3  Medium 

1220 

1 123 

1123 

3  Low 

Cathode  Dissipation 

104 

87 

90 

B  No  rmal 

(watts) 

92 

83 

78 

B  Medium 

112 

87 

87 

3  L  ow 

Min  RF  Drive 

175 

525 

23  60 

3  Normal 

for  star; 

262 

612 

2100 

3  Medium 

i  watts ) 

280 

5  77 

15  75 

B  Low 

*  _ 

B  Normal 

=  1150  Gaus  s 

3  Medium 

=  1050  Gauss 

B  Low 

=  950  Gaus s 
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FULL  POWER  OPERATION  DATA  OKS  13  19 


B.  Performance  above  20  amps  (peak)  @  .  01  du 


FI 

F2 

F3 

* 

UCMB 

30  (emission) 

>3  0 

3  0  (emiss  ion) 

B  Normal 

ib 

>30 

>30 

25  (M.B.) 

B  Medium 

22  (M.B.) 

20  (M.B.) 

25  (emission) 

B  Low 

Power  Output 

150 

140 

133 

B  Normal 

k',V  (peak) 

106 

114 

111 

3  Medium 

73 

78 

SO 

3  Low 

Pulse  Voltage 

11.0 

10.6 

10.9 

3  Normal 

(kV) 

o 

o 

—4 

9.  9 

10.  0 

3  Medium 

9.0 

3.3 

9.0 

3  Low 

Power 

Dis  s  ipation 

2140 

1926 

2  i40 

3  Normal 

( anode  1 

1534 

2  03  3 

2  097 

3  Medium 

watts 

1273 

1123 

1123 

3  L  ow 

Cathode 

Dis  s  ipacior. 

1  "3 

130 

113 

3  Normal 

'.watts  i 

107 

104 

93 

B  Medium 

In  addition  to  the  information  given  in  the  Raytheon 
report,  the  following  additional  tube  data  were  supplied  at 
the  June  29,  1978  meeting  at  Raytheon. 

Severed  Cathode 

The  cathode  is  in  two  regions:  an  emitting  region  and  a 
non-emitting  region.  The  emitting  region  covers  the  angular 
segment  -4.065°  _<  0  _<  308.69°,  where  0  =  0  is  defined  as  the  center 
of  the  input  anode  vane.  The  non-emitting  region  covers  the 
segment  312.86°  £  0  £  351.78°.  Two  4.17°  gaps  separate  the  emitting 
and  non-emitting  regions.  The  non-emitter  is  recessed  with  a 
radius  of  33.7  mm,  whereas  the  emitter  radius  is  34.86  mm. 

The  following  information  was  supplied  by  telephone  in  January 

1980. 

Modulator 

The  QKS1319  uses  a  hard-tube  modulator  with  an  open-circuit 
voltage  of  15  kv  and  a  series  regulator  tube  with  a  voltage  drop 
of  about  4  kV.  The  effective  open-circuit  voltage  seen  by  the 
CFA  without  RF  drive  is  between  10.5  and  11  kV  depending  on  the 
actual  equipment  used.  The  pulse  length  is  40  ps  (52,000  RF 
periods)  at  1%  duty  cycle. 

Measured  Power  Balance 

The  measured  power  balance  shows  a  discrepancy  of  about 
17  percent.  The  errors  are  probably  in  the  anode  dissipation; 

6  percent  is  a  more  typical  experimental  error. 

Magnetic-Field  Measurements 

The  magnetic  field  is  measured  with  a  gaussmeter  in  the 
absence  of  the  tube  both  with  an  ALNICO  permanent  magnet  and  with 
an  electromagnet  using  the  same  steel  pole  pieces.  The  electromagnet 
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current  values  are  tabulated  for  the  measured  fields.  These 
currents  are  then  used  to  estimate  the  field  when  the  tube 
is  inserted  in  the  magnet.  The  field  values  supplied  are 
reproducible  measurements  and  are  believed  to  be  accurate. 
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APPENDIX  E 


DATA  FOR  THE  RAYTHEON  QKS1705  CFA 


I .  3asic  Tube  Data 

1.  Manufacturer. 

2.  Tube  identification  number. 

3.  Forward-wave  or  backward-wave 
amplifier  (see  note  1)  . 

4.  Total  number  of  anode  vanes. 

5.  Number  of  active  anode  vanes. 

6.  Transverse  magnetic  field  (T) , 
treated  as  uniform. 

7.  Beam  width  (mm)  in  the  magnetic- 
field  direction,  treated  as  constant. 

8.  Cathode  radius  (mm). 

9.  Anode  radius  (mm),  measured  from 
the  axis  to  the  vane  tips. 

10.  Displacement,  if  any,  of  cathode 
axis  from  geometrical  center  (mm) . 

11.  Distance  (mm)  around  the  anode  from 
RF  input  to  R?  output. 

12.  Length  (mm)  around  the  anode  of 
the  RF  circuit  sever  between  the 
RF  output  and  RF  input.  (Items  11 
ana  12  should  add  to  equal  the 
anode  circumference.) 

13.  Pitch  (period)  of  the  RF  circuit 
measured  around  the  anode  surface  (mm) 

14.  Ratio  of  vane  spacing  to  pitch 
(period)  on  the  anode  surface. 

15.  Drive  frequency  (GHz). 

a.  Lower  end  of  operating  band. 

b.  Midband. 

c.  Upper  end  of  operating  band. 

(These  frequencies  are  used  in 
following  items . ) 


Raytheon 

QKS1705 


Backward 


Not  known. 


35 


0.36 


14.8 


SLJIQ 

11.94 


zero 


65.5 


9.5 


.1.829 


0.404 


9.5 


9.75 


10.0 
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16.  Anode-cathode  voltage,  as  a  function 
of  frequency: 


Frequency  (GHz) 


Anode-cathode  voltage. 


17.  RF  drive  power  (W)  for  peak 

output,  tabulated  as  a  function 
of  frequency: 

Frequency  (GHz) 


RF  drive  (W)  . 


a.  9.5 

b.  9.75 

c.  10.0 

a.  31,800 

b.  33,000 

c.  34,000 


a. 

b . 

c . 

a . 

b . 


9.5  d.  9.8 

9.6  e.  9.9 

9.7  f.  10.0 
30,000  d.  30,000 
30,000  e.  30,000 
30,000  f.  30,000 


Notice  that  a  separate  run  must 
be  made  for  each  frecuencv  soecified. 


18.  Cold-circuit  phase  delay  (degrees) 
per  vane  in  the  direction  of  power 
flow: 


Frequency  (GHz)  a. 

9.5 

a . 

9.8 

b. 

9.6 

e  . 

9.9 

c . 

9.7 

£ 

10.0 

Phase  delay  (degrees) .  a. 

43 

d. 

49 

b . 

45 

e  . 

51 

C  • 

47 

- 

53 
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LC 


B 


19. 


Interaction  impedance  ((2)  as  a 
function  of  frequency  (use  the 
definition  in  note  3  or  supply 
an  alternative  definition  in 
Part  III) : 

Frequency  (GKz ) 


Impedance  (£2) 


a. 

in 

• 

<7> 

d. 

9.8 

b . 

9.6 

e . 

9.9 

c . 

9.7 

f . 

10.0 

a . 

7.4 

d. 

9.2 

.Q 

8.0 

e . 

9.8 

c . 

8.6 

'  f . 

10.4 
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Cold-circuit  RF  power  attenua¬ 
tion  (dB)  per  anode  pitch  as  a 
function  of  frequency: 


* 

Frequency  (GHz) 

a. 

9.5 

d. 

9.8 

►  * 

b . 

9.6 

e . 

9.9 

c . 

9.7 

f  # 

10.0 

Atts 

nuation  (dB) . 

a . 

0.018 

d. 

0.018 

b . 

0.018 

0.018 

c . 

0.018 

• 

0.018 

21.  Type 

of  cathode  emitter  surface: 

r . 

a . 

Thermal  emitter  (yes/no) . 

»* 

b. 

Secondary  emitter  (yes/no) . 

• 

- 

c . 

Both  thermal  and  secondary 
emitter  (ves/no) . 

Yes 

• 

For 

thermally  emitting  cathode  (a) 

: 

4 

k- 

- 

(i) 

Current  density  that  would 
exist  under  temperature- 
limited  conditions  and  with 
no  transverse  magnetic  field 
(A/m2  }  . 

Not  known. 
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(2)  Average  energy  (eV)  of  thermal 

electrons  on  emission.  Not  known 

For  secondary  emitter  (b)  : 

(1)  Secondary-emission  coefficient  - 
(may  be  fractional)  as  a  function 
of  primary  impact  energy  (eV) .  A 
single  constant  value  will  suffice 
if  the  detailed  variation  is  not 
known . 


Average  energy  (eV)  of  secondary 
electrons  on  emission. 


Not  known 


t;  II.  Performance  Data  for  Comparison  of  Computation  and  Measurement 

I  ‘ 


Please  indicate  whether  the  value 

supplied  is 

a  result  of 

c 

measurement  or  an  analytical  estimate. 

please  write  "not  known." 

If 

no  value 

is  available 

»’ 

1. 

Drive  frequency  (GHz) . 

a . 

9.5 

d.  9.8 

• 

s 

b . 

9.6 

e.  9.9 

■ 

c . 

9.7 

f.  10.0 

» 

2. 

Anode  current  (A) . 

a . 

34 

d.  34 

1 

Measured 

b . 

34 

e.  34 

c . 

34 

f.  34 

3. 

RF  output:  power  (W)  . 

a . 

546,000 

d.  521,000 

Measured 

b . 

502,000 

e<  519,000 

t  * 

c . 

504,000 

fm  504,000 

4  . 

Efficiency  (see  note  4) . 

a. 

45.6 

d.  4*-7 

t 

Measured 

b . 

40.4 

e.  40*7 

40.5 

f  39.3 

• 

S 

5  . 

Measured  hot  phase  delay 
(decrees)  per  cavity  in 
direction  of  power  flow. 

a . 

Not  known. 

u 

•J  . 

. 

C  . 

'  < 

6  . 

Hoc  RF  phase  delay  (degrees) 
between  input  and  output. 

a . 

Not  known. 

• 

b . 

« 

c . 

7  . 

Total  anode  power  dissipation 
(W) ,  from  booh  RF  attenuation 
and  beam  interception. 

a  . 

488,000 

d. 552, 000 

« 

Estimated 

b . 

562,000 

e  .  577,000 

r 

^  • 

552,000 

f .  594,000 
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Total  cathode  power  dissipation 
(W)  due  to  backbcmbardment . 

Estimated 


a. 

54.400 

d. 56. 450 

b. 

54,900 

e. 57. 300 

c . 

55,750 

f .57,800 

Space  for  Additional  Comments 
Calibration  of  electromagnet: 
Magnet  current  (A) 


Magnetic  field  (T) 


.3250 


.3750 


Duty  factor  0.001 

Operating  pulse  length  0.44  us 


APPENDIX  F 


DATA  FOR  THE  RAYTHEON  QKS1840  CATHODE-CIRCUIT  CFA 


Basic  Tube  Data 


Manufacturer . 


Raytheon 


Tube  identification  number. 


)KS1840 


Forward-wave  or  backward-wave 
amplifier  (see  note  1) . 

Total  number  of  anode  vanes. 

Number  of  active  anode  vanes. 


Forward-wave 


Transverse  magnetic  field  (T) , 
treated  as  uniform. 

Beam  width  (mm)  in  the  magnetic- 
field  direction,  treated  as  constant. 

Cathode  radius  (mm) . 

Anode  radius  (mm) ,  measured  from 
the  axis  to  the  vane  tips. 

Displacement,  if  any,  of  cathode 
axis  from  geometrical  center  (mm) . 

Distance  (mm)  around  the  anode  from 
RF  input  to  RF  output. 

Length  (mm)  around  the  anode  of 
the  RF  circuit  sever  between  the 
RF  output  and  RF  input.  (Items  11 
and  12  should  add  to  equal  the 
anode  circumference . ) 

Pitch  (period)  of  the  RF  circuit 
measured  around  the  anode  surface  (mm) 

Ratio  of  vane  spacing  to  pitch 
(period)  on  the  anode  surface. 

Drive  frequency  (GHz) . 

a.  Lower  end  of  operating  band. 

b.  Midband. 


0.3  (approx.) 


14.48 


20.89 


24.13 


129.5 


22.11 


.  3.175 


0 . 304 


?.ll 


3.292 


c.  Upper  end  of  operating  band. 

(These  frequencies  are  used  in 
following  items . ) 


3.54 
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Interaction  impedance  (ft)  of  anode 
circuit  as  a  function  of  frequency 
(use  the  definition  in  note  3) : 


Frequency  (GHz) 


Impedance  (P.)  . 


a . 

b. 

c . 

a . 

b . 


c . 

Cold-circuit  RF  power  attenua¬ 
tion  (aB)  per  vane  as  a  function 
of  frequency: 

Frequency  (GHz)  a. 

b. 


Attenuation  (dB) . 


c . 

a . 

b . 


Cathode  rui  ter  i  a  1 . 


Fcrward-wave  or  backward-wave 
cathode  circuit. 


Total  number  of  cathode  vanes. 
Numuer  of  active  carhode  vanes. 


3.1  _ 

3.3  _ 

3.5  _ 

135 _ 

75 _ 

40  _ 

3.1  _ 

3.3  _ 

3.5 

0.05 

0.05 

0.05 

Platinum  (secondary 
emitter) _ 

Forward-wave;  wave 

is  parallel  to _ 

anode-circuit  wave. 

48_ _ 

40 


2.743 


25.  Pitch  (period)  of  the  cathode 
RF  circuit  (nun)  . 

26.  Ratio  of  vane  spacing  to  pitch 
on  cathode  surface. 

(The  vane  centers  on  the  cathode  and 
anode  have  no  angular  separation.) 

27.  Cold-circuit  phase  delay  (degrees) 
per  pitch  of  cathode  circuit  in  the 
direction  of  cathode  power  flow: 

Frequency  (GHz)  a. 

b. 

c . 

Phase  delay  (degrees)  per  pitch.  a. 

b . 

c . 

28.  Interaction  impedance  (3)  for 
cathode  circuit  as  a  function 
of  frequency  (see  definition  in 
note  3) : 

Frequency  (GHz)  a. 

b. 

c . 

Impedance  (.'.).  a. 

b . 

c . 


0.417 


3,100 

3,300 

3,500 

98 _ 

129 

147 


3 , 100 
3,300 
3,500 
136 

75 _ 

4  3 
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MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BUREAU*OF  STANDARDS- 1 963-A 


Performance  Data  for  Comparison  of  Computation  and  Measurement 


Please  indicate  whether  the  value  supplied  is  a  result  of 
measurement  or  an  analytical  estimate.  If  no  value  is  available 
please  write  "not  known." 


1.  Drive  frequency  (GHz). 


a.  3.110 


Efficiency  (%) . 


Hot  RF  phase  delay  (degrees) 
between  input  and  output. 


b.  3.292 


c.  3.540 


2.  Anode  current  (A). 


a .  49 


c.  45 


RF  output  power  (W) . 


a.  312,000 

b.  337,000 

c.  250,000 


Measured  hot  phase  delay 
(degrees)  per  cavity  in 
direction  of  power  flow. 


a.  Not  known. 


7. 


Total  anode  power  dissipation 
(W) ,  from  both  RF  attenuation 
and  beam  interception . 


Not  known. 


APPENDIX  G 


MEASURED  SECONDARY-EMISSION  CHARACTERISTICS 
OF  CATHODE  SURFACES 

The  data  here  show  the  secondary-emission  coefficient  as 
a  function  of  primary  impact  energy  for  normal  incidence  for 
(1)  platinum,  (2)  cermet  ( tungsten-thoria) ,  (3)  beryllium 
oxide,  and  (4)  gold  and  magnesium  oxide  surfaces. 


I .  Platinum 

Source:  Raytheon  Company. 

The  values  are  plotted  in  Fig.  G.l  and  tabulated  below. 

Primary  Impact  Secondary-Emission 

Energy  (eV)  Coefficient 


Secondary- 

Emission 

Coefficient 


Figure  G.l  Measured  secondary-emission  coefficient  versus 
primary  impact  energy  for  platinum  and  cermet 
cathodes . 


L  8  5 


Cermet  ( Tuners ten-Thona) 
Source:  Raytheon  Company, 
See  Fig.  G.l. 


Primary  Impact 


Secondary-Emission 

Coefficient 


III.  Beryllium  Oxide 


Source:  Mr.  Richard  Thomas,  Naval  Research  Laboratory, 
Washington,  D.C. 

Fig.  G.2  shows  the  experimental  values  obtained  at  NRL. 

The  following  compares  the  experimental  values  with  those 
originally  used  in  the  Varian  and  Harris  SAI  computer  simulations 
of  the  SFD-261  CFA. 


Figure  G.3  Secondary-emission  coefficient  of  magnesium  oxide  and  gold  surface. 
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APPENDIX  H 


CIRCUIT  THEORY  FOR  THE  DECFA 

This  appendix  derives  the  equivalent  networks  for  the  RF 
circuit  and  modulator,  and  gives  the  details  of  the  calculation 
of  the  applied  electric  field  and  the  induced  currents. 

The  three  basic  approximations  are  that 

(1)  the  electric  fields  are  static, 

(2)  the  modulator  voltage  varies  slowly  over  one  RF  period,  and 

(3)  the  instantaneous  average  of  the  vane  voltages  is 
equal  to  the  modulator  voltage. 

These  approximations  enable  the  vane" and  modulator  networks 
to  be  solved  independently.  A  further  approximation  of  uniform 
azimuthal  field  between  the  vane  tips  allows  a  cylindrical 
solution  of  Poisson's  equation  to  be  used,  with  the  corresponding 
analytic  solution  of  Laplace's  equation  for  the  applied  field. 

A.  Separation  of  Electric  Fields  due  to  Applied  Voltages  and 
Space  Charge 

With  the  electrostatic  approximation,  the  total  field  at 

any  point  is  separable  into  an  external  part  due  to  applied 

direct  and  RF  voltages  and  a  space-charge  part.  Each  field  is 

derived  from  an  electrostatic  potential.  The  space  charge 

potential,  $  ,  satisfies  the  two-dimensional  Poisson's  equation 

sc 

and  is  zero  on  tne  cylindrical  surface  of  the  cathode  and  on 

the  vane  tips.  It  is  also  approximated  as  zero  between  the 

vane  tips.  The  external  potential  V  satisfies  Laplace's  equation 

and  the  boundary  conditions  of  applied  voltages  on  the  vanes. 

The  total  potential,  $  +  V,  satisfies  both  Poisson's  equation 

s  c 

and  the  boundary  conditions  and  is  therefore  the  unique  solution 
to  the  static  problem. 


Green's  Function 

Suppose  that  vane  n  of  the  CFA  has  potential  Vn  and 


that  all  other  vanes  have  potential  zero,  relative  to  the 
cathode.  The  potential  at  point  (r,  0)  in  the  anode-cathode 
region  is  then 


$  =  vn  g  (r,  e-en) 


where  9n  is  the  angular  coordinate  at  the  center  of  the  vane, 
and  G  is  a  Green's  function.  Make  the  approximation  that  the 
angular  electric  field  is  uniform  around  the  circumference  in 
the  gap  between  the  vanes.  This  is  a  useful  approximation  at 
the  anode,  but  it  would  not  be  sufficiently  accurate  for  an 
emitting  cathode  circuit  where  electrons  move  more  slowly. 

Define  the  following  notations: 

0B  =  one-half  the  angle  subtended  at  the  axis  by  the 
gap  between  adjacent  vanes. 

0G  =  one-half  the  pitch  angle  subtended  by  one  vane 
and  its  adjacent  gap, 

r  =  anode  radius, 

a 

rc  =  cathode  (sole)  radius. 

Then  the  Green's  function  is  given  by  the  following 
solution  of  Laplace's  equation 


G (r , 9)  = 


a0  ln 


.  ?  . ,!V>« .  fet  ,v  -■ 


cos  k9 


0 


77  ln(ra/rc) 


in. 


and 


2  sin  kQg  sin  k0Q 


irk20/ 


L  v  * 


A2kl 

V  J 


(H.  4) 


for  k  21  1. 

The  coefficients  a.^  and  a^  are  derived  by  Fourier  analysis 
of  the  voltage  around  the  period  0  <  0  <  2ir,  with  G(r  ,0)  =  1 
on  one  vane  and  G(r  ,0)  decreasing  linearly  to  zero  in  the 

3 

gap  on  each  side. 


The  corresponding  electric  field  Green's  functions  are 


Gr(r'0)  =  -  37  G(r,0) 


and 


Gg  (r, 0) 


1  3G(r,0) 

r  30 


(H. 


(H. 
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Gr (r,6) 
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l  r  J  \  r  I 
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(H.  n 


*  -  “  * . 


G0 (r,8) 


where 


■  I,  *- [fer  -  ca>j"k*1 


sin  k8 
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(9b  +  9g) 

*  ln(ra/rc) 


(H.9) 


2  sin  k0R  sin  k0„ 

f  =  - pr-2- - (H.  10) 

,k6Gra  [l-(rc/ra)2kJ 

The  field  (G^.,  Ga)  has  zero  divergence  and  curl  to  satisfy 
Maxwell's  equations  for  a  quasistatic  field. 

In  the  model,  100  terms  of  the  anode  suffice  up  to  half 
way  from  the  sole  to  the  anode.  At  greater  radii,  500  terms 
are  used.  The  fields  are  evaluated  on  a  41  x  41  mesh.  Because 
of  the  symmetry  in  9,  only  the  fields  for  positive  9  need  to 
be  stored.  The  mesh  covers  values  of  6  up  to  2.5  vanes  distant 
from  the  reference  vanes,  as  trials  show  that  the  field  components 
are  less  than  six  percent  of  their  maximum  values  outside  this 
range.  On  the  anode,  the  radial  field  is  extrapolated  linearly 
from  the  two  nearest  mesh  points  and  the  angular  field  is 
known  exactly.  The  local  field  at  position  (r,6)  is  calculated 
for  a  given  particle  by  area-weighting  the  contributions  from 
the  four  nearest  mesh  points,  treating  the  mesh  as  locally 
rectangular. 

C.  Applied  Voltage 

It  is  convenient  to  separate  the  charge-free  potential  into 
a  slowly  varying  applied  part  V_,„,  (t)  due  to  the  modulator  and 


the  RF  part  vRFn(fc)  on  each  vane.  With  this  approximation  of 
slow  time  variation,  is  the  time  average,  over  1  RF 

period,  of  the  voltage  on  any  of  the  vanes. 

The  potential  Vn  on  vane  n  is 


V  =  V  +  V 
n  CFA  vRFn 


(H. 11) 


The  term  V^Fn  is  obtained  from  a  solution  of  the  network 
equations,  and  is  the  voltage  across  the  admittance  Y_  in  the 

a 

equivalent  network.  Figure  H.l. 

Summing  over  all  the  vanes  gives  the  total  potential 


V  (r ,  0 )  =  V 


NVANES 

cfa  2  G<r'0-y 

n=l 


NVANES 

+  2  VRFnG<r'0-V  ,H-12’ 

n=l 


Using  the  actual  expression  for  G,  Equation  H.2,  gives 


V(r,0)  = 


VCFA  ln  (r/rJ  ^ES„ 

ln(ra/rc)  z  RFn 

n=l 


G  (r,  0  -  6n) 


(H. 13) 


showing  the  separation  into  direct  and  RF  fields.  In  the  program 
the  two  terms  of  Equation  H.13  are  computed  separately. 


D.  Electrode  Current 

The  current  flowing  from  an  electrode  (cathode,  a  vane,  or 
the  drift  section  electrode)  into  an  external  circuit  is 


I 


+• 


(H . 14 ) 


cathode 


j  ICFA 


(t) 


Portion  of  CFA  equivalent  network. 
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where  Q  is  the  total  charge  on  the  surface  of  the  electrode. 

(Note  that  the  current  I  does  not  include  the  surface  "skin" 

n  - 

current  that  is  produced  by  the  magnetic  field  tangential  to 
the  surface.) 

Using  Gauss's  theorem,  an  equivalent  expression  for  this 
current  is 


n 


"  "  dt  J£ 0 


E  .d  S 


(H . 15) 


in  terms  of  the  outward  normal  component  of  the  electric  field 
E  (all  parts:  space  charge,  direct  and  RF)  on  the  surface  S  of 
the  electrode.  Equation  H.15  can  be  used  in  the  present  model 
to  calculate  the  total  cathode  current  but  not  for  the  RF  current 
on  an  individual  vane,  since  the  detailed  configuration  of  the 
slow-wave  circuit  is  not  included  in  the  calculation  of  the 
field  E.  In  general  a  three-dimensional  solution  of  the  wave 
equation  would  be  needed.  Separating  E  into  its  applied  component 
and  space-charge  component  gives  two  parts  of  the  current:  that 
due  to  the  electrode  static  potential  Vn  relative  to  the  cathode, 
and  the  induced  current  due  to  space  charge  in  the  device. 

The  current  can  therefore  be  written 

N 

I  =  induced  current  -  )  C  .dv./dt  (H.16) 

j  =  l 

where  N  is  the  total  number  of  electrodes,  and  C  .  is  the 

nD 

coefficient  of  capacitance  between  electrodes  n  and  j. 

Expressions  for  the  induced  current  are  given  in  Section 
F  of  this  appendix. 


E. 


Approximate  Equivalent  Circuits 

The  present  program  computes  the  induced  current  directly 
and  incorporates  the  capacitances  in  an  external  equivalent 
circuit. 

1 .  RF  network 

For  the  RF  network.  Equation  H.16  gives  for  vane  n 

a 

I  =  induced  current  -  -r-  (t) 

n  at  ufa 


N 


-  2  c  .  4 


j-i 


C  V 

nj  dt  RFj 


(H. 17) 


The  intervane  capacitances  Cn^  for  n^j  are  not  calculated 

directly,  but  are  included  in  the  admittance  Y^  of  Figure  HI. 

Here  the  capacitances  between  non-adjacent  vanes  are  neglected.’ 

The  capacitance  C  between  vane  n  and  the  cathode  is  included 

nn 

in  the  admittance  Y_ . 

a 

The  network  admittances  are  combinations  of  inductance, 
capacitance,  and  resistance,  which  cannot  be  computed  from  a 
purely  electrostatic  field.  Instead  they  are  chosen  to  fit  the 
experimental  phase  velocity,  interaction  impedance  and  attenuation 
that  are  obtained  in  cold  tests. 


Finally  the  approximation  of  slowly  varying  modulator  voltage 
allows  the  term  including  t0  536  ne9lecte<^  in  Equation 

H.17.  The  RF  network  (Section  VIII)  is  then  treated  as  independent 
of  the  modulator. 


2 .  Modulator  circuit 

For  the  modulator  current,  the  last  term  of  Equation  H.17 
is  neglected,  with  the  approximation  that  the  RF  voltages  cancel 


when  they  are  summed  at  any  instant.  (With  Equation  H.15 
this  approximation  would  be  unnecessary,  however.) 

Hence  Equation  H.17  becomes,  summing  over  n. 


XCT(t) 


induced  current  -  C^;VCFA(t) 


(H. 18) 


as  in  the  equivalent  circuit  of  Section  IX. 

The  capacitance  C  is  calculated  by  treating  the  anode 
and  cathode  as  concentric  cylinders,  and  is  then  added  to  any 
external  leakage  capacitance  that  is  to  be  included. 

F.  Derivation  of  Induced  Current 

1 .  Computation  from  charges 

The  total  charge-free  potential  at  position 

NVANES 

2  VnG<r'9  -  V 

n=l 

of  the  contributions  of  the  NVANES  vanes. 
(r,6)  induces  a  charge  Qn  on  vane  n  given 

-q  G  (r ,  0  -  0n)  .  (H.  20 ) 

Let  the  charge  q  have  velocity  v.  Then  the  induced  current  on 
vane  n  is^'^ 

dQ 

In  =  ^  =  -  q  v  •  VG(r,  9  -  0n)  (H.21) 


V(r ,  9)  = 
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Then  a  charge  q  at 
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(r , 0)  is 
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Notice  that  the  actual  rate  of  charge  collection  on  any  vane, 
or  net  emission  from  the  cathode  averaged  over  a  finite  time 
interval,  should  not  be  included  in  Equation  H.21.  The  long¬ 
term  averages  of  the  anode  convection  current  and  the  induced 
current  summed  over  any  set  of  vanes  are  equal.  For  each 
vane  the  instantaneous  current  is  given  by  summing  Equation  H.21 
for  all  the  charges  q between  the  anode  and  the  cathode. 

The  local  electric  field  En(r,0)  due  to  vane  n  without 
space  charge  is  given  by 


En  =  "  Vn  VG(r'9> 


Hence  combining  Equations  H.21  and  H.22  gives 

v  •  E 


In  “  “  q  V_ 


n 


n 


(H.22) 


(H.23) 


and  also 


V  I 


n  n 


q  v 


(H.24) 


Equation  H.24  shows  that  the  induced  current  on  vane  n  due  to 

charge  q  gives  the  rate  of  working  on  the  charge  of  the  field 

E_ ,  where  the  field  E  due  to  vane  n  is  computed  as  if  the  charge 
n  25  n 

q  were  absent. 

For  example,  the  induced  current  leaving  the  cathode, 
due  to  a  single  charge  q  at  radius  r,  is 

ICFA(t)  q  (dt)  r  In (r  /r  )  (H.25) 

where  the  anode  and  cathode  are  treated  as  concentric  cylinders 

of  radii  r=  and  r  . 

a  c 


111! 


Computation  from  fields 


Applying  the  divergence  theorem  to  the  volume  t  bounded  by 
the  anode  vanes,  the  gaps  between  the  vanes,  and  the  cathode 
(surfaces)  and  of  unit  width  in  the  magnetic- fie Id  direction, 
gives 


V*sc  *  W  dt 


«sc  IS  « 


*sc72v  dT 


(H. 26) 


-  /v  ^  ds  -  /v 
s  L 


7  *sc  dT 


(H. 27) 


where  3/3n  represents  the  derivative  along  the  outward  normal. 

This  is  zero,  because  $sc  is  zero  on  S  and  V  satisfies  Laplace's 
equation. 

Now  use,  in  place  of  V,  the  Green's  function  G(r,0  -  0r) 
with  unit  potential  on  vane  n  only.  Replace  the  volume  integral 
in  Equation  H.27  by  a  sum  over  all  charges  q.  Then  the  expression 
that  results  is 


-£g  G(r'  9'  V  =  so  / 


(H. 28) 


where  S  is  now  the  surface  of  vane  n.  This  is  an  alternative 
n 

expression  for  the  charge  induced  on  vane  n.  The  corresponding 
induced  current  due  to  all  charges  q  is  then 


If  this  expression  can  be  computed  with  sufficient  accuracy 
it  will  provide  the  induced  current  I  on  each  vane  more 
efficiently  than  direct  summation  of  Equation  H.21  over  all 
the  charges.  Notice,  however,  that  using  Equation  H.29  instead 
of  the  Green's  function  requires  two  solutions  of  Poisson's 
equation  per  time  step  in  order  to  separate  the  space-charge 
field  from  the  total  field. 


APPENDIX  I 


THEORY  OF  THE  RF  NETWORK 

Throughout  this  section  the  approximation  is  made  that 
the  modulator  voltage  varies  slowly  relative  to  the  RF  drive 
voltage.  Then  the  RF  network  can  be  treated  independently  of 
the  modulator. 


A  general  form  of  low-frequency  equivalent  network  is 
shown  in  Figure  7  (Section  VIII).  Nodes  1,  2,  3,  4,  . .., 
represent  the  vanes  of  the  CFA  circuit  and  their  voltages 
produce  quasistatic  RF  fields  seen  by  the  electrons. 

The  circuit  calculations  have  three  parts: 

1.  Derivation  and  solution  of  the  network  equations 
relating  the  vane  voltages  and  the  driving  current, 

2.  Calculation  of  the  induced  current  from  the  beam 
motion ,  and 

3.  Choice  of  the  network  parameters  to  represent  an 
actual  circuit. 

A.  General  Network  Equations.  Number  the  nodes  around 
the  tube  in  the  counterclockwise  direction  as  1  to  NVANES , 
where  NVANES  is  the  total  number  of  active  vanes  in  the  tube. 
For  a  forward-wave  circuit,  RF  drive  is  applied  at  node  1; 
for  a  backward-wave  circuit  at  node  NVANES.  Thus,  in  both 
cases,  the  RF  wave  and  the  synchronous  electrons  are  defined 
to  move  in  the  counterclockwise  direction  (9  increasing) . 

Define  the  following  notation: 

C  =  capacitance  (F)  between  node  and  ground, 

a 

r  =  reciprocal  inductance  (K  ^)  between  node 

cl 

and  ground, 

G  =  conductance  (Q  ^)  between  node  and  ground 

cl 

C.  =  capacitance  (F)  between  adjacent  nodes. 


Ji  1 

cos 

(u)t 

1— 1 

CM 

C 

•H 

H 

cos 

(wt 

-  y2> 

in  NV 

cos 

( 'Jjt  - 

ynv* 

in  NV1  COS  (ut  "  YNV1) 


reciprocal  inductance  (H  )  between 
adjacent  nodes, 

conductance  (Q  between  adjacent  nodes, 

capacitance  (F)  between  alternate  nodes, 

reciprocal  inductance  (H-1)  between 
laternate  nodes. 

conductance  between  alternate  nodes, 

load  capacitance  at  first  and  last  nodes, 

reciprocal  of  load  inductance  at  first  and 
last  nodes, 

load  conductance  at  first  and  last  nodes, 

load  capacitance  at  second  and  penultimate 
nodes , 

reciprocal  of  load  inductance  at  second 
and  penultimate  nodes, 

load  conductance  at  second  and  penultimate 
nodes , 

driving  current  at  first  node, 

driving  current  at  second  node  (normally 
zero)  , 

driving  current  at  last  node  (used  for  a 
backward-wave  circuit) , 

driving  current  at  penultimate  node  (normally 
zero) 


ill 


v^Ctj  =  instantaneous  voltage  at  node  k  (1  =  k 
=  NVANES)  , 

Iv(t)  =  induced  current  at  node  k.  (Appendix  H) 


Applying  Kirchoff's  Laws  gives  an  expression  for  the 
voltage  at  an  intermediate  vane  (3  =  k  =  NVANES-2)  as 


Ca  dt*  '  Cb  dl:^  (Vk-l'2Vk+Vk-l)  '  Cc  3F^(Vk+2 ' 2Vk*Vk- 2] 


j  2 

-2-^rv.  .  - 2 1 


uv,  , 

*  Ga  TT  -  Gb  K(Vl'2Vk‘vk-l)  -  Gc 


'aVk  ‘  rb(Vk  +  r2Vk+Vk-l)  '  :c(Vk^*2Vk+Vk-^ 


Similar  equations  apply  at  the  end  nodes  with  the  appropriate 
admittance  terms  replaced  by  zero  since  the  connecting  network 
elements  are  absent  there.  In  matrix  form,  these  equations 
become 


£  JtJ  &  3t  -v 


d 1 

dt  ’ 


where  V  is  the  vector  (V^  V^g),  the  vector  I  is 

(I1+Iin  1'  I2+Iin  2'  X3  "  ' '  INVANES-l+Iout  2'  INVANES+Iout  1* 
and  ='  £  and  L  are  sc3uare  matrices  of  size  NVANES  x  NVANES. 


B*  Numerical  Solution.  At  the  start  of  the  calculation 
the  capacitance  matrix  is  inverted.  With  this  rearrangement, 
the  node  equations  become 


This  is  a  set  of  simultaneous  second-order  differential 
equations  soluble  in  time  steps  by  the  predictor-corrector  method 
The  initial  conditons  are 


V  =  0 


(1.4) 


and 


dV 

dt 


0 


(1.5) 


The  network  equation  Equation  1.3  requires  the  time 
derivatives  dl^/dt  for  1  =  k  =  NVANES .  These  are  evaluated 
by  extrapolating  a  third-degree  polynomial  fitted  to  the 
four  previous  values  of  current  1^  computed.  Trials  have 
been  made  for  a  35-vane  tube  with  an  ideal  "spoke-of-charge" 
model.  The  results  show  that  calculating  the  induced  current 
at  eight  steps  in  the  RF  period  gives  a  steady-state  solution 
accurate  only  to  within  30  percent,  but  with  16  such  steps 
the  error  is  reduced  to  four  percent.  To  ensure  numerical 
stability,  the  network  equations  are  solved  on  a  smaller 
time  step  (usually  1/500  the  RF  period)  than  that  used  for 
the  beam  motion. 

C.  Choice  of  the  Network  Elements.  To  relate  the 
network  parameters  to  the  measured  quantities,  the  analytical 
steps  are  as  follows: 


1.  Specific  forms  of  the  admittances  Y  ,  Y.  and  Y  are 

a  t>  c 

chosen. 

2.  The  dispersion  equation  is  derived,  giving  the  phase 
delay  and  attenuation  per  vane  as  a  function  of  the  network 
elements  and  the  drive  frequency. 

3.  The  power  flow  for  a  single  traveling  wave  is 
computed,  giving  the  characteristic  impedance. 

4.  The  driving  current  1^  representing  a  specified  RF 
drive  power  is  derived. 

5.  The  vane  geometry  is  used  to  relate  the  characteristic 
impedance  to  the  measured  interaction  impedance  for  a  given 
space  harmonic  of  the  electric  field. 

Step  3  also  yields  the  load  admittance  required  to  match 
the  cold  circuit  (i.e.,  with  no  reflected  waves  at  the  ends). 

1.  Simple  Network.  Figure  8  (Section  VIII)  shows  a 
simple  form  of  the  general  network.  It  is  specified  by  only 
three  parameters  Ga,  C&,  and  F^,  which  are  chosen  to  give  the 
measured  values  of  phase  delay  per  sectional  pitch,  interaction 
impedance  and  attenuation  per  section  at  a  given  frequency. 

The  admittances  of  the  generator  and  load  are  assumed  to 
equal  the  cold  network  admittance  in  order  to  provide  a 
match  at  both  ends.  The  directions  of  power  flow  and  wave 
motion  are  the  same  in  both  forward-wave  and  backward-wave 
tubes  when  viewed  on  the  network,  but  in  the  backward-wave 
tube  the  driving  current  is  applied  at  the  right,  at  node 
NVANES . 


.-.  .V  -  .  - 


2.  Dispersion  Equation.  Define  the  following 
parameters  for  a  signal  at  frequency  oj/2tt: 

=  cold-circuit  phase  delay  between  adjacent  vanes, 

=  exponential  voltage  decay  constant  to  account 
for  attenuation, 

=  characteristic  impedance,  v£/2P  (real) , 

=  power  flow  along  the  network  for  a  single 
traveling  wave, 

=  complex  characteristic  admittance  of  the  network 
and  admittance  of  matched  load. 

Consider  a  single  traveling  voltage  so  that  the  forward- 
wave  voltages  are 

Vk  =  Re{ exp  j  [at  -  (k-1)  j> ]  }exp  -  (k-l)a 

an  a  backward-wave  tube  has  the  voltages 

V,.  =  Re  iV.  exp  j  [^t  +  (k  -  WANES)  p  ]  }exp  (k  -  WANES )  a 

rC  i 

for  1  =  k  =  NVANES .  Such  a  wave  propagates  with  a  matched  ]1 
load  and  in  the  absence  of  driving  currents  from  the  beam. 

Using  the  general  network  equation  1.3  for  a  node  away 
from  the  ends  gives  the  dispersion  equation  as 


a 


Jchar 

P 


cos  2  cosh  x 


b 


2 


3.  Characteristic  Admittance  and  Load  Admittance, 

The  load  admittance  Y^  required  for  a  reflectionless  termina¬ 
tion  is  given  by 

Fb  sin  4>e  *  jrb(l  -  cos 


This  expression  is  the  ratio  of  the  current  (in  the  direction 
of  power  flow)  through  the  admittance  Yb  to  the  voltage  at 
vane  k.  Then,  the  characteristic  impedance  Zcbar  is  given  by 


Vki2 


‘char 


2P 


where  P  is  the  power  flow  along  the  line  at  node  k  or 


'char  Re(TjT 


Using  Equations  1.8  -  1.11,  it  is  seen  that  the  values 
of  the  elements  in  Figure  8  are  as  follows: 


2Tb  sin  <f>  sinh  a /to, 

2rb(l  -  cos  <f>  cosh  a) /to2, 
toe  /  (sin  <£Z^ba^,)  ' 

1/Z. 


Jchar ' 


fb (1  -  cos  $e  a) . 


The  total  admittance  seen  by  the  beam  at  a  single  vane  is 


Y 


2  Y.  +  Y 


For  the  lossless  case  (a=0),  this  reduces  to 


Y  =  1  =  -I 

Zchar  Z0 


for  both  forward-wave  and  backward-wave  tubes,  where  Z_  is 

3  ^ 

the  impedance  parameter  used  by  Dombrowski.  This  identification 
enables  the  two  models  to  be  compared  for  the  case  of  ideal 
spokes  of  current.  Agreement  is  excellent,  with  less  than 
0.5  percent  error  for  the  synchronous  case,  thus  verifying 
the  numerical  procedure  used  here.  Of  course,  the  present 
network  equations  are  the  more  general  of  the  two  treatments. 


The  value  of  zchar  is  related  to  an  experimental  measure¬ 
ment  in  Section  F  of  this  appendix. 


4.  Power  and  Driving  Current.  The  driving  current 

is  I.  .  cos  cot  for  the  forward-wave  tube  (at  node  1)  and 
in  1 

I.  cos  cot  for  the  backward-wave  tube  (at  node  NVANES)  .  The 

m  nv 

amplitude  is  chosen  so  that  the  time-averaged  RF  power  flowing 
into  the  input  node  in  the  absence  of  the  beam  equals  the 
specified  RF  drive  power  excluding  the  power  dissipated  in 
the  load  conductance  . 

In  the  forward-wave  tube,  the  condition  that  node  1, 
with  voltage  V^,  "see”  an  infinite  line  is  that  the  driving 
current  is 


!m  i 


v  r  v 
1 1-  l 


*  Y.  f  exn 


:  exn  'j. 


-1"  1 
x  j  j 


(I 


in  the  absence  of  the  beam. 


The  RF  power  leaving  node  1  xn  the  direction  of  propaga¬ 
tion,  towards  node  2,  is  then 


drive 


-  7  Re  '"in  1  -  YlV  v] 


where  P(3r^ve  is  the  drive  power  given  as  supplied  to  the  CFA. 
Hence  the  voltage  at  node  1,  again  in  the  absence  of  the  beam 


Vx  =  {2Pdrive/Re[e3Cp  3  <t>  exP  a  "  1)J) 


(1.16) 


'2P.  .  Z  ,  exp  a 
drive  char  * 


(1.17) 


(The  general  network  of  Figure  7  has  been  treated  similarly, 
but  requires  two  driving  currents  of  different  phase  to 
represent  the  infinite  line.)  Equations  1.14  and  1.16  thus 
give  the  amplitude  Iin  ^  for  a  given  RF  drive  power  In 

the  backward-wave  tube  the  same  current  is  applied  at  node 
NVANES  instead  of  node  1. 


D.  Oscillation 

The  driving  current  I^n  or  I^n  is  kept  constant  in 

the  presence  of  the  beam.  Then  Equation  1.15  is  used  to  define 

the  effective  RF  drive  power.  In  general  this  can  be  less  than 

the  cold  RF  drive  power  because  of  power  traveling  backwards 

towards  the  input.  The  difference  between  the  computed  value 

of  P.  .  and  the  cold  RF  drive  power  equals  the  reverse-directed 
drive 

power  that  is  reaching  the  input  node. 


If  this  difference  is  zero  or  negative,  the  power  returned 
to  the  load  at  the  input  exceeds  the  total  RF  power  that  is 
delivered  by  the  driving  current  I^n  and  the  input  power 
from  Equation  1.15  is  printed  as  a  zero  or  negative  number. 

The  interpretation  is  that  the  network  is  oscillating.  However, 
the  numerical  solution  is  stable  throughout.  Notice,  too, 
that  this  reverse-directed  wave  which  produces  oscillations 
has  the  same  cold  phase  delay  <J>  as  the  cold-circuit  wave,  while 
the  directions  of  both  the  wave  and  the  power  flow  are  reversed. 
In  the  actual  CFA  such  a  condition  shows  as  a  mode  boundary 
and  may  induce  spurious  signals  at  other  frequencies. 

The  RF  power  at  intermediate  vanes  is  computed  as  vf/zcj1ar' 
and  the  RF  output  power  is  simply  GiVNVANES2 or  for 

forward-  or  backward-wave  devices  respectively.  The  bars 
denote  averaging  of  the  instantaneous  values  over  one  RF  period. 

E.  Power  Balance 

A  power-balance  equation  is  used  to  check  the  accuracy 
of  a  steady  state  solution  as  seen  by  the  network. 

Thus,  the  input  power  from  external  drive  and  the  beam 
is  given  by 

NVANES 

_  v  _ 

TOTAL  INPUT  =  V1 (t) Iin  1  COS  wt  +  Ik(t)Vk(t)  '  * 

k=l 


for  the  forward-wave  tube ,  or 


TOTAL  INPUT  =  (t )  I-n  Ny  COS  at 


NVANES 


I.  (t)V.  (t) 


,  *.  V- 


for  the  backward-wave  tube. 

The  gross  output  power,  including  the  power  dissipated 
in  the  generator,  is 

NVANES  _  _ 

TOTAL  OUTPUT  =  G.^  (t)  2  +  ^  GaVk(t)2  +  G1V  NVANES (t) 2 

The  program  compares  the  two  expressions  and  prints  the 
difference  as  a  percentage  of  the  total  input  power. 

This  difference  is  the  rate  of  change  of  stored  energy 
in  the  network  inductances  and  capacitances,  or  equivalently 
in  the  RF  field  between  the  vanes  and  the  cathode.  It 
is  zero  in  a  steady  state  but  in  general  fluctuates 
between  ±  80  percent  of  the  total  RF  power. 

F.  Phase  Delay,  Interaction  Impedance,  and  Characteristic 
Impedance 

The  vane  RF  voltage  VRF  corresponding  to  RF  power  P  is 
given  by  (Equation  I. 10) 

VRF  =  rJ2ZchaxP 

The  corresponding  interaction  impedance  Z^nt  for  the 
fundamental  space  harmonic  is  given  by 

E  2 

7  -  RF 

int  '  232p  ' 

where  ERR  is  the  fundamental  space  harmonic  amplitude  of  the 
electric  field  on  the  circuit  and  3  is  the  corresponding 
wavenumber . 

The  wavenumber  B  is  related  to  the  vane  period  p  and  the 
phase  shift  per  vane  (as  seen  by  the  beam')  by 


2 


t.'- 


► 


8p  =  4>b  .  (1.23) 

For  the  forward-wave  circuit,  the  angles  <j>  (Equation  1.6) 
and  ^  are  both  equal  to  the  total  externally  measured  phase 
delay  from  the  input  to  the  output  divided  by  the  number 
(NVANES-1)  of  pitches  along  the  vane  circuit.  For  the 
backward-wave  circuit,  the  angle  <j)^  equals  ir-<j>,  as  Section 
G  of  this  appendix  explains. 

Neglect  the  reentrancy  of  the  structure,  assume  that  the 

total  electric  field  is  uniform  between  the  vane  tips  and 

perform  a  spatial  Fourier  analysis  to  relate  the  terms  E_p 

30  Kr 

and  VRp.  The  result  is 


Jint 

Jchar 


''sin  4>^/2  sin  acp^/2' 


4>b7T 


a<j>b/2 


(1.24) 


where  a  is  the  ratio  of  the  gap  length  between  vanes  and  the 
vane  period,  and  <f>b  is  the  phase  delay  per  pitch  as  seen  by 
the  beam.  The  interaction  impedance  is  measurable  by  perturbing 
the  fields  in  the  beam  region  with  a  dielectric  sheet.  The 
computer  program  uses  this  measured  value  as  an  input  parameter 
and  derives  the  characteristic  impedance  zcjiar  using  Equation 
1.24. 


G.  Backward  Waves 

The  network  equations  used  here  make  no  distinction 

between  forward-wave  and  backward-wave  beam  interaction. 

The  network  elements  are  chosen  so  that  the  power  flow  and 

phase  delay  from  the  RF  input  to  the  RF  output  match  the 

measured  values.  A  strapped  vane  line  is  treated  as  a  parallel 

transmission  line  with  adjacent  vanes  connected  to  alternate 

31 

straps  as  in  Figure  1.1.  Because  the  vanes  are  portions  of 
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rv 

Phase  delay  from  A  to  C  =  29 
Phase  delay  from  A  to  B  =  9-tt 


Figure  II.  Strapped  vane  line  as  a  backward- 
wave  RF  circuit.  (After  Brown) 


resonant  cavities,  the  voltages  between  alternate  vanes  (at 
nodes  A  and  C  in  Figure  1.1)  are  in  phase  with  the  wave 
propagating  along  the  line  while  the  voltage  at  node  B  has 
a  phase  tt  radians  out  of  phase  with  the  average  of  the  phases 
at  A  and  C.  In  other  words,  the  experimentally  observed 
phase  delay  is  2 <f>  along  the  line  from  A  to  C  but  the  phase  at 
node  B  is  delayed  by  <J>  -  ir  from  that  at  A.  Since  <p  is  always 
less  than  v ,  the  vane  voltages  produce  a  phase  delay  in  the 
opposite  direction,  from  C  to  B  to  A,  of  magnitude  ir  -  <p  per 
pitch.  Hence,  the  beam  interacts  with  a  wave  traveling  in  the 
opposite  direction  to  the  power  flow.  In  the  model  the  voltages 
on  the  even-numbered  vanes  are  reversed  in  sign  when  the  beam 
RF  fields  are  computed.  This  change  is  valid  only  for  the 
time-harmonic  fields;  the  time-averaged  part  of  the  vane 
voltage  (the  direct  anode-cathode  voltage)  is  not  changed 
in  sign. 

Because  the  change  in  the  sign  of  the  vane  voltage  reverses 
the  rate  of  working  of  the  voltage  on  the  electrons,  the  induced 
currents  computed  from  the  beam  must  also  be  reversed  in  sign 
on  the  even  vanes  before  they  are  used  in  the  network  equations. 
Then  spokes  of  charge  traveling  with  the  beam  are  seen  by  the 
network  as  a  wave  traveling  in  the  opposite  direction,  from  RF 
input  to  RF  output . 


H.  Comparison  of  RF  Field  Calculations 

Here  the  anode  is  treated  as  a  set  of  segments  (vanes) 
equally  spaced  around  the  interaction  region.  The  RF  voltage 
and  phase  are  given  for  vane  j  as  VRF^  and  4k.  Then  the  RF 
field  (Er,  E^)  at  position  (r,9)  at  time  t  is  given  by 

N VANES 

(r ,  9 )  =  7  V__  .G_  [  r ,  (9  -  0J]  cos(~t  -  *  )  (1.25) 


Eq (r , 0) 


NVANES 

-  2  v 


RFjGe[r'(e 


0 j ) ]  cos  (wt  -  $j),  (1.26) 


where  9  ^  is  the  angular  position  of  vane  j. 

The  terms  Gr  and  GQ  in  Equations  1.25  and  1.26  are  given  by 
Equations  H.7-H.10  of  Appendix  H. 

The  9-independent  term  of  Equation  H.7  is  omitted  throughout, 
in  consistency  with  the  approximation  that  the  sum  of  the  RF 
voltages  over  all  the  vanes  equals  zero  at  any  instant.  If  the 
modulator  voltage  is  included  with  the  vane  voltage  Vn  this  term 
must  be  present. 

The  fundamental  RF  space-harmonic  amplitude  V  is  given  from 
Equations  I.  21-1. 23  as 


(1.27) 


Vs  =  VRF  /Z 


Zint 

char 


(1.28) 


In  Figures  1.2,  1.3  and  1.4,  the  total  fields  are  compared 
with  the  analytical  values  for  the  fundamental  space  harmonic 
in  rectangular  coordinates.  The  numerical  values  correspond  to 
the  QKS1842  tube  with  an  anode  radius  of  9.708  mm  and  a  sole 
radius  of  7.748  mm.  All  the  vanes  have  the  same  voltage  amplitude 
of  287  V  (corresponding  to  an  RF  power  of  5.5  kW) .  The  voltage 
on  successive  vanes  differs  in  phase  by  133  degrees  in  the  direc¬ 
tion  of  wave  motion  as  seen  by  the  beam  in  the  backward-wave 
tube.  Equivalently,  the  phase  shift  in  the  input-output  direction 
is  47  degrees  since  the  voltages  on  alternate  vanes  have  opposite 
signs  because  of  the  strap  connections. 
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b.  Radial  Fields 

Figure  1.2  Comparison  of  exact  values  and  single  space 
harmonic  of  RF-circuit  fields  in  QKS1842  at  radius  9.2  mm. 
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a.  Angular  Fields 
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b.  Radial  Fields 

Figure  1.3  Comparison  of  exact  values  and  single  space 
harmonic  of  P.F-circuit  fields  in  QKS1842  at  radius  9.6  mm. 
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b.  Radial  Fields 


Figure  1.4  Comparison  of  exact  values  and  single  space 
harmonic  of  RF-circuit  fields  in  sever  region  of  QKS1842, 
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The  single  space  harmonic  is  a  good  approximation  at  radii 
of  9.2  mm  and  less  (Figure  1.1).  Closer  to  the  anode,  the 
exact  field  shows  more  rapid  variations  than  the  simple  cosine 
curve  (Figure  1.2).  Both  fields  are  negligible  throughout 
70  percent  of  the  sever  region  (Figure  1.3). 

Similar  close  agreement  has  been  demonstrated  to  verify— _ 

the  RF-field  calculations  in  the  SFD-261  and  QKS1319  tubes. 

Of  course  the  Green's  function  method  is  more  generally  applicable 
to  arbitrary  sets  of  vane  RF  voltages,  including  standing  waves 
of  any  wavelength  or  frequency. 


APPENDIX  J 


THE  CYLINDRICAL  BRILLOUIN  STREAM 

Estimates  of  the  circulating  current  and  the  size  of  the 

t-fltf  charge .  around  the  cathode  of  a  cutoff  magnetron  or 

distributed-emission  CFA  are  readily  obtained.  Following,  for 

32  33  34 

example.  Slater,  Osepchuk,  and  Smith,-3  assume  that  all 
electrons  move  in  circles  concentric  with  the  cathode  with 
constant  angular  velocities,  and  assume  that  the  stream  is 
in  stable  equilibrium. 

Define  the  notation  as  follows,  with  MKS  units  throughout 

e  electron  charge  (in  magnitude) 

m  electron  mass 

n  e/m 

rc  cathode  radius 

r  anode  radius 

a 

ru  radius  of  outer  edge  of  charge  hub 

n 

(in  the  absence  of  RF  fields) 

X  rH/rc 

V  anode  voltage 

B  magnetic  field 

Using  Busch's  theorem  and  conservation  of  energy,  and 
solving  Poisson's  equation  gives  the  result 


34 

This  is  quoted  from  Smith's  article,  but  with  the  correction 

33 

of  a  typographical  error.  Osepchuk  gives  the  correct 
expression. 

The  computer  program  solves  Equation  J.l  using  Newton's 
method  with  the  initial  estimate  r„  =  1/2 (r  +  r_) .  Convergence 

n  a.  C 

to  within  0.1  percent  is  attained  in  three  iterations. 

32 

The  electron  velocity  v^  at  radius  r  then  is  given  by 


v 


3 


the  potential  $  by 


ri  > 


J 


$ 


1  ,2  '  ( 

8  'B  r~ 


the  local  charge,  density  p  by 


c 

'0 


1_  d  rd  $ 
r  Jr  or  * 


and  the  circulating  current  I  by 


( J.  2) 


(J.3) 


( J.  4 ) 


-  v^dr 


( J.  5) 


where  h  is  the  width  of  the  beam  in  the  magnetic  field 
direction. 


The  last  integral  evaluates  as 


For  the  QKS1300  amplitron,  the  circulating  current 
computed  in  this  way  is  0.606  A  and  the  hub  radius  rH  is 
1.36  mm.  With  330  simulation  particles,  the  DECFA  model  with 
a  primary-emitting  cathode  and  no  angular  space-charge  field 
gives  a  circulating  current  varying  between  0.598  and  0.619  A 
as  particles  are  collected  or  emitted  at  the  cathode.  The 
computed  hub  radius  is  approximately  1.42  mm. 

For  the  QKS1842,  the  circulating  current  in  the  particle 
model  varies  between  42.2  and  51.3  A  over  1.18  cyclotron  periods, 
with  a  mean  of  47.0  A.  The  analytical  value  is  46.5  A  for 
the  same  anode  voltage  (25.5  kV)  and  magnetic  field  (0.4114  T) . 

The  computed  hub  radius  lies  between  8.30  and  8.42  mm,  and 
the  analytical  value  is  8.36  mm.  Thus  the  hub  extends  from  the 
cathode  a  fraction  of  0.31  of  the  anode-sole  distance. 

Of  the  three  tubes  shown  in  Table  J.l,  the  QKS1319  has 
the  smallest  hub  because  its  anode  voltage  is  the  least  relative 
to  the  cutoff  voltage  (at  which  the  hub  would  extend  to  the 
anode) .  The  low  ratio  of  operating  current  to  circulating  current 
in  the  QKS1842  is  due  to  a  low  RF  field. 


TABLE  J.l 

CHARGE  HOBS  AND  CIRCULATING  CURRENTS 


Anode 

Voltage/ 

Cutoff 

Voltage 


Hub  Thickness 

Anode-Cathode  Circulating 
Distance  Current  (A) 


Operating  Anode 
Current  (A) 


QKS1842 

0.57 

0.31 

46.5 

SFD-261 

0.49 

0.28 

0.0216/D1 

(normalized) 

QKS1319 

0.33 

0.23 

8.79 

APPENDIX  K 
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PROGRAM  OUTPUT  INFORMATION 

Two  programs  are  run:  the  primary  DECFA  program  and  the 
post-processing  AVERAGE  program  to  summarize  the  results. 

A  third  program,  ERF,  generates  the  RF  field  tables  but  its 
output  is  not  shown.  Appendix  L  gives  examples  of  the  output. 

The  following  information  is  printed  by  DECFA.  Much  of 
it  is  self-explanatory. 

•  Data  and  initial  conditions. 

•  At  each  time  step: 

Step  number. 

Total  rods. 

Number  of  active  rods  (NPCPRE) , 

Total  charge  units, 

Induced  current  (IDCIND) , 

Modulator  current  (IDCTOT) , 

Anode  voltage , 

Mean  radial  electric  field  around  cathode  (ERSOL) , 

Collected  charges  (DQSOLE) , 

Maximum  secondary  charges  (QNOLIM) , 

Emitted  charges  (DQSEC) , 

Charges  collected  in  emission  step  and  suppressed  (DQSAME) , 
Anode  and  sole  impact  energies  (eV) , 

Emission  charges  suppressed  as  rods  of  less  than  h  their 
initial  charge  (QSMALL) 

At  selected  time  steps: 

•  (Cathode  electric  fields) / (applied  direct  field)  (%) 

at  each  emission  site,  (charge  density) / (Brillouin  density) 
(labeled  as  " (RHO/RHOB%) " ) ,  and  numbers  of  emitted  rods. 


•  Plotted  charge  distribution  around  the  tube 

•  Currents  and  powers  averaged  over  the  preceding  RF  period 

•  Vane  voltages  and  currents  ** 

* 

* 

•  Power  balance  in  network  and  in  entire  device  (but  for 
a  single  RF  period  the  steady-state  power  balance  does 
not  hold) 

•  Averaged  results  over  one  transit  around  the  tube' 

(periods  76  to  90  in  the  example) 

•  Plots  of  the  RF  power  and  phase  distributions  over  the 
vanes  (averaged  over  1  transit  around  the  tube) 

The  post-processor  program,  AVERAGE,  summarizes  the  results 
for  all  the  RF  periods,  prints  the  means  over  1  transit  around 
the  tube,  and  prints  the  power  balance  from  these  mean  values. 

The  columns  P-ATTEN,  P- ANODE,  and  P-CATH  are  the  power  losses 
due  to  RF  attenuation  and  anode  and  cathode  (sole)  backbombard- 
ment.  C-IMP  and  EV-IMP  are  the  cathode  impact  current  and  mean 
backbombardment  energy;  ASEC  is  the  effective  secondary-emission 
coefficient  after  space-charge  limitation. 
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EXAMPLE  OF  PROGRAM  OUTPUT 


Results  of  Program  DECFA 
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APPENDIX  M 


PROFESSIONAL  PERSONNEL,  INTERACTIONS  AND  PUBLICATIONS 

The  principal  investigators  for  this  study  of  the 
distributed-emission  CFA  were  Dr.  Joseph  E.  Rowe,  Vice  President 
for  Technology,  Harris  Corporation,  and  Dr.  Donald  M.  MacGregor, 
Senior  Engineer,  Harris  SAI,  Inc. 

The  basic  data  and  much  useful  information  about  the  five 
Raytheon  CFA's  were  supplied  by  Mr.  W.  Griffin  and  Drs.  G. 
MacMaster,  L.  Nichols  and  J.  Skowron  of  Raytheon  Company  Microwave 
and  Power  Tube  Division,  Waltham,  MA. 

The  SFD-261  information  was  provided  by  Mr.  H.  McDowell 
of  Varian  Associates,  Beverly,  MA,  with  whom  many  useful 
discussions  were  held  at  Varian  on  June  22,  1978  and  subsequently 
by  telephone. 

Dr.  A.  Drobot  of  Science  Applications,  Inc.,  McLean, 

Virginia  provided  helpful  advice  about  the  modelling  of  the 
cathode  emission  and  the  modulator  in  the  DECFA. 

Harris  SAI  presented  a  paper  describing  the  results  of 
this  contract  at  the  Tri-Services  Cathode  Workshop  at  Rome 
Air  Development  Center,  April  15-17,  1980.  In  April  1980 
a  short  paper  entitled  "Computer  Simulation  of  the  Backward-Wave 
Distributed-Emission  Crossed-Field  Amplifier"  was  submitted 
for  publication  in  IEEE  Electron  Device  Letters. 
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